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ABSTRACT 

A theoretical study of a collisional-radiatlve model of 
low temperature (kT^ = 1-20 ev) heliiim plasmas has been carried 
out. The model treats Hel and Hell levels together/ includes the 
processes of d.trect ionization-excitation/ heavy x^^i^ticle colli- 
sions and the effect of optical pumping by a laser pulse. Gcanpar- 
iaons have been made between calculated population densities of 
excited states and those available in the literatxire frcan experi- 
ments of other workers. 

A comparison of the theoretical results of the present 1 
work with the results of the steady state experiments on helium 
plasmas shows that agreement is good at low electron temperatures.) 
The direct ionization-excitation process has been found to be 

12 ■ 

important in populating Hell levels at kT^ > 2 eV and n^ 10 
-3 

cm / A procedure for simplification of the elaborate model has ; 
also been worked out/ and by application of this procedure it has 

i 

been shown that with suitable neglect of even 70% of the atonic 

processes/ population densities can be predicted within 20% of the i 

! 

results obtained from elaborate models. This sirr^lified model 1 

I, 

6 12 I 

works satisfactorily in the range kT = 4-20 ev and n_ * 10 -10 i 

—3 ‘ 

cm . Application of the model used in the present work to i 

several recent pulsed laser excitation st tidies cm helium plasmas 
(at tone range gas pressures) indicate a need of lower values of 
heavy particle rate coefficients and higher values of optical 
escape factors at torr range presstores than those used in recent 



xvi 


literature. On the basis of calculations carried out on the 
method of determining electron temperature from line intensity 

9 

ratios, it has been concluded that at electron densities > 10 

>i«| ■ ' ' „ 

cm" , the line intensity ratio method is not suitable for 
determination of electron temperatxire . 



CHAPTER 1 


INTRODUCTION 

The Study of instantaneoiis population densities of 
eaoiited atoms and ions in gaseous plasmas is of considerable 
importance from the viewpoint of spectroscopic diagnosis. 
Emission of line radiation associated with cooling of plasma 
constitutes an important souirce of infortnaticai on the various 
collisional~radiative processes occxirring therein. These 
processes include electron impact excitation, deexcitation, 
ionization, recatibination, neutral-neutral collisions etc. 
Besides low-temperature laboratory plasmas, there is also much 
interest in astrophysical plasmas , fusion plasmas , afterglow 
plasmas and laser plasmas . A study of laboratory plasmas can 
provide insight into the various processes that occur in astro- 
physical plasmas . In almost all types of plasmas the study of 
the populaticMn densities of excited atomic and ionic species 
are of much importance. 

As a result of various excitation processes in a plasma 
the higher energy states of atoms and ions are populated which 
in turn depopulate by vario\iS mechanisms e.g. spontaneous 
emisslcn, coUisional transfer etc. Emissicn of radiation is 
perhaps the most easily noted plasma ptenomenon. Plasma einiss- 
Ivity can be used for identification of plasma constituents as 
well as determining pliu*m-,oonditions in laboratory as well 
as in remote astrophysical plasmas. 



The concept of complete/local thermodynamic equilibrium 

1 2 

(CTE/LTE) is of great importance in plasma spectroscopy. * In 
complete thermodynamic equilibrixim, all the properties of a 
plasma, enclosed in a black body of wall temperature ^ can 
be uniquely described by this temperature. In CTE, the kinetic 
energy distribution of all plasma pairticles is given by Maxwell- 
Boltzmann distribution function of kinetic temperature 
the Planck's law determines the distribution of the radiation 
field with a radiation temperature excited states are 

populated according to the Boltzmann distribution law with the 


excltaticai temperatxire and the atoms, molecules, ions, 

and electrons are distributed according to the mass action law 


with reaction temperature In CTE, all these temperatures 

are equal to one another, hence 


sat ^ ' ' ' ssa ' ' ass ■''*?* ' ■ ass HT*' ss' * 

^ Vail ^kin Vad ^excit /react* 

The condition of GTE is realized under very special 
conditions. The core of plasmas foimed during the short time 
Interval of strong explosions, the state of plasma in interior 
of stars are examples of plasmas in GTE. It is practically very 
difficult to realize plasmas in GTE under laboratory conditions. 
However, a large class of plasmas exists in what is known as 
local thermodynamic equilibrium (ITE) . Except Planck's radia- 
tion law, the plasmas in LTE obey all thermodynamic distribution 
laws applicable to GTE plasmas. However, the role of radiaticn 
in the establishment of local thermodynamic equilibrium is 
inportant at low electron densities. Here, the radiation field 
produced by the plaama (given by Kirchoff 's law) interacts with 



3 


the plasma constitxjents and this interaction contributes to 
•radiative transfer processes. LTE plasmas exist only when the 
following condition^ is satisfied 

Uq V 1.2 X 10^^ (1.1) 


where n^ is the electron density y the electron temperature in 
®K and E . is the excitation energy from level p to level q in 

Jr /St 

electron volts. For helixom the values of n above which LTE 

e 

will exist at electron temperatxares of 1 y 10 and 20 eV are 

respectively 1.9 x 10^^ cm”^/ 6.08 x 10^^ cm*’^ and 8.59 x 10^^ 

•*3 
cm . 


In LTE plasmas y the Saha-Eggert equation is valid and 
correlates the population density of a species having 

Charge z-lto the populaticn density of its ficoic species with 
charge z as follows^ 


n^(l)n^ 

z „e 


n, 


E 




( j) 


2g^(l) 


2iT m_M 3/2 E' . (J) - AE^ 

( l-S) e=cp(- 


h 


3cT, 


( 1 , 2 ) 

where n (1) 'and n_ (J) are respectively the population densi- 

Z ' «#• ’ ' ' ■ 

'Z— 1 

ty of the ground state of of charge z and Saha 

equilibrium population density of the jth level of species 

of charge z-iy n^ and are respectively the electron density 

and electron temperature. is the ideal icaiization 

energy of level j and A is the lowering of ideal ionization 

energy due to electric micro fields y g^Cl) and g^^j^(j) are the 

multiplicities of the levels. The Boltzmann equation correlates 

Eu (j) with n_ (1) as follows* 

'“'z-l z-1 



( 1 . 3 ) 


(j) 

2-1 


n. 


E, 


liT 




“FJ 


/ <4 \ 


Under laboratory and many native conditions# however# 
most plasmas neither follow the classical thermodynamic eqxiili- 
brixjm relations^ nor satisfy GTE conditions. They are even far 
from being in LTE and represent an open system having temperatxire 
and density gradients. The reabsorption is very weak and most 
of the photons may escape. Under sxich conditions the popula- 
ticsns of discrete levels are not described by T^ alone and cannot 
be calculated simply by applying Boltsmann distribution law. 

This class of plasmas #whiah does not obey classical thermodynamic 
equilibrium relations is termed as non-LTE plasmas. 

The properties of non-LTE plasmas can only be described 
by taking into consideratlcai the varioxas processes involving 
(a) electron impact excitation# deexcitation# ionization# 
recombination (b) photo-excitation# deexcitation# ionisation# 
and (c) diffusion etc. By considering all the processes respon- 
sible for population and depopulation of all individual levels# 

coupled, rate equations describing the rate of change of popula- 

1 

tioas of the levels with respect to time can be written in 
terms of Einstein coefficients for spontaneous emission and 
absorption# A (p#q) and B{p#q) besides various collisional excit- 
ation and deexcitation rate coefficients K‘s, The collisional 
coefficient K of particles of kind k and 1 is given by 



1 

where represents the velocity distribution functicai of 

the relative velocity v^ ^tween the particles of Icind k and 1. 

is the croas section as a function of v^. The distribu- 
tion function f (v) is obtained by solving the Boltzmann collisicai 

■ 2 ' 

equation with the proper initial and final conditions. 

In non-LTE plasmas the excited state populaticxis deviate 
from the Saha equilibrium populations due to radiation losses 
and diffusion flxjDces, The departures of the populaticn densities 
of excited states from Saha population densities are generally 
expressed in terms of reduced population coefficients (Saha 
decrrements ) p(j) and defined as 



( 1 . 5 ) 


where n( j) denotes the non-LTE population density of the jth I 

level and njd) represents the Saha equlllbrinm population density j 
Thus in non-LTE conditions the population densities of j 
excited states can be calculated only by simultaneously solving j 
the coupled rate equations obtained on the basis of a collisional- | 

f i 

radiative model. The role of electron impact processes depend on t 
the magnitudes of electron density^ electrcxi temperature and 
electron impact cross sections. Similarly cross sections of 
pShcton induced processes and the magnitude of the photon f lux in 
the plasma decide the role of photon-induced processes. Pron the 
viewpoint of ihotce induced processes, one can classify plasmas 
in two categories , When the plasma conditims are such that ail 
radiations can freely escape from the plasma without absorption, 
the plasma is called optically thin plasma. In contrast, when 



the plasma reabsorbs the escaping photons either ccsnpletely or 
partially, one has either optically thick or partially optically 
thick plasma. This absorption phenomenon is quantitatively 

incorporated in the collisional-radiative joodel through optical 

ft Q 1 0 

escape factors, *• In optically thin case, the optical escape 
factors of all radiations can be taken as vinity but in partially 
optically thick case, the optical escape factors are less than 
one, the magnitude depending on the nattire of excitation transi- , 

i 

t ion, electron density, gas pressure etc. In case of canplete ; 

I 

absorption,, the optical escape factor is taken as zero, 

, ' ' , I 

A collisional-radiative model is a versatile and convenient 
aid in understanding the nature of non-LTE plasmas. It may be 

. I" 

used for prediction of population densities of various levels vinder ! 

different plasma conditions. When applied to an experiment, it i 

11 ! 
can be used to interpret experimental data and for deriving i 

12 13 

plasma parameters. ' A collisional-radiative model simulation 

14 

for a laser Induced selective excitation experiment can also 

help in estimating cross sections of elementary processes. In a 

15 

laser generated plasma , relaxation times can also be estimated 
with the help of a collisional-radiative model, A collisional- 
radiative model along with exper^imental plasma mlssivities can be 
used for determining plasma parameters of remote astrophysical 
plasmas. 

There has been considerable amount of work on collisional- 
radiative models applied to various systems. Bates, Kingstcxi and 
McWhirter proposed a colllsicxial-radiative model for hydrogenic 
plasmas and calculated the red\iaed population coefficients p(p) 
for hydrogen atom and hydrogenic icns in optically thin plasmas. 



The processes considered include electron Impact excitation, 

deexcitation, electron impact ionization, three -body recombinaticai, 

radiative recombination and radiative deexcitation . The coupled 

differential equations describing the time development of the 

population densities of the excited states were solved nxjmerically 

under the quasi-steady state approximation. The asstimption of , 

■ . ^ 
quasi-steady state is valid because of the fact that for a wide 

range of plasmas, the Saha equilibrium population densities of 

excited levels are much less than the electron density and the 

population density of the bare nxiclei. This implies that for a 

still wider range 

n(p) < < (p 1) (1.6) I 

where n(p) is the population density of the pth level and is 
the electron density. This condition is violated only in very 
dense plasmas. Moreover/ if the mean thermal energy is much 
less than the first excitation energy, n(p) is much less than 
n(l), the ground state population density, when the steady state 
is attained. In such plasmas, the excited states attain quasi- 
steady state population densities almost instantaneously and the 
electron density and the density of bare nuclei are not signifi- 
cantly changed. The relaxation times for the excited levels are 
significantly shcarter compared to the relaxation times of the 
ground level or of the free electrons. Thus all the rate equations, 

except that for the ground state, can be set equal to zero and can 

16 

be conveniently solved numerically. Bates et al suggested that 
the number of coupled differential eqtiations could be limited in 
number by maJclng xise of the fact that for higher lying levels. 



Golllsional processes are much more important than radiative 

processes and n(p) satisfies the Saha equation. The levels with 

p greater than some value s may be grouped together such that 

p(s) is vezry close to mity and the finite nuinber of coupled 

differential equations may be solved to obtain p(p)'s. If s is 

chosen correctly and is large enough, then the line p(p) - 1 will 

be tangential to the curve obtained by plotting P(p) against p. 

17 

McWhirter and Hearn presented values of reduced popu- 
lation coefficients and relaxatiai times for hydrogenic systems 

2 7 

and used the reduced variables 9 = T/z , ti(c) = n^/z and 

n(p) *s 3c n(p)/z^^ where z is the nuclear charge and ^ is the 

ratio of electron density and the population density of singly 

18 

charged ion of charge z. Drawin solved a system of coupled 

collisional-radiativB rate equations for a hcmogeneous and steady 

state hydrogen plasma. It was found that hydrogen plasmas are 

totally optically thick towards all Lyman lines and partially 

optically thick towards ; equilibrium populations are establish 

16 -3 19 

at electron densities greater than 10 cm , Orawin also 
calculated the relaxation times for the ground and excited states 
of hydrogen atoms and hydrogen like ions. The relaxaticm times 
were found to be sensitive to the number of reabsorbed photons 
fr can the resonance lines and the resonance continuiam in collisicai 
dominated plasmas. Even under LTE conditions, in the presence of 
strcxig resonance absorption, the relaxation times of hydrogen^^^^^^^^^^^^^^^^ ^ 

were found to be of the order of 1-50 ts for n_ # 10 can"* at 

3 4 ' 

temperattures lying between 8 x 10 and 1.6 x 10 •K. 

20 

Johnson and Hinnov used a new set of transition rates 
and calculated the populaticai coefficients, ionization and 



recombination coefficients for the excited levels of hydrogen in 

S ' 1. 0 ' 

the electron density range 10 —10 cm” and in the electron temp— 

eratxjure range 2.5 x 10^-8. 2 x 10^ ®k. 

21 

nrawin and Emard investigated the role of particular 
collision and radiative processes in the final solution of coupled 
rate equations for atomic hydrogen levels. The results were'foxnd 
to be sensitive to the processes as well as number of levels 
indicated in the model. They showed that the critical quantum 
nxamber n defining the limit of partial lte could not be consi- 
dered as the critical nxiiriber of rate equations above which the [ 

population coefficients become independent of p. Drawin and J 

22 

Emard also discussed the role of cross section data and recalcu- 
lated the population coefficients, collisional-radiative reccmbi- i 
nation and ionization coefficients for hydrogen atom and i 

hydrogen like ions for optically thin and three different sets of 

23 : 

partially optically thick conditions. They also obtained the ; 

homogeneous stationary state solution and calculated the ground 

state populations of atomic hydrogen and hydrogen like ions in [ 

i 

non thermal plasmas and presented the Saha decrement for ground ■ 

' 

state population, ionization and recombination coefficients in 

8 18 -3 ’ 

the electrcai density range 10 -10 cm and electron temperature | 

3 5 ! 

range 8 x 10 °K-2*56 x 10 for optically thin and partially 

optically thick conditions. | 

The present work is on the collisicaaal-radiative model i 

of helium. Bef ore we discuss collisional-radiative models of 

helium, we might mention the collisional-radiative models of 

other polyelectronio systems. In polyelectronic systems lack of 

cross secticai and other data and large number of energy levels 



10 


make computations more difficult. Here/ we mention a few models 

applied to non hydrogenic systems other than helixim. 

24 

Park has applied a collisional -radiative model to nitro- 
gen plasmas to calculate spectral line intensities. He used 41 
groups of the bo\and states of atcamic nitrogen and calculated the 
population densities of excited states as a function of electron 

density and electron temperature following the method of Bates et 

16 

al. The objective of this work was to use intensities of two 
atcamic nitrogen lines to determine electron temperatures in nitro- 
gen plasmas. He fotond that if the equilibrium relation is used 
for determination of T_/ non eqirilibrium effects cause the appa- 

rent T to be higher than the true in an expanding plasma. 

25 

Cartwright , using a collisional-radiative model/ calcu- 
lated vibrational populations of triplet excited states of 
under auroral conditions and predicted the relative contributions 

of varioiis cascade and intersystem cascade processes. Ashraf et 

26 

al developed a collisional-radiative model to calculate the 

population densities of the triplet states of using quasi-stead; 

state conditions and Maxwellian electron energy distribution. 

27 

Glannarls and Incropera explored the effect of colli- 
sional and radiative processes on electronic state population 
densities in a cylindrically confined argon arc plasma on the 
basis of a collisional-radiative model following the method of 
Bates et al.^^ They coalesced the allowed energy levels having 
same electronic configuration and assumed levels with p > IS 
to be in equilibrium with free electrons and compared their 
calculated populatioi densities with the experimental results- 



28 29 

Van Der Sijde et al ' developed a collisional-radiative model 
of the argon ion laser system and compared the calculated popula- 
tion densities of 4 s and 4p groups with the experimental results 

of their continuous hollow cathode arc discharge experiment. t 

* 

Helium plasmas have been of much interest for a nuiriber 
of reasons. Among the polyelectronic atcans/ helium is the simplest 
one and has attracted considerable theoretical attention. Higher 

levels of helixan can be approxiirated as hydrogenic levels. It j 

1-3 i 

has two metastable states 2 S and 2 S, which are not radiatively i 

coupled to the ground state and require a collisional mechanism 

for depopulation to the ground state. The study of helium 

plasmas is also interesting dxie to the fact that helium shows 

30 

largest deviaticn from LTE amongst all neutral particles tmder 

given values of electron temperature and electron density. The 

large energy gap between the groxmd level and the first excited 

3 

State, the non -emitting nature of the 2 S state, and the presence 
of autoionizing states causing dielectronic recearibination at high 
electron temperatures and low electron densities are responsible 
for this deviation of He from LTE. Experimentally, because of 
the mcaaatemic nature of helium, it is quite single to ascertaih 
neutral concentration in plasmas. Further, collisional-radiative 
model already worked out for hydrogenic icais can be easily 
applied to He ions in a helium plasma. 

It will be appropriate to mention briefly a<.#8f«r ,pKpein±*% 
ments in helium plasmas before discussing in detail the theore- 

ticals^calculaticais that have been carried out . Hinnov and 

' '31 "32'' " 

Hirschberg ' studied the rate of disappearance of electrons 
and the oharacteristios of the afterglow spectrum in a quiescent 



helium plasma produced In the B-i Stellarator. They determined 
the electron density and electron temperature from absolute 
measurements of line intensities. They also made brief calcula- 
tions on the population densities and recombination coefficients 
in which only collisions with An = +1 were considered; collisions 
involving optically forbidden transitions were neglected. Also# 

energy sublevels were not considered separately. Motley and 
33 

Kuckes studied plasma loss in helium afterglow discharges in 

B-1 Stellarator and found that the three-body recombination 

reaction is the principal mechanism for the charge removal in i 

low temperature highly ionized helium plasma of the experiment - 

34 

Johnson compared the results of his experiment on 
afterglow and ohmic heating of discharges with results of theore- | 
tical calculations from a collisional-radiative model. The i 

, I 

12 13 1 

measxxrements were done in the electron density range of 10 -10 I 

-3 

cm and electron temperature range of 0.25-14 ev. The 

collisional-radiative model used energy levels upto n *= 25 with 

separated sublevels upto n = 8. He assiamed that for the levels 

35 

with n > 26 Saha equilibrium is valid. Johnson and Hinnov 
measured electron densities and population densities of excited 
states of neutral helium as a functicxi of time in afterglows in 

1 '1 ’ ■ "t' ' '*31 

the C Stellarator in the electron density range 10 -3 x lO'^ c^ 

and electron temperature range 0.04-1 ev. They also carried out 

calculations on popxilation densities and derived semiempdrical 

cross section formulas for transitions among excited levels by 

ccHi^Bring experimental and theoretical population densities. 

36 

Ikee and Takeyama determined n^ and T_ from measurements 
of spectral intensities of the continuum spectrum following the 



3 3 

series 2 s-n p in a brush cathode plasma in a magnetic field of 

3 

1.3 X 10 Oe and showed that for brush cathode discharge # in a 
magnetic field^n^ increases as the current density of the 

37 

plasma increases. Ikee and Takeyama also determined popula- 
tion densities of various helium levels by measxiring spectral 
line intensities emitted by a brush cathode helium discharge in 
a magnetic field. They also calculated excited state population 

densities for their experimental conditions by interpolating the 

30 

results of Ettrawin et al and compared the calc\alated and 

experimental population densities, 

38 

Otsuka, Ikee and Ishee have made studies on helium 

plasma using their TPD (Test Plasma by Direct Cxirrent Discharge) 

machine. They measured Hel and Hell line intensities at various 

points of the magnetised plasma colxjmn and also estimated plasma 

parameters at the points of emissic«i intensity meastirements , 

14 —3 

The T^ range of their experiment was ''' 0,2-8 eV/ n^ 10 cm 

' 4* 

and the gas pressure ''' 10 torr. The estimated values of n 

++ 14—3 11 

and n densities were of the order of '^10 cm and vv 10 - 

12 —3 

10 cm respectively. From emission intensities^ absolute 
population densities of Hel and Hell (1. = 4) levels were deter- 
min®d|. Ihe estimates of T^# n^^ n"^ and n''"^ were applied to 
predict, \asing a collislonal-radiative model, the expected 
population densities for i * 4 state of He|I. fkv&y also 
proposed a mechanism for the e^erved phenomena in TPD machine, 
according to which Hs reconft)lnation is the major process 
upstream in the plasma colxmai causing fiell excited states to be 
populated by gaseading of higJ^r states. At down stream condi- 
ticais, with lower electron temperatiures , recombination of He'*' 



ion constitutes the major process causing higher levels of 
neutral helium to increase in population density which results 
in enhanced emission from such states. Because of the reduced 
recombination entire plasma column was f oiand to be in quasi- 
steady state -with respect to atcmic processes , 

Hegde and Ghosh^^ carried out experiments and also made 
calculaticHis on electron beam excited helium plasmas . They 
measured emission intensity enhancements of Hel and Hell lines 
frcan a helium plasma column in a longitudinal magnetic field in 

“1 >1 O 

the range 0-700 Oe at a neutral particle density of 10"^ cm” . 

The measxired Hel lines were n^D-2^P (n = 3-5 )/ n^S-2^P (n « 3-5), 

n^P-2^S (n s 3, 4), n^D-2^P (n = 3-6), n^P-2^S (n =* 3, 4), 

11 

n P-2 P (n « 4, 5) and that of Hell was i « 4 toi ■=» 3. They 

also measured electron density and electron temperature as a 

function of magnetic field and f omd that on increasing the 

field from 0 to 700 Oe, the axial electron density increases from 
10 ’“‘S 12 '""3 

10 cm to 10 cm and the electron temperature falls from 

11 eV to 6 eV. The experimental emission enhancement profiles 
were compared with those calculated on the basis of a collisional- 
radiative model. 

Several laser indiaced selective excitation studies have 

39 

also been made in helium plasmas. Catherinot et al optically 
1 111 

pumped 2 S-3 P and 2 P-4 D transitions of Hel with a txinable dye 

laser of power 500 W/pulse at 5000 A with a pulse width 5 ns and 

spectral width 0,1 A in a helium glow discharge. The measurement 

of temporal decay of the fluorescence revealed a large excitation 

1 3 

transfer between the n D and n D levels. They estimated a cross 


. -14 2 1 

section o = 5 x 10 cm for the reaction He (4 D) + 

He (l^S) He (4^D) + He (l^S). 

40 

Dubreuil and Catherinot used a two laser absorption- 
perturbation technique to study the poptilation relaxation of 2^P 

level in a helium glow discharge. They pumped the transition 

3 3 3 

4 D-2 and studied the population relaxation of 2 P by a probe 

3 3 

laser beam tuned to 3 d- 2 P transition. Using the laser pertiir- 

41 

bat ion method^ Dubreuil and Catherinot determined the optical 

11 11 
escape factor for 3 P-1 S transition. They p'umped the 2 s-3 P 

1 

transition with a tunable dye laser and analysed the 3 P level 

population relaxation for different experimental conditions in 

a helium capillary glow discharge . The optical escape factor 
11 

for 3 P-1 S transition at 350 "K and 0.4 torr pressure was 
estimated to be 0.014. 

42 11 

Yasumaru et al optically pumped 2 S-3 P transition 

of Hel with a tunable dye laser ^ PWHM 5 ns^ spectral width 0.1 A 

11 

and monitored the population densities of 3 P and 3 D levels at 

two presstires 0,4 torr and 5 torr^ and estimated cross sections 

for the excitation transfer processes, 

43 

Thoma measured the absolute helium emlssicn coeffi- 
cients in the range X « 109-540 nm from the radiations originating 
from axis of a cylindrical helitjm arc at 1 atmospheric pres sxire 

with T„ in the range 2.5 x 10^ ®K-2.65 x 10^ ®K and n_ in the 

"IS ■ ■ 16 ■ ■ ■ ■ 

range 3 x 10 -4 x 10 cm « This measurement was claimed to be 

helpful in the use of hellwn continuian radiation for calibration 
purposes. 

The above examples provide an idea of the experiments 
done in hell\an plasnias. Besides the theoretical calculations 



described in scxne of the above studies, there has been consider- 
able work solely on collisional*^a.diative models, we describe 

here some representative ones. 

44 

Bates and Kingston calculated radiative reccmbination 

coefficients, quasi -equilibrium electron temperature and ion 

temperatxare using an approximate expression for ion-electron 

45 

elastic scattering cross section. Bates, Bell and Kingston 

studied theoretically the decayjLng properties of an optically 

thick helixjm plasma in the electron temperatxare r^ge 0,09-1.4 ev 

and electron density range 10 -10 ■ cm” . Considering Penning 

ionization, they found that the collisional-radiative reccmbi- 

nation is a significant process in populating the first 

excited state. ' 

Drawin and co-workers, using their collisional-radiative 

model, carried out extensive calculations on helium plasmas and 

presented nxjimerical values of population coefficients in a wide 

range of plasma parameters, in their collisional-radiative 
30 46 

model ' they separated the level system of helium into two 

parts, the singlet system termed as the X system and the triplet 

system ternrad as the Y system and the two systems were separately 

taken into accoxmt. The separated sublevels were used only up 

to the principal quantum number n = 2, and for levels n ^ 3 tte 

singlet and triplet stiblevels of a particular principal quanttam 

number were grouped separately. The Jiighest level taken into 

account was that of principal quantum number n « 25. Thus a 

total of 51 levels were tised in the collisicnal -radiative model. 

' 47 

In an earlier work , Drawin Included the atom-atom exchange 

' 30 

collisions between the singlet and triplet levels but later on 


# 



heavy particle collisions were dropped in order to get a linear 
system of equations. To compensate the error ^ they incorporated 
electron atom exchange collisions between all levels of triplet 
and singlet systems. The rate coefficients were calculated 

I 

using a Maxwellian velocity distribution. The 51 linear simul- 
taneous equations thus obtained were solved nimeric ally both for 
hcanogeneous stationary state and inhomogeneous transient state 
solutions. The numerical values of collisional-radiative volume j 
recombination and ionization coefficients for quasi-stationary 
helixim plasmas in n^ range 10^-10^^ cm”^ and T^ range 0,25 x 
10^ °K-2.56 X 10^ "K were presented. They^*^ also calculated 
Saha decrements of ground state population densities in the 
n_ range lO^^-lO^^ cm**^ and T^ range 8.0 x 10^-1.28 x 10^ ®K 
for optically thin, slightly optically thick and strongly 

absorbing plasmas. | 

48 ! 

Drawin and Emard listed extensive tables of reduced j 

population coefficients for various levels of Hei in a wide range | 

of plasma parameters using six sets of optical escape factors . j 

22.49 ! 

For hydrogen like ions also^ they calculated the values of | 

( 

population coefficients for optically thin and paartially opti- 

cally thick conditions in reduced electron density « n^/z ) 

4 IS ^3 ' ' ' 

range 10 -10 cm and reduced electron temperature “ 

' '2 ' ^ 3 ' S 

tVz ) range l x 10 -2,56 x 10 ®K, In all the above calcula- 

e ; 

t ions, however# the separated stiblevels are used cxily up to n » 2 
ajxl the Hel and Hell levels are treated separately.. Mostly 
empirical or semiempirical cross sections were xased. The term 
system of energy levels is largely simplified and applies to 
plasmas of siofficiently high electron density. Moreover the 



model does not take into account the processes like ionizaticai 

due to collisions of two metastable states^ 

1 1 SO 

Hegde and Ghosh * developed a collisional -radiative 

model which considers both Hel and Hell levels simultaneously. 

They used completely separated stiblevels up to n = 5 ^ with 

sublevels of principal quantum ntannbers 6-12 grouped together. 

In addition to these levelS/13 levels of Hell were also Included. 

They interpreted the results of their emission enhancement 

experiments on the basis of this model. They also calculated 

relaxation times for establishing steady state population in 

optically thin helium plasmas. However mostly semiempirlcal 

analytic forms of cross sectiais weire used in this work. 

51 

Hess and Burrell in their collislonal-radiative model 
considered completely separated sublevels up to n = 4 and 
included the effect of oftical pumping. They presented results 
for both steady state and quasi-steady state solutions . The 

4 

enhancement of populatiohs on pumping with a laser of power 10 
2 «> 

W/cm -A having various rise times ^ were also presented. The 

enhancement was found to be larger at low electron densities. 

In their model, Hess and Btirrell did not include Hell levels. 

Neutral-neutral collisions and dielectronic reccanbinatlons were 

also not Included. 

52 

Pujimoto applied his collisional-radiative model to a 
discharge plasma. This model takes into account completely 
separated sublevela up to n = 7 with an upper limit of n ** 26 and 
deals with time-independent solitt;iariS, Two types of formulations 
were emphasized in this work. In the first case, the population 
densities of the levels located abcwe a sufficiently high lying 



level V were assumed to be given by Saha-Boltzmann equilibrium. 

Puirther# the rate equations of the levels p < v, except fpr 

groiind state and two metastable states were set equal to zero. 

The rate equations for the ground state l^S, and the two meta-* 

stable states 2^s and 2^S were solved independently. In the 

second case, when the change in population densities is slow 

enough as compared to the relaxation times, the rate equations 

of all the level p v, except that of the groiind state, were ! 

made equal to zero and population densities of the excited states 

were calculated in terms of ground state population density. The 

4 

results were presented at T = 1.6 x 10 and in the n range 

0 0 

•5 2.8 **3 

10 -10 cm . The popiolation density distribution among excited 
levels was also interpreted in terms of excitation-ionizaticai | 

mechanism in the plasma. This model used mostly semiempirical | 

cross sections and did not Include Hell levels along with the 1 

53 1 

Hel levels. Pujimoto also discussed the Icinetics of Icaiizationi* I 

' ' ' ' ! 

recombination of a plasma. He showed that the populaticai density ! 

of ions in a low density plasma could be described by corona | 

■ i 

eqxjiillhrium or by the capture cascade scheme. In case of a | 

! 

Optically thln^high temperature plasma for hydrogen like ions 

the population density per unit statistical weight of a given 

level was foxind to be proportional to inverse of the square root 

of the principal quantum number of that level. 

14 

Dubrexiil and Catherinot studied quenching and excitation 
transfer of n * 3 helium sublevels in a low pressure 
discharge by means of a laser perturbaticMi method. They calcu- 
lated the pert-urbations in population densities induced by a 
laser by solving the coupled rate equations as a function of time* 



By comparing the theoretical resxilts with the experimental ones, 
they deduced radiative coefficients and collisional cross sections.; 
At a pressure p < 10 torr and discharge cxarrent i < 50 mA, they 
found that only radiative and atcan-atom collision processes 
contribute to quenching and excitation transfer in the n = 3 ‘ 

levels of Hel. Also, by using this technique, they^"^ estimated a 
set of thermally averaged cross sections for quenching and 

excitation transfer processes for n = 4 level of Hel. 

55 ■ 

Dothan and Kagan proposed a scheme for calculating 

populations of excited levels and intensities of spectral lines in ; 

the positive column of helium discharges at intermediate jaressxire 

(p > 5 torr). They used a non-Maxwellian electron distribution ^ 

(for Eg > 19.8 eV) and found the, -excitation from 2^S level to be 

very Important for popxilating higher singlet levels. i 

The collisional-radiative models mentioned above differ 

in the number of energy levels used, choice of the processes 

included, and the nature of atomic parameters employed. All the 

collisional-radiative models cited here, except that of Ref, 11, 

do not consider Hel and Hell levels simultaneously. Whereas in 

isolated ranges several of these models do predict populatican 

densities satisfactorily, due to their peculiar choice of atomic 

processes and parameters, none is applicable over a wide range of 

electron density and electron temperat\ire , Also there has been 

no concerted effort on inclusicxi of experimental cross sections, 

and this could be one reascai wl^ agreement between theoretical and 

experimental populaticai densities is unsatisfactory in several 

cases. Further, in the above stxjdies, there has been no attempt 

to explore the quantitative role ©# a particular rate process in 



21 


populating a given level. There seems to be thus a need for an 

improved model which eliminates some of the above deficiencies. 

The purpose of the present work is to evolve an improved 

collisional-radiative model of helium plasmas. This is carried 

56 

out by (a) first incorporating in the model, as many as possible 
experimentally obtained excitation cross sections and then examin- 
ing predictions from the model against populatioi densities 

observed in steady-state experiments (a detailed examination of \ 

57 

the role of direct ionization-excitation process is also made ) , 

^ , I 

(b) quantitatively examining the role of a given process in i 

. . I: 

populating a particular level and thus evolving a detailed popu- : 

58 

lation mechanism, and (c) critically examining the results of 

a n\amber of recent laser induced selective excitation experi- i 

ments. On the basis of theoretical results, a brief examination I 
■ ■ ' [' 

is also made of the line intensity ratio method for estimating I 

t: 

I 

electron temperatures. i 

We first present methods of calculations and choice of j 

atomic parameters. This is followed by results of the collisional-| 
radiative model calculations, and then finally a discussion of all j 

j 

the results. . ' J 



CHAPTER 2 


THE MODEL, . ATOMIC PARAMETERS AND METHOD OF CALCULATION 

In this chapter, we first present the collisional- 
radiative model used in the present work, then details about the 
cross sections employed^ and finally present the procedwe adopted 
in the calculations. 

! 

2.1. THE COLLISIONAL-RADIATIVE MODEL 1 

I 

The various collisional-radiative models suggested for 
the helium system differ in the nxomber of atomic levels used, 5 

mode of grouping of sxiblevels, choice of atomic parameters and 

methods of solving the linear simultaneous equations . The | 

56 i, 

collisional-radiative model used in this work consists of i 

completely separated Hel sublevels up to principal quantum nximber ; 

n ** 5. The sublevels of principal quantum numbers n « 6-12 are | 

grouped together and are treated as single levels . This makes a j 

total of 32 sublevels for Hel. Besides these levels^lB levels of | 

■ ■ ■ ' ' f 

Hell have been incliKied, thus making a total of 47 levels. The j 

-1 59 ' . I 

energies (in cm ) and the weight factors for these levels are ; 

given in Table 2.1. In this section the levels of Hel are repre- i 

i 

sented by p and q while the levels of Hell are represented by . i | 

and J. ^ • I 

■ '■ ' ■ '' ' ' ' ' ' , , ' i' 

The following processes are included in the model. Here 

K's represent electron impact excitatien and deexcitaticHi rate 
coefficients, A's represent spcajtaneous transition jarobabilities 
and B's represent radiative recanbination rate coefficients. 



Table 

2.1. 

Energy levels of 

Hel and 

Hell 

used 

, in the model 

bevel 

No, 

(P) 

I' 

State 

1 

; Energy 
(cm~^) 

i 

1 s 

1 1 

! ! 

! P i 

1 1 

1 I 

1 f 

1 1 

1 1 

i 1 

Level 

No. 

(p) 

"T 

1 

J State 

1 

1 

1 

Energy 

(cm"^) 

E 

P 

T — — — 

1 

1 

l^s 

0 

1 

17 

■ 4^F^ 

Ar. 

4 F 

191452 

28 

2 

3 

2 S 

159856 

3 

18 

4^P 


191493 

3 

3 

2^S 

166278 

1 

19 

5^S 


193347 

3 

4 

2^P 

169078 

9 

20 

5^S 


193663 

1 

5 

1 

2 P 

171135 

3 

21 

3 

s-^p 


193801 

9 

6 

3 

3"’S 

183237 

3 

22 

3 

5-^0 


193917 

15 

7 

3^S 

184865 

1 

23 

5^D 


193919 

5 

8 

3^P 

185565 

9 

24 

5^F, 

5^F, 

193921 

64 






5^Gv 

5^G 



9 

3^D 

186102 

15 

25 

5^P 


193943 

3 

10 

..1 

3D 

186105 

5 

26 

6(n) 


195251 

144 

11 

3^P 

186210 

3 

27 

7 


196070 

196 

12 

4^S 

190298 

3 

28 

8 


196595 

256 

13 

1 

4 S 

190940 

1 

29 

9 


196954 

324 

14 

4 P 

191217 

9 

30 

10 


197213 

400 

15 

3 

4 D 

191445 

15 

31 

11 


197398 

484 

.16 

4^D 

191447 

5 

32 

12 


197543 

576 






TL 

OD 

198311 



mm ■»m mm m-m ■mm 
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The energy levels of Hell used in the model for i = 1 to 15 were 
calculated from the expression 

» 109677 X 4 (1 - ■^) + 198311 cm*"^^ ** 2i^ . 




( 1 ) 


Electron impact excitation and deexcitation : 

He(p) + e He(q) + e (2,1) 

(ii) Electron impact ionization: 

He(p) + e > He'^(i) + e + e (2,2) 

This process represents the direct ionization-excitation 
of Hel levels (p = 1-32) to excited Hell levels (i = 1—15) 

(iii) Three -body recombination: 

Hfe’*'(l)+e +e ■> He (p) + e (2.3) 

Since recombination of higher excited states of Hell is 
not important in populating Hel levels, only the recombi- 
nation from the ground state of Hell has been taken into 
account; K(i,p) represents the rate coefficient for this 
process, 

(iv) Spontaneous radiative transitions: 

He(p) » He(q) + hv (2,4) 

(v) Radiative recombination: 

He"*^(l) + e > He(p) + hv (2.5) 

(vi) Heavy particle collisicaas involving excitation, deexci- 
tation, ionization and reccMtibination I 

He(p) t He(l) He(q) + He(l) ' (2.6) 

KN(q.p) 

He(p) + He(l) He’^(l) + HeCl) + e (2.7) 

KN(i.^p)i„l 



At high gas pressures (^1 torr) the neutral-neutral 
collisions become significant. 

(vii) For Hell/ the processes similar to processes (2,1)/ (2,2)/ 
(2,3)/ (2,4) and (2.5) for Hel, are included. They are 
represented by the following equations: 


He''^(i) 

K(lr jX ^ 

^ ® JcTj'/i) ^ ® 

(2.8) 

He'^(i) 

+ e He'^’*' + e + e 

(2.9) 

He^"^ + 

e + e He'‘'(i) + e 

(2.10) 

He'^(i) 

i> He'^(j) + hv 

(2.11) 

He'^^ + 

e -» He’^(i) + hv 

(2,12) 


In the above equations/ K(i,/j) represents the rate coeffi- 
cient for the electron impact excitation from ith level to 
jth level and K(j/i) is the rate coefficient for the 
reverse process, K(i/C) denotes the rate coefficient for 
the electron impact ionization from ith level of Hell to 
the continuum, c represents the continuum of free elec- 
trons. K(C/i) is the rate coefficient for the three-body 
recombination . A ( i / j ) and 6 ( i ) respectively represent the 
spontaneous transition probability and radiative recombin- 
ation coefficient for Hell, 

In partially optically thick plasmas / the reabsorption of 

radiation is taken into account through optical escape 

83 ' 

factors and A.,. In such a situation A (p/Zj) and 

A(i/j) in Eqs. (2,4) and (2,11) are replaced by x 

■ ,P^ 

A(p/q) and A^j x A(i/J) respectively/ where A^^ and 



26 


are optical escape factors for the corresponding transi- 
tions of Hel and Hell. 

(viii) Optical pumping: If a laser beam of radiation density 

p(v) (energy per unit volume per \anit frequency interval) 
is tuned to a particular transition p q^ the population 
of upper level q could be significantly enhanced. This 
laser induced selective excitation may be described as 
follows: 

He(p) + hv He(q) (2.13) 


where ^(p/q) and W(q/p) are the rates of induced absorption 
and emission given by W(p,q) = p(v) x B(pyq) and W(q/p) = 
p(v) B(q^p) where B(p/q) and B(q/p) denote the Einstein 
coefficients defined by the following equations (see, 
however Section 2.3,5) 

_3 

B('i/^i?) = — g- . A(%,^) (2.14) 

8 ir h V 

^ . B(^,1P) (2.15) 

where c is the velocity of light and h is the Planck's 
constant . v represents the frequency of the transition 
p q and g^ and q^ represent the weight factors of levels 
p and q respectively. 

In the model, autoionization processes are not taken into 
account and the losses* due to diffusion to the walls for 
all the excited states of Hel and Hell are also assumed 
to be negligible. 



27 


2.2. ATOMIC DATA 
2.2.1. Cross Sections 

While choosing cross section data/ emphasis has been given 
on the use of experimental cross sections. When more than one 
experimental cross section Is available, the one which Is more 
recent has been used. Theoretical cross sections have been used 
only when experimental cross sections are not available. For a 
given transition. If more than one theoretical cross section are 

I 

' ' ' ■ 

available In the literature, the one which leads to calculated * 

* 

population densities closer to experimental results is used. 

(a) optically allowed transitions: Among the optically i 
allowed transitions, the experimental cross sections of Westerveld 

ot al have been used for the cases 1^S-2^P and 1^S-3^P. For ' 

11 * 
1 S-4 P, the data have been taken from the experimental work of [ 

61 

Donaldson et al. The cross sections used for the transition i 

1.^1 62 * 
1 S-5 P are from Scott and McDowell, Theoretical values of the 

. ■ I 

63 

cross section of Flannery and McCannS have been used for the | 

transition 2^S-3^P. For the transitions 2^S-2^P, 2^S-3^P, 2^S-4^P, | 

3 3 ' ! 

and 2 S-5 P, the data have been taken frcm the theoretical calcu- 

64 ! 

lations of Ton-That et al. The cross sections for the remaining 

65 

optically allowed transitions have been taken from Drawin. 

Following is the expression used by Drawin for electron impact 

' ' , . '■ ' ■" i 

excitaticai cross section „(E_ ) I 

e” 2 ■ 

9<Vq> 

Where E represents the electrcai kinetic energy, a_ is the Bokr 
radixis, E^ is the ionization potential of hydrogen equal to 



zo 


13.59 eV, is the difference in the energies of levels q and 

P« f(P/g) is the absorption oscillator strength. g(U ) is given 

P/<3 

by ' 

U - 1 

g(U ) = a (-£^2 ) ln(1.25 6 u^. ) (2.17) 

pq P/q p.q 

P/q 

Here U _ = = ; a and 3 „ are constants for a given tran- 

P^q E P#q P/q 

sition. In general o _ is about 1 and 3 lies between the ; 

P/q P/q 

values 1 and 3, The oscillator strength values have been taken 

fifi 

from the tables of Wiese/ Smith and Glennon, Rate coefficients 
are obtained by integrating the cross section expression over 
Maxwellian electron energy distribution using Bq. (1,4) ,. 

(b) Optically forbidden transitions without change in ‘ 

11 i 

multiplicity; The cross sections for the transitions 1 S-2 S* 
1111 

1 S--4 S and 1 S-5 S are taken from the experimental results of 

67 6 A 6Q 

Rice et al . Raan et al and Pochat et al . respectively. 

11 11 I 

The cross sections for the transitions 1 s-3 D and 1 S-4 D have i 

been taken/ respectively, from the expexo-mental results of St. John 

70 71 1 

et al and Moussa et al. The theoretical values of Flannery 

■ 63 ' 'll' 1 1 '■ ^ 

and McCann have been used for the transitions 2 S-3 S and 2 S-3 D. 

72 

Theoretical cross sections of Khayrallah et al have been taken 

for the transition 2^S-3^S. For the transitions 2^S-3^D/ 2^S-4^S/ 

2^S-4^D/ 2^S-4^F/ 2^S-5^S. 2^S-5^D/ 2^S-5^P/ and 2^S-5^G/ the 

64 

theoretical cross secticais of Ton-That et al have been used. 

For the cross sections of other optically forbidden transitions/ 
without change in multiplicity/ that are not available in recent 

• 65 : ^ 

literature, following expression due to Erawin has been used. 





4'^ a? Q (1 - u“^ ) 

o P/q p^q p^q 


(2.18) 


where U „ has the same meaning as in Eq. (2.17). Q is a 
constant for a given transition and its values for various transi- 
tions in helium are tabulated by Earawin.^^ 

(c) Optically forbidden transitions with change in 

multiplicity; Among these transitions, the experimental cross 

73 

sections of Jobe and St. John have been used for the transition 
1 ■ • 3 13 

1 S-2 P. For the transition 1 s-2 s, the theoretical cross sections 

74 

of Lins de Barros and Brandi have been used. For other such 

65 

transitions, empirical cross sections from Drawin have been 

taken.. For such transitions which involve the ground state, 

65 

Drawin has used the following expression 


q ^(E^) = 4 ir at Q 

P/q e o p,q 


(2,19) 


where n « 5 and the other symbols retain their meanings menti<xied 

earlier. For transitions between excited states i.e. not involving 

3 X 3 X 65 

the ground state (e.g. 2'^S-2'^S, 2'^S-3‘^S etc.), Drawin® suggested 
the following expression for the cross sections 


q (E ) a 4 '»r a^ Q u”^ 
^p,q^ er o '^p,q p,q 


( 2 . 20 ) 


The values of Q_ are given by Drawin. The cross 

Jbr Si 

sections for the optically forbidden transiticais with change in 

65 

multiplicity not given by Drawin are taken frcati Moisewitsch and 
75 ^ ^ 

Smith. 

(d) Electron Impact ionization; For the cross sections 


of electron Impact ionization processes involving the excited 



states of neutral helixim and ground state of helitira ion, following 
expression due to Gryzinski"^^ is used. 

£i 

p 

where E* is the ionization potential of the pth level and E is 

Jr' 0 

the electron kinetic energy. Putting E /E' = x, g (x) is given by 

0 

^p^^^ ^ X Cl + |(1 -1^) ln(2.7 + (x-1)^/^] 

( 2 . 22 ) 

(e) Direct ionization-excitation: For the direct 

ionization-excitation from the ground and excited states of Hel 

to the excited states of Hell/ following empirical formula for 

57 

the cross section has been used: 

E^ 2 

q .(x) = 46 X 10*^ ira^ [=X-] ^ ln(10 x) (2.23) 

^ ° -^pi x^ 

H ^e 

Where E^ » 13,59 ev, x = = — / E j = E.-E and E. and E respec- 

X px X p X P 

tively/ are the energies of ith level of Hell and pth level of 
Hel, The cross sections given by this empirical formula agree 

77 

well with the experimental cross sections of Weaver and Hughes 

78 

and Sutton and Kay for level i = 4 of Ifell. A comparison 
between the empirical and experimental cross sections is given in 
Figure 2.1. 

(f) Heavy particle collisions: For neutral-Hrieutral 

excitation/ deexcitation/ iaiization and recombination processes 

79' 

the expressicais for rate coefficients given by Drawin and Emard 
have been used. For neutral-neutral excitation/ the rate coeffi- 
cient KN(p/q) is given by 
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KN(p^q) = 


2 

128 IT a? £(p,q)(~^) 


kT 




O'E 


P/q 


1/2 


4* 




(2.24) 




where a - E- , E and £(p,q) retain the same meaning as in Eq. 

U X P/H 

(2.16). Tg represents the gas temperatxare., m^^^ and m^, respec 

tively stand for the mass of helixim atari and mass of hydrogen 

atom, W represents the excitation energy in vinits of kT^ and 
P/q g a 


is given by (W_ _ ) is defined as 

g 

2 


p.q kT_ - m^^-p.q- 


1 + 


V (w ) = 

"Vie 


w 


2x3^ 


2 m, 


1 + ( 




%e ^ "'e^^p.q^ 


(2.25) 


The recombination rate coefficient KN(q,p) is given by 


KN(q^p) « ^ . KN(p,q) exp(W ) 

yp P/q 

where g^ and g^ are the weight factors for levels p and q 
respectively. 

The rate coefficient KN(p^i) for neutral-neutral ioniza- 
tion is given by 




kT„ 1/2 


KN(p,i) = 128 tr aQ(^) 

o Sjj, P rritij^Q 




“--"lie 


(Wi> 


e P 


(2.26) 


Where W is the ionization energy of level p in units of kT and 
^ E ' ^ 


is given by = k^ ' ^p nuniber of equivalent electrcns 

in shell p, is given by Eq. (2.25) with W^',, replaced by 

P/q 

Tl® rate coefficient KNU,p) for the reverse process Is 
given by 





KNd^p) 


g 


^ itL.^ + 1/2 , ; 

£_ (Jfe h 

' wi nn • 


2g (2 tr 


He 


T ' ”.T^ 

g e 


KN(pd) exp(-U + 2(U 

P P P 


( 2 . 27 ) 


e 


(g) Hell transitions: Rate coefficients, for electron 
impact excitation and ionization for Hell are calculated by using 
the expressions given by Bates et al.^^ For electron impact 
excitation/ the rate coefficient K(i/j) is given by 


K(i/j) = 4.75 X 10 


-5 _(- 

9 1/9 ^P'^- ) 




z2 tV2 


(2.28) 

Where z is the nuclear charge and E. . is the excitation energy 

X/J 

i.e. “ Ej-E^. The rate coefficient for electron impact 

ionization K(i/C) is given by 

K(i/C) = exp(- 5 ^) (2.29) 

■*■6 

where E^ is the ionization energy of level i. 

2,2.2. Other Data 

Spontaneous transition probabilities for Hel are taken 

AO 

from Wiese et al and those for Hell from Green et al.°^ 

Radiative recombination coefficients 8(p) for Hel are obtained 

81 ' 

by Integrating the following expression due to Kramers 


100 „4 ^10 

cr(v) B 

qiv; 2 . 2 3 

c h V n 
e 


g(v) 


( 2 , 30 ) 



12 

where v is given byhv=hv^+^m^v , v is the frequency of 
radiation^ the threshold frequency and g(v) is the free-boxand 
Gaunt factor, z is the nuclear charge. 

QO 

3(i) for Hell are calculated from the following Seaton's 
formula for hydrogenic ions 

0(i) = 5.197x10"^^ z x-Y^ S^(r>) (2.31) i 

where D s= 157890 z^/T^ and x, = . S (D) values are obtained by 

Q X n ■ 

interpolating Seaton's values. i 


2 •2.3* Optical Escape Factors 

Optical escape factors have been calculated using the 

83 

treatment of Drawin and Emard. They considered stark broadening 
and used a Voigt profile to calculate the optical escape factors. 
The expression for the escape factor for ifel and Hell is as 
follows; 


1 + 


ii. 


2 + 

AV « ' 

1/J 1 + [trlnd + 


X 


J72 


-y(0. 

“ij 

1 

A 


HP 

2 + X 


1 + 6 1 + ( 


M— 

T72 


(2.32) 


where 



35 


f . s 4.11 X 10”^ ... j- J X X 

±j _l/2 ±j 

g 


“I7277T TT 

IT (1 + 


.(1 - 


, )] 


(1 + ir a t . ) 


ij 


and 


3.72 X 10"^^ j^(j-l) 

*ij - — - - ^ 

ij g 


In the above expressions , ^ = 100, n(l) is the ground 
state population density, 1 is the diameter of the plasma t\ibe, 
f (i,j) represents the oscillator strength, the wavelength and 
T denotes the gas temperature (or the ion temperature in case of 

•I , ' , 

Hell). 


2,3. METHOD OF CALCULATION 


2,3,1. Rate Coefficients 

As mentioned above, rate coefficients are obtained by 
integrating the cross section over Maxwellian velocity distri- 
bution of electrons. Two cases may arise — cross sectiai is 
given either in the form of an algebraic expression or in the 
form of data (as in the experimental cross section case). Methods 
of calculations in both the cases are described below. 

(a) When an empirical facmula for cross secticai is given; 
For example, in ease of optically allowed transitions, the empiri- 
cal expression given by Drawin^^ is as follows 


^P/g e 


4w a: 


eH 2 

(^~^) f (p,q) g(U cj) 
p^g 



where g(U ) 
^ P/q 


U 


= a 


£2- 

P/q 


ln(1.25 3^ „ 

P/q P/q 


pq 


On integrating the above expression over Maxwellian electron 

it it 

energy distribution ^ in the range of velocity v to oo, where v 
is the threshold velocity, we get the rate coefficient K(p,q) as 
follows: 


K(p/q) 


OD 


E- 2 


f 4'n’ a ( y ■ — ) f(p,q) g(u ) f (v),v.dv 

* P/q 


(2.33) 


and 


f (v) 

1, 1/2 
= <|) 

it 

V « 

2E 1/2 

( £*2) 

m 


e-V2 n-"ATe ,2 ^ 

■^e 


« k is the Boltzmann constant and T is 

e 


the electrcai temperatvire , m and v respectively are the mass and 
velocity Of the electron. ie the difference in the energies 

of levels q and p. 

„ 2E dE 

Putting -r mv = E , which gives v dv = 
z © 

Eq. (2.33) can be rewritten as 


m 


K(p.q) = 4'Ta^(E“)2 2- J® _|_f(p.q) ftW ) 


-E^/kT^ 

e ® ® E_ dE 

e e 


(2.34) 


E 


Sxibstituting 

gets 




E 


P/q 


K(p,q) =C* a 


U and putting the value of g(u _) , one 


P/q 


P/q 


P/q i 


r u - 1 

— In (1.25 

" P/q P/q P/q 


."^P/q’^P/q^^ 


dU 


P/q 


(2.35) 



where C * 




4^a2(E^) (|) (j^) 


m 


Substituting U -1 = Y/ and 


E 




- kT, 


= a, Eq, (2.35) reduces to 


00 


K(p/q) = ^*‘%q £(P/«l)e"“ / In {1.25 e_ _(1+Y)} e““^ dY 


p^q 


(2.3&) 


further substituting oY = one gets 


-a 00 


K(p.q) = C f(p,q) V J TTz 

(2.37) 

Sxibstituting the numerical valxje of c',Eq, (2.37) is reduced to 

“ ^ f (p#q) e””® 


CD 


K(p^q) « 5.436272598 / f(Z) e“^‘ dZ (2.38) 

m / Q 


a 


This equation can be integrated n-umerically using Gaussian 
quadrature method involving Ldguerre polynomials 
00 .7 n 

/ f(Z) e"^ dB = I (2.39) 

o i=l 

where H^'s and a^'s are respectively the weight factors and zeros 

84 

of nth order Lagxaerre polynomial. The accuracy of the results 
depends upon the proper choice of the order n of the polynomial. 

(b) In case of experimental cross sections, when the 
numerical values of cross secticai are available at discrete values 
of E_# following procedure for integration is adopted. 

If Q(E^) is the cross secticm as a function of E^|K(p/q) 


is given by 



K(p,q) 


(2.40) 


OD 

= / Q(E^) f(v) v.dv 

V 

where f (v) and v retain their old meaning. Siibstituting for f (v) 
12 

and putting mv = the above equation reduces to 

K(p.q) = (~) .%/ Q(E^) e ® ®E dE^ (2.41) 

^■^e E ® e e 

P 


where E is the threshold energy. 

Jr 

This equation is integrated numerically using either 
trapezoidal method or Simpson's rule of integration. Results 
obtained with Simpson's rule are more accurate. 


2.3.2, Population Densities 

Applying the principle of detailed balancing in the 

n __ n 

processes used in the model and putting " - t ' ~ — = sind ~ = Y; 

+ ++ n (1) n 

where n (1) and n ^ respectively^ are the population densities 

of the groxind state of Hell and that of doubly charged helium icai , 

the following equilibrium relations are obtained. 


K(p/q) nj,(p) 

s= K(q/P) nj,(q) 

(2.42) 

n 

^ K(i/p)j^^^ 

“ K(p/i)^_j^ 

(2.43) 

K(i/j) ng(i) 

= K(j/i) ng(j) 

(2.44) 

n ■ ' ' ■ ■ 

rr K(c/i) 

== K(i/c) ng(i) 

(2.45) 


Here ng(p) and ng(q) are/ respectively/ the Saha equilibrium 
population densities of the pth and qth levels of Hel; nj^(i) and 

r ■ ■ ■ ■ ■ 

njjCj) are the Saha equilibriim population densities of the ith 
and jth levels of Hell. 
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For the processes described by Eqs. (2.1-2.12)/ the time 
derivatives of the Hel and Hell levels, after incorporating the 
Saha equilibrixjm population densities, may be written as follows: 




15 32 

p(p) { n [ I K(p,i) + I K(p,q) ] 


32 


i=l 


q/p 


+ n(l) I KN(p,q) + I A(p,q) + n(l) KN(p,i) . 

q/p p>q 

32 n^Cq) 32. 

+ Z p(q) KN(q,p) n(l) - ;-v + I p(q) K(p,q) n 




and 


ng(i) 


where 


ng(q) 


n 


q>P ng(p) "^e ^^P'^^i=l 5jng(p) 

(2.46) 


15 


= - P (i) {njK(i,c) + I K(i,j)] + I A(i,j) } 


i^j 


15 

I 


+ I p‘'‘(j) K(i,j)n + I p'''(j) A(j,i) 


±>i 

ng(j) 


j>± 


ngd) 


P < S(iV 

4- I p(p)n -r- — K(p,i.) + n K(i,G) + — . .-XJ : 

p=l ® n^(i) ® y r^d) 


(2.47) 


P(p) 




E 


P) 


and P^(i) *= ^ 

n^(i) 


(2.48) 


. mjfm 

n(p) and n (i) represent the population densities of pth level of 
Hel and ith level of Hell respectively. In the form of matrix. 



the rate equations (2.460 and (2.47) can be represented (Pigxire 
2.2) as 

= A p(p) + B + C (2.49) 

where is a coiximn matrix containing 47 elements. The 

values of p from 1 to 32 refer to Hel levels and from 33 to 47 
refer to Hell levels, p = 33 represents the ground level (i = 1) ; 

of Hell while p = 47 indicates i = 15 level of Hell. A is a \ 

. , f! 

47 X 47 matrix containing the coefficients of p's, B and C are | 

I 

column matrices^ each containing 47 elements which are indepen- [ 

dent of p's. ' 

i 

Under the quasi-steady state conditions^ it may be ! 

! 

i! 

assumed that the time derivatives n(p) for the excited states are j 
negligible compared to those for the ground states of Hel and ; 

Hell. Hence an inhomogeneous stationary state solution can be ! 

obtained by equating the right hand sides of Eq. (2.46u) for p 1 
and Eq. (2,47) for i 1 to zero. Thus equation (2.49) reduces to i 

ti' 

' : 1 
I 

A P(p) = - B - C (2.50) 1 

I 

' 

where A now is a 45 x 47 matrix^ the elements of which correspond 
to p a= 2 top “ 32 and p = 34to p = 47. It may be mentioned here 
that A is obtained by dr oping the first and thirty third rows of 
the old matrix. Similarly B and C, now contain only 45 elements 
each. 

Cto transferring the first and thirty third colxirans of A 
to right hand side and multiplying both sides by -1# one gets 
(Figure 2.3) 


- A P(p) 


B + C + D P(l) + E p(33) 


(2.51) 



o 
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, ■ ' f' ', ■ 

Figure 2.2. Mat^-ix representation o£ collisional-radiative 
' : rate processes, 















Figure 2, 3* Diagrammatic representation of the proceciure of 
calculating po^julation densities of He. I and Hell 
levels . ' ■ ' ■ ■ ■ 





The elements of A (a 45 x 45 matrix) correspond to levels p = 2 
to p = 32 and p = 34 (i = 2) to p = 47 (± = 15). The colxann 
vector p (p) contains p's corresponding to p = 2 to p = 32 and 
p s= 34 to p == 47. p(l) is given by n(l)/nj,(l) and P(33) by 
n (D/n^jd) where n(l) and n (1) represent the population densities 
of Hel and Ifell ground states. nj,(l) and nj,{l) are corresponding 
Saha equilibrium population densities . D is a column matrix whose ; 
elements aare those which are transferred fran the first column of 

A and E contains the elements of the thirty third coluinn of A. j 

1 

Eq, (2.51) can be rewritten as | 

i 

■ . ' : , 



p(p) = 

a“^ [ b + C + D P(l) + E P(33) ] 

(2.52) 

or 

P(p) * 

rQ(p) + r^(p) P(l) + r 2 (p) p(33) 

(2.53) 

where 

^o = 

a'“^(B + C)^ and r 2 = A~^E. 


where 

-1 

A 

Eq, 

represents the inverse of matrix (-A) . 

(2.53) represents the solution of Eqs, (2.46) 

and 


(2.47) in terms of local ground state population densities of Hel 
and Hell, Eq. (2.53) can be written separately for Hel and Hell 
levels in the following form 


P(p) = P(l) (2.54) 

p'*'(i)>= r^(i) + r^(i) P (1) + r^d) P'^(l) (2.55) 

Diagramatically this is shown in Figure 2,3, Eq. (2.54) gives 
the reduced population coefficients for Hel levels. It may be 
mentioned that this equation contains no term involving r 2 becaijse 
first 31 elements of E are zeros. r^(p), r^Cp) and r^ are termed 


as the 


population coefficients. represents the 

contribution from the first continuim towards p(p) and rj^(p) p(l) 
is the contribution from the groxind state of Hel. Eq. (2.55) 
gives the reduced population coefficients of Hell excited states. 
P”*^(i) stands for the values of P for He^ levels, r^z r^ and r^ 
correspond to r^^ and r 2 values of Eq. (2.53) for He"^ levels 

* 4 " 

and P (1) corresponds to P(33). is the contribution from 

the second continuum towards P^(i); p(l) is the contribution 

from the ground state of Hel i.e. the contribution from the 
direct -ionization-excitation process and r^Ci) P^(l) is the 
contribution from the ground state of Hell towards p"*^(i). It is 
to be noted that 


ng(p) 

njd) 


iE nlOi 


) 


(2.56) J 


Here g_ and g. are, respectively/ the weight factors of the pth 
P 

state of Hel and of the ith state of Hell, Prom Eq. (2.56) it is 

clear that the value of the ratio must be known in order to 

^ + 

obtain the final solution. The final solution for P(p)/ P (i) 
and population densities of excited levels of Hel and Hell are 
thus obtained in terms of n / T / gas temperat\xce , n(l)/ n'^d) 

n ^ ^ 

£-1^ (foa;; conputer programme see Appendix iv) . 


2.3,3, Role of Direct Ionization-Excitation 

The rate coefficients for direct ionization-excitation 
process (Eq, 2.2) are obtained by integrating the cross section 
expression (Eq, 2,23) over a Maxwellian velocity distribution of 
electrcns in the following way; 


K(p^l) 


-4 2 „ 1 2 m 

/* 46 X 10 -a^C^) 24 in(l 0 x) (|) (jJ-) 

V P-^ X © 


-1/2 mv /kT 


e 2 , 

. V .vdv 


(2.57) 


Here the terms involved have the same meaning as described in Eqs. 
(2.23) and (2,33). Following the procedure described in Eqs. 
(2.33-2.37)^ one obtains the following expression for K(p,i) 

00 


K(p,i) >= 6.25 X 10 


-3 


,3/2 


~a 


- / 

Ct J 


2 , 

a 


+ 1 


In {10.0(1 + |)} 


e"^ dZ 


(2.58) 


where 


kT 

e 


and 


Z = a (x-l ) 


Eq. (2,58) is solved numerically by using the Gaussian Quadrat \ire 
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method involving Laguerre integration. 


15 


m Eq. (2.46)^ the term - P(p) n 1 K(p/i) represents 

. , . ® i=l 

direct ionization-excitation and depopulates the excited states of 

Hel. In Eq. (2.47) the direct ionization -excitation process 

populates the excited states of Hell through the teirm 
32 

I p(p) n — 37 - — K<p,i). On dropping these terms from Eqs. (2.46) 
i=l ® bE(^) 

and (2.47) and solving the equations simultaneously^ the role of 
this process can be explored. 

2.3,4. Mechanism of Population of Excited States 

Quantitative calculaticais based on the collisional-^^^ 
radiative models use cross sections of various processes involved. 
Due to a pa\K:ity of experimentally determined cross-sections / 
collisional radiative models often lean heavily on cross section 
data which are either theoretical or empirical. Whereas# some 



times there are no alteimatives^ it is clear that use of these 
must be kept to a minimum; because, not infrequently, theoretical 
population c3ensities wherever they can be compared with experi- 
mental results are off by as much as 100% and little is 
achieved by including processes which are not of outstanding 
importance and whose experimental cross sections are unknown. 

As mentioned earlier, the various collisional radiative 
models, employ a very large nvimber of processes. Therefore, it j 
is worthwhile to neglect the processes whose role is insignificant 
in populating the excited states . On this basis a simplified 
model can be evolved. It will be shown in the next chapter that 
(adopting the procedure outlined below) even more than 70% of 

I 

the atCHtiic processes used in some of the recent models can be i 

neglected and yet obtaining theoretical popxilation densities 
substantially close to the values predicted from the models whicdi 
include all the processes. The following method also permits an 

I 

insight into the mechanism of the populaticxi of a particular level, i 
For a quasi-steady state solution with pop\ilation densities 
of the excited states derived using the population densities of 
ground levels of ffel and Hell as parameters, the first and the 
thirty third rows in matrices of Figure 2.2 are removed, and then 
the left hand side vector is made equal to zero. This results in 
45 l i near simultaneous equations. The first and the thirty third 
columns, which respectively involve transitions f ream the ground 
states of Hei and Hell, are then transposed (Bq. 2.51) . 

To identify the important elementary processes and thus 
to evolve the popvilation mechanism, we proceed with the matrix 
Eq. (2,52). On the left hand side we have matrix A of dlmensioh 


45 X 45 and on the right hand side we have four colximn matrices 
B, C/ and Ef each of dimension 45 x ly which for convenience 
may be represented by a single matrix F of dimension 45 x 4 
(Pigxare 2.3). The first column of this matrix P contains the 
three-body recombination coefficientSy the second column contains 
the radiative recombination coefficients. The third and the fourth 
columns are those transferred from the left side; these contain 
transitions from the ground levels of Hel and Hell. A three 
dimensional matrix G of dimension 45 x 45 x 4 is then formed in ‘ 


the following 

manner: 


G(iyjyl) 

= A“^(iyj).F(jyl),ng(i+l) 

(2.59) 

G(iyjy2) 

= A“^(iyj).P(jy2).ng(i+l) 

(2.60) 

G(iyjy3) 

= A'^Ciy j) .F( jy3) .n(l) .ng(i+l)/ng(l) 

(2,61) 

G(iyjy4) 

= A“^(iyj).F(j,^).ng.ng(i-fl)/ng(33) 

(2.62) 


where both i and j vary from 1 to 45 and represent both Hel and 
Hell levels, represents the inverse of matrix (-A)y and nil) 

and n(33) respectively represent the population densities of ground 
states of Hel and Hell. For population densities of Hell levels 
(l.e. 34 i 1 47) the factor (i+l) in Eqs. (2.59-2.62) is to be 
replaced by (i+2) , Diagramaticallyy this is represented in the 
last equation of Figure 2.3. 

The sum of the elements of the ith row of all the four 
planes of G Gonstitutes the population density of the 1th level. 

The elements on the first plane of G represent contributions from 
three body recombination and those on the second plane represent 



contributions from radiative recombination. The elements of the 
third and fourth planes represent contributions from the groiand 
states of Hel and Hell respectively. In. order to identify the 
dominant terms, say for level 1 population density, we start 
scanning from the ith row of the first plane of G and scan all 
the elements of the corresponding rows on the four planes. We 
set all those elements equal to zero whose contributions to popu- 
lation density is less than 5%, This sorts out those terms of F 
and those of the inverted matrix which are important and also 
provides relative contributions of the four segments to population 
density of a particular level. The insignificant elements of F 
(corresponding to small A^^.F values) are then set equal to zero. 
This is done for all the levels and this process identifies the 
important elements of F. 

For sorting out those rate coefficients which make 
dominant contribution to a particular element of the original 
matrix A, the following procedure is adopted. The diagcxial element; 
of the original matrix contain those colli si onal rate coefficients 
which contribute to depopulating a level. We scan all the rate 
coefficients contained in a diagonal element and retain only, 
those rate coefficients which contribute significantly. For ncoi- 
diagonal elements of the original matrix, we start elementwise 
fron the first row of a, a chosen non-diagonal element is then 
set equal to zero., matrix A is again inverted and a corresponding 
matrix Gl of dimension 45 x 45 x 4 is formed after Eqs. (2.59- 
2.62). A comparison of the non zero elements of G with the 
corresponding elements of G1 brings out the quantitative role of 
the element which is set equal to zero in the original matrix. 



If any element of G1 differs from the corresponding element of G 
by more than 15% (an arbitrary limit )y it is assumed that the 
element of A which was made equal to zero is impoartant and is 
reinstated for all further operations ^ otherwise it is left as 
zero. This process is repeated over all the elements of A. The 
individual rate coefficients and transition probabilities cont- 
ained in the remaining non -diagonal elements are further scanned 
and only those which contribute significantly are retained. A 
is then left with only important rate coefficients and transition 
probabilities in diagonal and non-diagonal elements. In order to 
identify t^ processes which are important in determining the 
population densities of individual levels ^ we proceed as follows. 
In matrix A (now containing only the important terms ) ^ one parti- 
cular rate coefficient is set equal to zero and the resulting 
niatrlx is inverted Those elements of the inverted 

matrix which are Important for populatlai of a particular 

level are compared with the corresponding elements of the inverted 
matrix A (with no rate coefficient equal to zero). If the 
inverted matrix elements of a”^ show significant departure (say 
more than 10%)^ we assume that the particular rate coefficient, 
which was set to zero, to be important for all those levels whose 
corresponding matrix elements show considerable depart xires. This 
is repeated with all rate coefficients and transition probabi- 
lities and this yields the mechanism of population. 

2.3.5, Laser Induced Selective Excitation 

If the plasma is stibjected to a tunable laser radiation 
resonant between a lower state p and an upper state q of Hel, the 
radiative processes induced can be written as 



He(p) + hv He(q) 


(2.63) 


where W(p,q) Is the rate of induced absorption and W(q,p) j^epre- 
sents the rate of induced emission. The induced transition rate i 
may be obtained by multiplying the radiation density per unit 
frequency interval p(v) with the Einstein absorption coefficients 
corresponding to the transition and integrating the prodvtct over 
the line shape functions of the transition as well as that of 

! 

the laser pulse . We must mention here that due to non eqtrillbEi.um 

I 

situation# the commonly used A# B (Einstein coefficients) rela- 
tionship is not strictly applicable# though it has been losedk | 

51 » 

in recent literature. In a later section we shall discuss tl^ 

possible nature of the departures. 

(a) Laser power and radiation density; Radiation 

density P(v) is defined in terms of the energy per unit volume i 

2 

per unit frequency interval. Let I be the intensity (in watt/cin ) 
of the laser pulse. This means that I amount of energy is passing 

through 1 sq. cm. area in 1 sec. Therefore in 1 sec# I energy will t 

-3 i 

contained in volume c cm (1 sq. cm x c)#,. where c is the velocity , 

of light. Therefore density is given by | 

P = J (2.64) i 

I 

/Av' 

Ptirther if the full width at half maximum of the laser pulse ^is 
larger than the width of the transition Av # the energy density per 
unit voltme per imit frequency interval may be given by 

1/Av* 
c 


P (v) 


(2.65) 



(b) Rate of Induced transitions: The Einstein coeffi- 


cients for induced absorption B(p,q) and for induced emission 
B(q,p) are given by 

. c^ 

B(q.p) = ^ . A(q,p) (2.66) 

8 TT h ^>P 

% 

and B(p,q) = ^ . B(q,p) (2.67) 

where v is the frequency of the transition q -► p and g and g 

P 5 

are the weight factors of levels p and q respectively. 

■ ■ . i 

The induced transition rate W (q^p) due to interaction | 
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with the field is given by : 

W (q#p) = B(q,p) p (v) 

, j 

^3 

= . A(q.p) . p(v) (2.68) | 

8 ff h v‘ f 

t 

TO derive the induced transition rate W(q/p) due to a single | 

j 

frequency field, it is postulated that the transition rates due j 

■ i? 

to different frequencies in a multifrequency radiation field are 1 
additive so that Eq. (2.68) is a special case of a more general | 

form given by 

W (q,p) = / -K p(v') g(v' - v_) dv' (2.69) 

-00 Sir h ° 

where g(v' - v^) is the normalized line shape function (centred 
at v^) for the transition, is the frequency corresponding to 
the peak of the radiation field. 

In a monochromatic radiation field of frequency v and 
radiation density p the energy density per unit frequency P(v' ) is 

r.v .. 


p ( V ) 


5 ( v' - V ) 


(2.70) 


where 6 ( v' - v ) is a Dirac delta fianction. Eq. (2,70) when used 
in Eq. (2.69) gives 

3 

W(q/p) = — — T A(q,p) g( v- v_) (2.71) 

Stt h v-^ ° 


The rate of induced absorption similarly is given by 


g 3 

W(p^q) (^) 2 — « , A(q,p) .p 

^p 8Tr h V 


g(v - v^) 


(2.72) 


It is assumed that the laser pulse is of Gaussian shape. 
The Gaussian line shape results most often in absorption or 
emission from gaseoxis pairticles. The resulting radiation is 
Doppler shifted owing to the velocity spread of the emitting 
particles. For a gas in thermal equilibrium^ the velocity distri- 
bution function is Maxwellian and the resulting frequency distri- 
bution fvinction g(v - v^) is a Gaussian. The normalized Gaussian 


can be written as 


g( V - v^) 


2(ln 2) 


V - V 2 

1/2 -4 (in 2)( -■) 


1/2 , 
t A\) 


(2.73) 


Where Av is the full width at half maximum (FWHM) and is given by 

/21c T 

Av = Av ^ ■ = 2v / ^ In 2 

Doppler o ^^^2 


where T^ is the gas temperature and m^^^ the mass of helium atom. 

(c) Enhancement in population densities: When a laser 
pulse is timed to a particular transition (say p' q')/ the 


population densities of most of the other levels do not change 



much. Only the population densities of the upper level q' and few 
other levels are enhanced. The tinne derivatives of the excited 
states other than p' and q' can still be described by Eqs. (2.46) 
and (2.47). The time derivative of the lower level p' is given 
by the following equation: 


n(p' ) 

h;tf) 


15 32 

P(p') {n [ I K(pM) + I K(p',q)] 


i=l 


q?^P* 


32 

+ n(l) I KN(p'^q) + I A(p',q) + n(l) KN(p‘ ,i) 
qf^p' P'>q 

+ W(p',q') L(t)} +p(q') W(q',p')/L(t) + I p(q) ■ 

P'^ 

n-(q) 32 

KN(q,p') n(l) r-V ' Y + I o(q> K(q,p')n 

+ I p(q) A(q,p') n 
q >P E 


n. 




(2.74) 


where L(t) is the shape function of the laser pulse. Similarly/ 
the time derivative of the upper level q' can be obtained by 
replacing p' by q' and q' by p' in Eq. (2,74). 

Eq. (2.46) (for p ^ p' and q')/ Eq. (2.47) and Eq. (2,74) 
(for p* and q') form a set of 47 linear first order differential 
equations. On multiplying these equations with the correspojiding 

n_(p) values/ one gets the following equation in the matrix form 

E 

^ n(p) e= A . n(p) +B (2.75) 

where n(p) is a vector of population densities containing 47 
elements. A is a 47 x 47 matrix/ which is a function of time and 


B is a coluinn matrix containing 47 elements and^is independent 
of time. 


Assuming that the electron density does not change with 
time, the population densities as a function of time can be obtained 
by solving the system of equations described by Eq, (2.75). In 
this work, the system of 47 coupled linear first order differ- 
ential equations is solved numerically by using Hamming's modified 
predict or -correct or method. It is a stable fourth order inte- 
gration procedure and obtains a fairly accxjrate solution of the 
system of coupled differential equations with given initial values. 
The initial population densities of all excited states are taken 
as the quasi-steady state population densities. Hamming's 
predictor -corrector method is not self starting i.e. the functional 
values at a single previous point are not enough to obtain the 

functional valves ahead. Therefore, to get the starting values, 

87 

a special Rvinge Kutta procedure followed by one iteration step 
is added to the predict or-correct or method. The solution of Eq. 
(2.75) gives the population densities of all excited states 
(including those optically pvomped) as a fvinction of time. The 
enhancement of the population density n ,of level p' is defined 

Jbr 

as follows: 


E ,(t) 
p' 


(t) - n^, (t = 0) 

Up, (t =0) 


(2.76) 


where n , (t = 0) is the population of level p' before the laser 
pulse is injected. This initial state is assvimed to be a quasi- 
steady state solution of Eqs.( 2. 46 d) and (2.47). 
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In the time dependent solution of Eq. (2.75), only the 
laser flux is assumed to be a time dependent quantity. The shape 
of the laser pulse L(t) is given by 

L(t) =1 for 0 < t < T (2.77) 

L(t) = T<t (2.78) 

T 

where ^ is the FWHM (nsec) of the laser pulse. The subroutine 
used to solve the linear system of first-order coupled differ- 
ential equations is given in Appendix V. 



CHAPTER 3 


RESULTS OP CALCULATIONS 

3.1. RATE COEFFICIENTS 

Typical values of rate coefficients for various types 

of transitions of Hel and Hell calculated cai the basis of procedure 

described in Chapter 2 for a wide range of kT^ are presented in 

Tables 3. 1-3. 4. Electron impact excitation rate coefficients for 

optically allowed and optically forbidden transitions of Hel are 

presented in Table 3-1 in the kT^ range 1-20 eV. Table 3.2 

0 

presents electron impact ionization^ direct ionizaticai-excitation 
and radiative recombination rate coefficients for Hel. The rate 
coefficients for heavy particle excitaticxi^ deexcitation/ ioni- 
zation and recofribination are listed in Table 3.3. Table 3.4 
presents rate coefficients for electron impact excitation / ioniz- 
ation and radiative recombination for Hell in the kT range of 

0 

1-20 eV, 

3.2. OPTICAL ESCAPE FACTORS 

Values of optical escape factors for the resonance 
transitions of Hel and Hell calciilated using Eq. (2.32) are 
presented in Table 3.5. The values corresp<and to = 300 
n(l) e 2.0 X lO^'^ cm“^^ 1 = 0.84 cm and n^ = 10^°, 10^^ and 10^^ 
cm . The corresponding values of a and t are also given in 
Table 3.5,: . 

3.3. POPULATION DENSITIES 

The computed values of reduced population coefficients 
(results partially presented in Ref. 56) r^ and r^ for Hel are 



Table 3.1. Electrcai impact excitation rate coefficients for Hel (cm"^ sec 
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Table 3,2. Electrai impact ionization/ direct ionization-excitation and radiative 
coefficients for Hel (cm^ sec“^) 
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Table 3.3, Heavy particle collisional excitation # deexcitation 
ionization and recombination rate coefficients 
(cm^ sec“^) 

(kT = 10 eV, = 300 °K and 1000 °K) 

® Q 


Type of Transition 

1 

1 

f 

1 Transition 

1 

1 

i T = 300 °K 
. Q 

1 .. - 

1 

; T = 1000“K 
I 9 

1 

.1 - - - - 


3^D-3^P 

2.23-10* 

1 

1.66-09 

Neutral-neiiitral 
excitation KN(p,q) 

3 1 

3-^D-3 P 

1.09-10 

8.12-10 


1 1 

3 S-3 P 

4.24-14 

1.09-11 


3^P-3^D 

6.15-10 

3.22-09 

Neutral-neutral 
deexcitation KN(q,p) 

1 3 

3 P-3"^D 

9.10-10 

4.74-09 


3^P-3^S 

8.94-12 

2.53-11 


n=6 i=l 

2.34-16 

2.02-11 

Neutral-neutral 
ionization KN(p,i)j^_^ 

n=8 i=l 

9.38-13 

1.15-09 


n=10 i=l 

4.47-11 

8.78-09 


i=l n=6 

7.44-36i 

1.79-31 

Recombination 

KN<i'P>i=l 

i=l n=8 

5.59-32 

2.58-29 


i=l n=10 

3.11-30 

2,62-28 


-10 


* Read 2.23-10 as 2.23 x 10 


Table 3.4. Rate coefficients for Hell (cm' 
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Read 1.2-25 as 


T^le 3.5. optical escape factors for rescMiance transiticais of Hel and Hell 

(n(l) » 2.0+14 = T. = 300‘*K, 1 = 0.84 cm) 
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presented in Figures 3. 1-3, 8 and the values of r^/ r^ and r]t 

o 1 2 

presented in Figures 3.9—3,14. The data presented in these 

figures refer to kT = 1.38 ev (1.6 x 10^ °K) or kT =10 ev, 

Tg = 1000®K^ (n /n^^) = 1000, n(l) = 1.0 x 10^"^ cm”^, and n^ n^. 

The electrcn temperature kT = 1.38 eV and 10 eV were .chosen to 

0 

allow a comparison of calculated results with the results of Refs. 

51 and 52. In Ref. 51, the results are presented at kT^ = 10 eV 

while in Ref. 52 the results are given at kT =1,38 eV. 

e 

Figure 3.1 shows a plot of log r^ vs electron density at 

kT^ - 1*38 eV for optically thin conditions. The dotted lines 

52 

represent the results of Fujimoto^, . Figure 3.2 presents a plot 

of log against electron density at kT^ = 10 eV under optically 

thin condition. The r^ values for the 3^s, 3^p and 4^F levels 

51 

are compared with the values of Hess and Bxirrell. Using 

partially optically thick conditions, the values of r^ increase 

considerably. This is evident from Figures 3.3 and 3.4 where 

log r^ values have been plotted against n^ at kr^ = 1.38 eV and 

kT « 10 eV, respectively. Figure 3.5 presents log r, values 

plotted against n at kT = 1,38 eV for optically thin conditions. 

At kTg = 10 eV, log r^^ values are plotted against n^ for the 

o;ptically thin case in Figure 3.6. The dotted lines correspond 

1 51 

to the values of 3 D frcxn Hess and Bxjrrell. Figures 3.7 and 
3.8 are plots of log r^^ vs n^ for the partially optically thick 
case at kT^ =» 1.38 and 10 eV, respectively. The values of r^^ 
also increase in the transition f ran the optically thin case to 

■ "V' 

the partially optically thick case. 

Figures 3,9-3.11 present reduced population coefficients 
of Hell levels at kT^ « 1,38 eV under optically thin conditions . 




Figure 3 . 1 , 


Plot of log for Hel as a function of electron 
density at kT^ » 1-38 ©V fc the optically thin 
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Pigur« 3.3. Plot of log for HeT vs electron density, at 

kT “ 1.38 eV for the partially optically thick 
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Figure 3.4. Plot of log for Hel vs electron density at 
kT » 10 eV for the partially optically thick 









Figure 3.6, Plot o£ log for Hel vs the electron density 

* 10 eV for the optically thin case? 

IT , , R©£ • 51#, 
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Figure 3.7* Plot of log for Hel vs electron density at 

jcT^ * 1*38 ev for the partially optically thick 


case . 
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Pigur« 3,10* Plot of log Jtoll vs electron density at 

JcT* “ 1-38 eV for the optically thin case. 
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rigur* 3.13. Plot of log for Holl vs ©loctroo danslty at kT ■> 1.38 ®V for 

partially optically thick caaa©. 
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Pigure 3.14. Plot of log rj for Hell v» «l«ctraci density At *• 1*38 sV for 
; 'pArtlally optically thick case, ^ ■ 
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Figure 3.9 presents values of r^ for the Hell levels ± = 4-15, 
the r^ values are shown in Figure 3.10 and r^ is plotted in 
Figure 3.11, For partially optically thick conditions, the 
values of r^, r^ and rt for Hell levels at kT =1.38 eV are 
plotted in Figxores 3.12-3.14. 


3.4. ROLE OF DIRECT-IONIZATION EXCITATION 


In Ref. 11, possible importance of the direct ionization 
excitation process in populating the Hell levels was indicated- 
Here we make a detailed calculation on the role of this process. 
In Eq. (2.55), at low gas pressures (< 1 torr), the first term 
rQ(i) is essentially the contribution from the second continuxam 

4* -h 

towards p (i) , the second term, r^ix) p(l), is the contribution 

from the ground state of Hel, i.e. the contribution frcan the 

direct ionization-excitation process. The third term, ^ 

represents the contribution from the ground state of Hell . 

57 

Rewriting Eq, (2.55), one gets 


n‘^(i) 


r^(i) ng(i) + r^(i) 


n(l) 


. ng(i) + r2(i) n'^'d) 


n^(l) 


(3.1) 

where n"'^(l) is the population density of the ground level of Hell 
which is approximately taken equal to the electron density 

In the results presented in Figures 3.15-3.19 the plasma 
parameters used are n„ = lO^-lO^"^ cm"*^, kT_ = 0.1-18 eV, n(l) = 

^ ^ 4- ' 

1.0 X 10^^-2.25 X 10^^ cm"^ and ^ = 1-10^. Figure 3.15 shows 

n 

the individual contributions of the first, second and third terms 
of Eq. (3.1) for i = 4 at n^ = 10®, 10^°, 10^^ and 10^^ cm“®. 

c!' 

■ 8 "^3 

At n^ =10 cm and kT^ < 3.0 eV the contribution of the third 
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Figure 3.15, 


Plot of contributions front the first/ second and 
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Wiqnvm 3.17, Piot of contributions from th« socond and third terms of 

m* (3.1) for^i » U as a function of electron temperature 
n, • 10 , 10 ^ , lO^ Vand a « second tom, t « 

: third ter». 
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Eq. (3.1) for A •» 15 as a function of electron temperature 
at n_ * 10®, 10^*^, 10^^ and 10^^ » ■ eecond term, 

t • third terra. 
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term is larger than the contribution of the second term. For 

kT > 3 eV , the reverse is true. At kT = 5 eV, the second term 
© © 

is about two orders of magnitude and at kT^ = 15 eV it is about 

three to four orders of magnitude larger than the third term. 

10 ~3 

At n^ =10 cm , the direct ionization-excitation is important 

only above kT^ = 6.3 eV. At this higher electron density# the 

direct ionization -excitation becomes relatively less important 

8 "“3 

compared to the case at n^ = 10 cm . This is evident from the 
fact that at kT^ =15 eV# the second term is only about one to 

11 

two orders of magnitude larger than the third term. At n = 10 

-3 

cm # ionization-excitation process is important above kT^ = 10 

eV. At kT^ =15 eV# the second term is about one order of 

12 -3 

magnitxjde larger than the third term. At n = 10 cm # the 

third term dominates over the second term up to 18 eV# the upper 

12 -3 

limit of kT used in the computation. Thus at n = 10 cm # 
the i = 4 level of Hell is populated mainly from the ground 
state of Hell and not from the ground level of Hel. At low 
electron temperatures ( <2 eV) # the first term of Eq. (3.1) 
dominates over the second and third terms# which means that the 
contribution from the second continuum is the main factor in 
populating Hell excited levels at low electron temperatures. 

Pig\ares 3.16-3.18 present similar plots for Hell levels 
i = 7# 11# and 15. log of second and third terms of Eq. (3,1) 
for Hell level i = 7 are plotted against kT^ in Figure 3.16# for 
level i =11 in Figure 3.17 and for level i =15 in Figure 3.18. 

From Figures 3.15-3,18# it is evident that for a parti- 
cular Hell excited state a range of electron density exists where 
the role of ionization-excitation is considerable and above a 
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particular electron temperature this process plays a dominant 

role- These findings have been summarized in Figxxre 3.19. This 

figure shows ^ fori=2toi=15 and in the electron density 
8 1 3 ••3 

range 10 -10 cm ^ the electron temperatiire above which the 

direct ionization-excitation dominates. At low electron densities 

(10 -10 cm ) and at low electron temperatures, the direct 

ionization-excitation becomes important. At higher electron 
10 -3 

density (10 cm ), the ionization-excitation is important at 

.ij*3 

relatively higher electron temperature. Up to n^ = 10 cm'^'^, 

the direct ionization -excitation is the main process for populating 

all Hell excited states from i=2toi=15. The electron 

temperatxare above which the direct ionization-excitation is 

important is lower for higher excitation. At n^ = 10 cm’" and 

for the level i = 2, the ionization-excitation process is not the 

major population mechanism up to kT^ = 18 eV, the highest value 

for which the computations were carried out- At still higher 

electron densities and kT < 18 eV, the direct ionization- 

e 

excitation is important only for higher Hell levels. It is to be 

noted that at all electron densities, the most of dcminance o| the 

ionization -excitation process occurs for the higher levels at 

relatively low electron temperatures. 

Table 3.6 presents the population densities of i = 4 level 

at n_ = 10®, 10^®, 10^^, 10^^ and 10^® cm"®, calculated both with 

and without direct ionization-excitation process at various 

electron temperatures. The data given in this table correspond 

to n(l) = 2.25 X 10^^ cm"® and - - = 1000, From this table it 

n 

is clear that the population densities increase considerably (one 
to three orders of magnitude) on adding the direct ionization- 
excitation processes. 
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N5 {^,11 

/ I 




wlgnm 3.19^ 


. Plot of the electron teinperatmre's ■above which .tlte ' 
direct ionization-excitation process is dominant, for 
various Hell levels, as a function of electron density, 





Table 3,6. Population densities (cm ) for i = 4 level of Hell 
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3.5. MECPIANISM OF POPULATION 

Using the procedxare outlined in Chapter 2^ the collisional- 

radiative model is optimized at kT = 12 eVy n = 10^^ cm”^/ n(l) = 

^ © 

2.0 X 10^"^ cm”^ and T = 1000 °K. On optimization^ it is foiind 

y . ■ 

56 

that only 30% of the atomic processes of the elaborate model are 

sufficient (results partially presented in Ref. 58) to yield 

population densities within 15% of the results of Ref. 56. 

For most of the levels, the major contribution to 

population densities is found to originate from the Hel ground 

state through electron impact excitation process. 

Individual rate coefficients and transition probabilities 

of the processes affect the elements of inverted matrix A 

(Eqs. 2.59-2.62) considerably. Table 3.7 shows the quantitative 

effect of a few rate coefficients on the elements of the inverted 

matrix. On using the important rate coefficients and transition 

probabilities (obtained by the procedure outlined in Chapter 2) 

the calculated population densities obtained are compared with 

the population densities obtained from the elaborate model in 

Table 3.8. To test the validity of this simplified model (having 

only 30% of the total processes) in a range of plasma parameters/ 

the population densities were calculated using the simplified 

model in the electron temperature range 4-20 eV, electron density 
,6 ' .12 "^3 

range 10 -10 cm and pressure range 0, 01—1 torr. The percent 

change in population densities compared to those of the elaborate 

56 ^ ^ ^ ^ 

model is presented in Table 3,9. From this table it is evident 
that the simplified model containing only 30% of the total terms 

can predict the population densities within 20% of those obtained 

■' 56 ' ' ' 

frcxn elaborate model in a wide range of electron density and 


e lectron temperature . 



Table 3.7. Effect of dropping some electron impact rate coefficients on elements of inverted 
matrix 
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2(21) means that the value of the element 2/2 of the inverted matrix is increased by 21% 
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Table 3.8. Population densities of Hel and Hell levels obtained 
from elaborate^^ and simplified models 

(kT^ = 12 ev, n = 1.0+12 cm"^, n(l) = 2.0+14 cm“^) 

0 0 


Level 

I Population Densi- 
I ties (Gm~3) from 

I Elaborate Model 

-L_ 

1 

i 

; Population Densi- 
I ties (cm'*^) from 

I Simplified Model 

1 

1 . - 

Percent Change 
in Population 
Densities 

.... - 

2^S 


1.854+11 . 

1.674+11* 

-9.7 

2 S 


4.591+10 

4.182+10 

-8.9 

2^P 


9.721+07 

8.392+07 

-13.6 

3 S 


4.183+08 

3.897+08 

-6.8 

3^P 


9.300+08 

7,926+08 

-14.8 

3^D 


2.104408 

1.913+08 

-9.1 

3^D 


1.658+08 

1.397+08 

-15.7 

3^P 


7.732+07 

6.811+07 

-11.9 

4^S 


2.345+08 

2.224+08 

-5.2 

4^S 


3.917+08 

3.573+08 

-8.8 

4^P 


2.759+08 

2.695+08 

-2.3 

4^D 


1,159+08 

1.027+08 

-11.4 

4‘^F, 

4^P 

6.092+07 

5.380+07 

-11.7 

5^D 


2.359+07 

2.351+07 

-0.3 

5^P 


3.396+07 

3.366+07 

-0.9 

If 

a 


5.751+07 

5.708+07 

-0.7 

11 

o 


2.503+07 

2.868+07 

14.6 

n=12 


3.121+07 

3,712+07 

18.9 

Hell 

i=2 

8.512+04 

8.417+04 

-1.1 


i=5 

4.302+03 

4.292+03 

-0.2 


i=10 

4.161+03 

3.918+03 

-5.8 


i=15 

7,071+03 

6.809+03 

-3.7 


* Read 1.674+11 as 1.674 x 10^^. 
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The population densities of 18 lowest levels of Pfel were 
also calculated by truncating the elaborate model to contain 
only 18 levels of Hel besides the ionic ground state. The 
mechanism of population, i.e., the processes which contribute 
significantly to the population of a particular level, evolved 
from this f\irther simplification is presented in Table 3,10. 
Absolute population densities and percent change in population 
densities compared to those of the elaborate model are also given 
in Table 3,10. The population densities of most of the levels 
are again within 20% of those predicted by the elaborate model. 

It may be added here that 18 levels of Hel would still provide a 
range of lines for spectroscopic measurements. This model, 
containing only 18 levels, consists of only 6% of the total ntomber 
of processes of the elaborate model. The size of the model can 
ftirther be reduced, without sacrificing much the accuracy in 
calculation of population densities, by taking into account only 
10 levels. The results are shown in the enclosure in Table 3.10. 
This latter model consists only 2.4% of the total terms of the 
elaborate model and yet predicts the populations of first ten 
levels quite closely. 

3.6. LASER INDUCED SELECTIVE EXCITATION 

As discussed in Chapter 2, the process of optical pumping 
by a laser pulse is included in the collisional-jradiative model 
and absolute population densities as a function of time are 
calculated. A typical plot of the enhancement and decay of 
absolute population densities against time for the pumping of 
2^P-4^D level by a laser pulse of power lO^-lo"^ W/cm^, FWHM = 5 
nsec, spectral width 0.1 A is shown in Figure 3,20. The data 
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12 —3 

presented in Figure 3.20 correspond to T =10 eV, n = 10 cm , 

0 0 

Tg = 300 °K and optically thin conditions. Continuous curves 

3 2 

correspond to calculations involving laser power = 10 W/cra and 

n(l) = 1,0 X 10^"^ cm Dotted line curves are obtained by using 

4 2 14 

laser power = 10 W/cm and n(l) = 1.0 x 10 cm while the 

dashed and dotted line curves correspond to data obtained by 

using laser power = 10^ W/cm^ and n(l) = 1.0 x 10^^ cm”^. On 
3 3 

pumping 2 P-4 D transition^ the population densities of levels 
3 3 1 

4 F, 5 P and 4 D also change considerably. An increase in the 

laser power shifts the enhancement curves towards left on the 

time axis. From the figure, it is evident that the excitation 

3 1 

transfer between 4 D and 4 D level is more pronounced at higher 


gas pressvire. 


CHAPTER 4 


DISCUSSION 

The procedure of calculations and the results of calcu- 
lations have already been presented in Chapters 2 and 3 respect- 
ively. In this chapter we critically examine the results obtained 
in the present work in the perspective of results of other workers^ 
and derive conclusions. The order of presentation is the same as 
the order of results presented in Chapter 3. First a discussion 
of cross sections and rate coefficients and their effect on calcu- 
lated population densities is presented (Section 4.1); and this is 
followed by a discussion of the optical escape factors in the 
light of calculations of other workers (Section 4.2). Then the 
calculated reduced population coefficients and populatiOTi densities 
are compared with the calculated and experimental results of other 
workers (Section 4.3). Section 4.4 examines varioias factors which 
affect the role of direct ionization excitation process. In 
Section 4.5, a discussion on the simplified model is presented. 

This is followed by (Section 4.6) an examination of results of 
the laser induced selective excitation experiments and finally 
(Section 4.7) a discussion on the line ratio method of determin- 
ation of electron temperature. 

4.1. CROSS SECTIONS AND RATE COEFFICIENTS 

In the present work, while choosing the cross sections, 
emphasis has been given on vise of experimental cross sections; 
the theoretical cross secticais are used only when the experimental 
cross sect icns are not available. Scanetijmes there is considerable 
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difference between cross sections obtained by various workers. 

Figures 4. 1-4.3 present a comparison of cross sections for electron 
. Ill 

xmpact excitation of levels 2 P, 3 P and 4 S from the groxand state. 

In Figure 4.1 the electron impact cross sections for the transition 

1 1 65 62 61 

1 S-2 P of Drawin , Scott and McDowell / Donaldson et al , 

Moustaffa Moussa et al^^ and Westerveld et al^*^ are compared. The 

62 

peak value of the theoretical cross section of Scott and McDowell 

is larger than the experimental cross section of Donaldson et al^^ 

71 

and Moustaffa Moussa et al ^ but is smaller than the experimental 

cross section of Westerveld et al^^. The empirical formula of 
65 

Drawin for the cross section of optically allowed transiticxis 
gives the largest peak value. The experimental cross section of 
Westerveld et al^^, being the latest one^ has been used in this 
work. Figure 4.2 compares the cross sections of various workers 
for the transition 1^S-3^P and Figxire 4.3 presents a comparison 
for the optically forbidden transition 1^S-4^S. In Figure 4,3^ 

62"' 

the peak of the theoretical cross section of Scott and McDowell” 

appears at about 30 eV impact energy, whereas in the experimental 

70 

curve of St. John et al it appears at aboxjt 45 eV, However, the 

magnitude of the cross section at peak is not much different. The 

68 ' ■ ' 

peak of the experimental cross section of van Raan et al (xised 
in the present work) is in between the above two cross sections. 
This cross secticai yields a better value of the population density 
of level 4^S, 

The choice of cross sect icais determines the value of rate 
coefficients and population densities of levels involved. Table 
4.1 compares rate coefficients and population densities for some 
electron impact processes obtained on using different cross 
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Figure 4.1. Total tarosa Motions £or tbs slsctron liqiMOt - sxxjitMian 

1- 1 2 . s* ■ 
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sections. All rate coefficients listed in Table 4,1 have been 

obtained by integrating cross sections over Maxwellian electron 

energy distribution. For a given transition, with a larger area 

under the cross sectional cxirve, the rate coefficient obtained is 

higher and consequently the population density of the upper level 

involved in the transition is also high. The population densities 

12 

given in Table 4.1 have been calculated at kT = 10 eV, n = 10 

0 0 

cm ^ and n(l) = 2.25 x 10^^ cm""^. This table examines various 
cross sections and their effects on rate coefficients and popula- 
tion densities. For optically allowed excitations from the ground 
state and for optically forbidden transitions without change in 
multiplicity, various cross sections are available in the liter- 
atiore and one can choose the most suitable cross section by 
comparing them critically. However, for electron impact transi- 
tions between excited states and also for optically forbidden 
transitions with change in multiplicity, the choice is very 
limited and often one has to depend on empirical or semiempirical 
cross sections. Similarly, for electron impact ionization, 
radiative recombination and heavy particle collisions , one has to 
depend on empirical, semiempirical, or theoretical cross sections. 

4.2. OPTICAL ESCAPE FACTORS 

8 ' 

The optical escape factors calculated from Holstein's 

expression are small as pointed out by otsxdca et al^*^. In the 

83 '■ 

present work the expression given by Drawin and Emard has been 
used for the calculation of optical escape factors. This expre- 
ssion considers Stark broadening and uses a Voigt profile to 
calculate the optical escape factors. The optical escape factors 
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calculated by this methcxi agree well with the values obtained by 
Otsuka et al^*^ both for Hel and Hell transitions iinder various 
experiinental conditions. In Ref. 10 optical escape factors have 
been calculated by nxjunerical integration of the expression for an 
infinite cylindrical geometry and the Doppler profile. The 
particle distribution is assumed to have "a bell-shaped form" 
proportional to that of electrons. For a cylindrical geometry, 
Otsuka et al^*^ give the following expression. 


^ 00 t / 

^k^'^o^ = ^ exp[- j- 


^ exp(- x^) ] 


/rTT2 

A 

For larger t this expression approximately behaves as 


(4.1) 



0.746. 
T ( Tf logx 


A 

where t is the optical depth. 

83 10 

A comparison of Drawin and Emard's amd Otsuka et al's 

A 

values of OEF for cylindrical geometry and for given values of 

is given in Table 4.^. Otsuka et al’s^^ values of optical escape 

factors are lower boxinds, and an use of additional broadening 

(e.g. Stark broadening) causes an increase in the escape factors. 

83 

Following Drawin and Emard , the Voigt profile was used in this 

work and optical escape factors were calculated for the various 

38 

conditions of TPD plasma machine experiment. These values are 
compared with those of Otsuka et al^*^ in Table 4,3. Although the 
present values are somewhat higher (due to inclusion of additimal 
broadening) than the corresponding values of Otsuka et al^^, yet 
the overall agreement may be considered to be reasonably good. 
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83 10 

Table 4.2, Comparison of Drawin and Emard's and Otsuka et al 
values of optical escape factors for a cylindrical 
geometry and Doppler profile 


1 

t 

4*^ * 

T ; 

1 

I 

O.E.F. from Drawin and 1 

Emard83 ! 

1 

- ,, . . ,....- 1 . , 

O.E.F. from 

Otsuka et al^® 

1.0-01 

9.320-01* 

9.007-01 

1.5-01 

9.000-01 

8.573-01 

1.0+00 , 

5.138-01 

4.211-01 

2.0+00 

2.894-01 

2.197-01 

1.0+02 

2.531-03 

1.951-03 

1.0+03 

2.083-04 

1.615-04 


* Read 9.320-01 as 9.320 x lO""^. 






Table 4.3. C^tical escape factors for Hel and Hell transitions for the TPD plasma experiment 
conditions, calculated by using Eq. (2.32) 
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Read 7.35-3 as 
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4.3. POPUIATiai DENSITIES 

The values of reduced population coefficients plotted 

against electron density at kT = 1.38 eV and kT = 10 eV under 

© © 

optically thin conditions are compared with the corresponding 

52 51 

values of Fujimoto and Hess and Bxarrell (shown by dotted 
lines in Figures 3,1 and 3.2 respectively). In general, the 
values are about 1-2 orders of magnitude higher than the corres- 
ponding values of Hess and Burrell^^. The r^ values for 4^F are 
10% larger than those of Ref. 51 at n^ = 10® cm~®. The differ- 
ences between the r results at lew electron densities of Refs. 

o 

51 and 52 and those of the present work may be traced to the 

choice of radiative recombination coefficients 6(p) used in the 

calculations. Whereas authors of Refs. 51 and 52 use empirical 

and semiempirical cross sections, the present work uses theore- 

81 

tical cross sections from Kramers , Experimental iresults on 
helium plasmas that are available at present refer to sufficiently 
high ground state densities and electron temperatures that the 
contribution from the second term of Eq, 2.54 dominates. Experi- 
mental values of Hel population densities from plasmas at very 

9 —3 

low ground state densities e.g. < 10 cm or at extremely low 
electron temperatures (where the first term of Eq, 2,54 dominates) 
can help settle the question of r^ values at low electron densi- 
ties. The comparison of r^^ val\ies of the present work and those 

of Ref. 51 presented in Figures 3.6 show that the r^^ values of 

. ' 51 

the present work agree well with those of Hess and Burrell . 

' ■*!- * 4 * 

The values of r_, r, and r^ for Hell levels were shown in 

Figures 3.9-3.14. For lower levels of Hell, the values of r"*^ 

o 

' 12 "' ' " — 3 '' ■ 

(Figure 3,9) remain constant up to n_ = 10 cm and then rise 
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sharply. In going from lower levels to higher levels the electron 

density at which this sharp rise takes place decreases* The r^ 

values (Figure 3.10) remain almost the same for all the levels of 

13—3 

Hell up to n^ =10 cm , beyond which the Values are lower for 
higher levels. Similarly r^ values (Figure 3.11) are lower for 

higher levels at all electron densities. In contrast to r^ and 

“ 4 " 

r« values, the values of r (Figxire 3.9) are higher for higher 
levels at high electron densities. Thus it can be concluded that 
at higher electron densities, the contribution from the second 
continuum (through r^) daninates over the contribution from the 
groxmd state of He I (through r^) and frcxn the ground state of 
Hell (through r^ ) . 

It may be pointed out here that, strictly speaking, r^ 
values for Hel levels such as those presented in Figures 3. 1-3. 4, 
are functions of n(l), the Hel ground state population density, 
as the heavy particle collisions have been taken into account in 
the model. Contributions to r^ are thus not exclusively from 
the first continuum. The contribution from the ground state of 
Hel is through both r^ and r^^ terms, but the dominant contribution 
is through r^. The dominant contribution to rQ(p) is thus from 
the first continuum. The values of r^ presented in Figxires 3.1- 
3.4 have been computed at low He pressxire corresponding to a^^^ 
ground state population density n(l) = 2.25 x 10 cm . At 
this low pressure , the heavy particle collisicais are not very ^ ^ 

effective and do not contribute significantly to r^ values. Hence 
r^ values given in Figures 3.1— 3,4 may be taken to be independent 
of the gas pressure. 
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Perhaps the best way of testing a collisional -radiative 
model is to apply it to an experiment and to compare the calcu- 
lated and experimentally observed population densities. The 

present collisional radiative model explains quite well the 

34 35 

experimental results of Johnson / Johnson and Hinnov / Otsuka 

33 11 

et al f and Hegde and Ghosh . 

We have made comparisons of our calculated population 

densities with the experimental values of Johnson for eight 

sets of conditions. Results for two sets of conditions y which 

show the best and the worst agreements are shown in Table 4.4. 

In the first set the plasma parameters are kT =7.9 eVy n = 

fc# W 

13 “3 13 —3 

3.3 X 10 cm and n(l) =1,8x10 cm . In the second sety 

the values of plasma parameters are kT^ = 14.4 ev, n^ = 5.1 x 

12-3 11-3 

10 cm and n(l) = 4,1 x 10 cm . The agreement is better 

at low electron temperatiures . In the worst case the population 

densities agree within one order of magnitude. 

Figxire 4.4 presents the population densities of Hel 

35 ■ 

excited states from the experiments of Johnson and Hinnov along 

51 

with the values calculated by Hess and Burrell and also those 

obtained from the present calculations . The plasma parameters 

are kT^ = 1.1 eV, n^ = 2.4 x 10^^ cm~^y and n(l) = 5.75 x 10^^ 

-3 

cm . Major discrepancies (by about a factor of 2) between the 
experimental values and those calculated in this work are observed 
for 3^P and 3^D. The population densities of the remaining states 
are in fair agreement with experimental values y including the sta- 
tes 5^Sy5^Py5^Dy 5^Sy 5^Py sud 5^D f ot which results are not 

51 

available frcxn Hess and Burrell , as they consider separated^^^^^^^^^^^^^ 
sublevels only up to n =4* 
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Table 4,4. Comparison of Johnson experimentally observed 

values of Hel population densities n{p)/g with pre~ 
seitt . calculations ^ 


1 

1 

’l 

1 

States 

kT =7.9 eV 
e 

n^ = 3.3+13* 
e 

n(l) = 1.8+13 

-3 

cm 

-3 

cm 

kT^ =14.4 eV 
e 

n = 5.1+12 
e 

n(l) = 4.1+11 

-3 

cm 

-3 

cm 


—3 

n(p)/gp (cm 

) 

n(p)/gp (cm"^ 

) 

r" 

1 

I 

f 

^ I 

Expt. I Present 

( Ref . 34 ) J Calculations 

1 

Expt. I Present 

(Ref. 34) ; Calculations 

J ^ ^ 

4^S 

4.80+6 

3.27+7 

9.30+5 

1.64+6 

5^S 

1.18+6 

3.33+6 

4.10+5 

1.60+5 

3^ 

1.91+7 

2.23+7 

1.11+6 

3.73+5 

4lp 

4.10+6 

8.17+6 

5.90+5 

2.21+5 

5^P 

1.09+6 

1.43+6 

3.20+5 

9.40+4 

3^D 

1.00+7 

1.95+7 

4,80+5 

3.34+5 

4 D 

3.00+6 

2.62+6 

4,10+5 

7.28+4 

5^D 

7.40+5 

4.12+5 

1.85+5 

2.04+4. 

3^S 

2.60+7 

3.53+7 

2.30+6 

1.22+6 

4^S 

4.40+6 

1.58+7 

7 .00+5 

5.60+5 

5^S 

9.70+5 

1.11+6 

1.76+5 

6 .76+4 

33p 

1.89+7 

2.29+7 

2,20+6 

6.86+5 

4^P 

3.40+6 

1.82+6 

5.00+5 

1.10+5 

5^P 

1.04+6 

3.28+5 

2.10+5 

1.81+4 

3^D 

1.47+7 

7.60+6 

1.11+6 

1.46+5 

4^D 

2.90+6 

1.83+6 

3.60+5 

5.85+4 

5^D 

7.50+5 

3.53+5 

1.22+5 

1.75+4 


* Read 3.3+13 as 3.3 x 10^^. 
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calculated values. 
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The experimentally observed population densities for 

3 1 38 

levels 3 D and 4 D of Hel of Otsuka et al and those calculated 

for their experimental conditions using che present model are 

shown in Table 4,5. The calculated population densities of Ref. 

88 are also presented in this table. In this experiment, n(l) = 

2.25xlO^'^cmr^,the values of electron temperatures at three windows 

are respectively 3.17, 0,88 and 0.34 ev, and those of electron 

14 14 14 —3 

densities are 2.68 x 10 , 1.83 x 10 and 1.06 x 10 cm 

respectively. For window 3, the results of present calculation 

are worse than those reported in Ref. 88, but for windows 4 and 5, 

the agreement between the experimental and calculated population 

densities using the present model is good and better than those 

presented in Ref. 88. 

The emission enhancement factors (R„ = ) for Hel and 

o 

Hell at various magnetic fields based on ratios of population 
densities, as reported in Ref. 11, were also calculated using the 
present model. Whereas the ratios of population densities are not 
sensitive probes to test a model, it was fovind that the theore- 
tical enhancement profiles for all of the Hel lines are in agreemsttt 
within 25% of the experimental profiles. A typical result for the 
transitions 3^D-2^P is shown in Figxare 4.5. The Hell enhancement 
profile, however, was fo\and to be quite sensitive to the ioniza- 
tioi-excitat ion cross section used. This may be expected since 
at zero magnetic field, the dominant term determining the popula- 
tion density of Hell is the ionization-excitation process, whereas 
at all other magnetic fields, the dominant term results fr can the 
population of Hell excited states from tte ground state of Hell. 



Table 4.5. Ccsnparison of experimentally observed populatic«i densities (Ref. 38) for 3^D and 4 
levels of Hel with the present calculations 



lower value to partially optically thick case 
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The calculated population densities of n = 4 level of 
Hell for the conditions of TPD machine are also compared with the 

30 

corresponding experimental results of Otsuka et al and are 
presented in Table 4.6. The following equation has been used for 
calculating the population densities of Hell levels. 


n''’(i) 


r^(i) 


n_ n 


++ 4,14 X 10 


-16 


^72 


exp(Ej, A^T^) 


e 

r, (i) n(l) 4i —7 exp((E^ - E')AT^) + 

^ n-'d) ^ 1 ® 

r 2 (i) 1^ n'^(l) exp((E* - Ep/kT^) (4.2) 

■' - 4 “ + 

where E^ and E^ are respectively the ionizaticai potentials of the 
ith state and ground state of Hell and E^ is the ionization poten- 
tial of the groxind state of Hel to the first continuum. In Table 
4,6, the results of Ref, 88 are also given. For window 3 and for 
the partially optically thick case, the dominant contributicd is 
made by the third term of Eq, (4.2), in contrast to the calculated 
results of Ots\aka et al^^ where the first term is dominant. For 
windows 4 and 5, present calculations and the calculations of 
Otsuka etal^^ agree about the dominance of the first term. The 
n(4) population densities of Hell using the present calculations 
for window 3 are not in as good agreement with experiment as the 
correspcaiding value of otsijka et al^*^, but those for windows 4 and 
5 are in better agreement with the experimental results than those 
of Refs, 10 and 88 . 


Therefore it may be concluded that at low electron temper- 
atures, the results of present calculations agree with the experi- 
mental results but at high electrcaa temperatures the matching is 



Table 4.6. The first# second and third terms of Eq. (4.2) for the i =4 level of Hell 
comparisons of calculated population densities (cm"3) for i = 4 with those 
Refs, 10# 38 and 88 
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not good. The calculated results are in near agreen^nt with the 
experimental results, if lower values of radiative recombination 
coefficients and higher values of election impact excitation rate 
coefficients are _ used, 

4.4. ROLE OF DIRECT IONIZATION -EXCITATION 

Figures 3.15-3.19 show that the direct ionization- 

excitation process plays a significant role in determining the 

population densities of Hell levels, particularly at high electron 

12 

temperatures ( > 2 eV) and at electron densities less than 10 

••3 ^ j. 

cm . Here we comment on the role of the ratio n (l)/n , that of 

n(l) and the contribution of the excited states of Hel. 

The second term of Eq. (3 . 1) presented below, 

n'^(i) = ^i^i^ nH^ * ”e^^^ 

E 

r 2 (i) n ( 1 ) 

n^(l) 

is a f motion of the ratio n^(l )/n'*^'^. It is fomd that for 
kT^ > 5 eV, the variation of this ratio does not alter the popu- 
lation densities of Hell levels more than 2 % in the range n^(l)/ 
-f-j. 2 4 

n = 10 -10 . The values of population densities of Hell levels 

-f- ; ' ' 4.4. 4 

at n (l)/n = 10 are about a factor of 3 larger than those at 

. 4 44 . . 

n (l)/n =1. This is clear from the values of population 

densities of Hell levels shewn in Table 4.7, This table presents 
the population densities of Hell excited states at kT^ = 5, 10, 

15 and 18 eV, n^ = 10^^ cm”^ and for the ratio =1 , 100, 

: 14 :; -3^' ■ 

1000 . The ground state density of Hel is 2.25 x 10 cm and 

. 4 ' ■ : ■ ■ . ’ ■ 

n (1) % n_. The values of population densities of level i = 4 at 

w' 



+ ++ 

Table 4.7. Effect of the variation of ratio (n (l)/n ) on the population densities of Hell 

excited states 
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various electron temperatures and electron densities corresponding 
1*1 

to — j— - = 1000 were already presented in Table 3,6. 
n 

The second term of Eq. (3.1) is a function of also the Hel 

ground state population density n(l),For all the levels / on taking 

lower values of n(l) the value of electron temperature , above 

which direct ionization-excitation becomes important/ increases. 

8 -3 14 

For example/ for i = 4/ at n = 10 cm and n(l) = 1.0 x 10 

e 

-3 

cm / the direct ioniz at ion -excitation is important above kT^ =; 

3.3 eV/ while at the same electron density and for n(l) = 1.0 x 
12-3 

10 cm this process is important above kT^ = 6,6 eV. A 
decrease in Hel ground state population density results in a 
reduced contribution of direct ionization -excitation and in an 
increase in the electron temperature above which this process is 
important. 

Present calculations also indicate that/ for the condi- 
tions of calculations used in the present work/ the major contri- 
bution to population densities of Hell levels is mainly from the 
ground state of Hel. The contribution from the excited states of 
Hel to the population densities of Hell levels is much less. 

This is because of the f act that although the ionization-excitation 
rate coefficients for the metastable states of Hel are much larger 
than the rate coefficient of the ground state of Hel (Table 3.2) 
but the population densities of metastable states of Hel are 
atleast three orders of magnitude smaller than the populaticai 
density of Hel ground state. The contribution from the higher 
states of Hel is even less in populating the excited levels of 
Hell. The present calculations provide an estimate of the role 
of direct ionization -^excitation in populating Hell excited states. 
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Non availability of the exact simultaneous ionization-excitation 

cross sections for the higher excited levels constitutes a major 

shortcoming. However, it was found that the results obtained 

about the range of electron temperat\are where ionization-excitation 

is important are not very sensitive to considerable variations in 

the cross sections. Hence, the results obtained in the present 

work should give a fairly good idea about the quantitative role of 

direct ionization-excitation processes in spite of xancertainties 

in the cross sections. Thus, it may be concluded that direct 

ionization-excitation process dominantly populates Hell levels at 

12-3 

kT_ > 2 eV and n <10 cm and must be an integral part, of the 
model under such conditions. 

4.5. MECHANISM OP POPULATION OF EXCITED STATES 

The procedure outlined in Chapter 2 and the results 
presented in Chapter 3 for evolving the mechanism of population of 
excited states give an idea about how the various excited states 
of Hel and Hell are populated and which rate coefficients and 
transition probabilities have major significance. The present 
work, perhaps for the first time, focusses a method which systema- 
tically examines the importance of individual atcmic processes. 

It is true that the dominant mechanism that finally evolves ulti- 
mately depends on the choice of cross sections; however, for any 
choice of cross sections, the above procedure can be used, which 
is of general applicability. 

It may be added here that whereas the above numerical 
proced’ure can be used for identification of the important elemen- 
tary processes, an idea about the dominant processes of the 
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mechanism in populating various levels can be had from an analy- 
tical consideration of the form and constituents of the elements 
of the inverted matrix. This can be done in the following manner. 
The elements of a® inverted matrix X (X = ) pertaining to a 

collisional-radiative model are constructed out of the various 
elements a^^j of the original matrix K, which themselves are 
ccxribinations of rate coefficients of elementary processes as 
described in Figiare 2.2. A peculiar property of the inverted 
matrix is that an element x. . of the inverted matrix X of dimension 
m X m can be written in terms of the corresponding element x^j of 

inverted matrix X' of dimension (m-1) x (m-1) as x. . = xf . + an 

J J 

additional term- For example, the terms of inverted matrices 
of dimensions 2x2, 3x3, and 4x4 (obtained by using SNOBOL 
programme; see discussion later) are as follows: 

2x2 


^ 21*^12 


(4.3) 


3x3 


= X, 


2x2 ^131*^351 


(4.4) 


4 X 4 


^ 131 *^ 41 . 

^31 ' ^ 

^441 


^43i'^351 

^31 

^ 431 *^341 


(4.5) 
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Population of a level is given by the product of matrix 
multiplication of a given rcw such as i of the inveirted matrix 
with Goiiartm matrices B, C, D, and E (as shown in Figure 2 .2) . 

That particular term x. . which leads to the most significant value 
of the product ma]<es the dcminant contributicai to the mechanism. 

A detailed analysis of the term structure of an inveirted 4x4 
matrix (along the lines of Eq. (4.5)) and a comparison with the 
ntmerical valvies of the elements of 45 x 45 inverted matrices (for 
various values of plasma parameters in the collisional-radiative 
model) indicate that^ for the range of values of rate coefficients 
in the present work^ the first term of all the elements of the 
inverted matrix is the dominant one. Therefore/ all those 
processes whose rate coefficients make significant contribution 
to the matrix elements (of the original matrix A) which 

appear in the first term of x^j would be impoirtant. 

The process of sorting out the important terms is perhaps 
more conveniently accomplished by having the inverted matrix in 
algebraic form. For this a programme written in SNOBOL was used 
(see Appendix III), where the output provides the inverse of a 
matrix in the algebraic mode as is given in Eqs. (4.3-4.5). From 
this one could directly estimate the role of a particular rate 
coefficient in populating or depopulating a particular level. 
Without repetitive inversion of the original matrix, the method 
using SNOBOL permits direct estimates of contributions of various 
processes, but the ccxnputer memory requirement to completely write 
out the inverse of a 45 x 45 matrix in algebraic mode is substan- 
tial. The method however is potent as it permits a direct insight 
into the complicated processes. 
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From this study in simplification of the model we conclvide 
that more than 70% of the processes used in the model, could be 
neglected and yet population densities could be obtained within 
20% of those from the elaborate models. 

This method focusses attention on those atomic processes 
which make substantial contribution to the population densities, 
hence identifies the processes where knowledge of experimental 
cross sections would be of great value, if they are not already 
available. 

4.6. lASER INDUCED SELECTIVE EXCITATION 

In this section we critically examine the results of seme 

recent laser induced selective excitation experiments against the 

predictions of the collisional-radiative model of the present 

39 

work. We consider the experiments of Catherinot et al , Yasvimaru 

42 90 

et al , and Gauthier et al . 

39 

In the two sets of experiments of Catherinot et al , 

two different transitions are pumped. In the first set, 2^S-3^P 

transition (5016 A) is pumped by a pulsed laser (500 Watts/pulse) 

at P = 1 torr, T = 400°K, n_ = 10^^ cm”^, kT. = 5.17 eV, laser 

^ © © 

O 

pulse base width 5 nsec, spectral width 0.1 A, and a pulse time 

delay of 2.1 nsec. The calculated enhanced population densities (in 

1 1 3 3 

terms of arbitrary \anits) for 3 P, 3 D, 3 D and 3 P states vinder 
the above conditions along with the experimental results are 
presented in Pig\are 4.6. In the second set, the transition 2^P-4^D 
(4922 A) is pumped with a laser pulse, time delay of 11.2 nsec. 
Other parameters are the same as for the first set. Experimental 
results for 3^P, 4^P and 4^D levels along with the calculated 


results are shown in Pigxire 4.7. 
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Time (ns) 

Figure 4.6, finhancement and decay of population densities of 

'1 3 3 39 ■ 

3 D, 3 D. and 3 P states r experiment ? ■ 

calculated, partially optically thick and low KN(p>q) 

values; calculated, partially optically thick 

79 ’ ' 

and high i<M(p,q) valiies ; calculated, optically 

values; calculated higher 

^ ^ optical escape factors and lew KN(p,q) valu 


Enhanced Population Densities (arb. units) 



Fig\ire 4.7, Enhancement and decay of populations of 4^P and 

■ ' ' 3 39 ’ 

4 D states j , experiment ; calculated, 

partially optically thick and low KN(p,q) values; 

-calculated, partially optically thick and higl 

79 ^ , 

KN(p,q> values , ■ ■ 
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Figure 4.8 presents the calculated time dependent popula- 
tion densities compared with the corresponding results of Yasumaru 

■ -42 42 1 1 

et al . In the experiment of Yasumaru et al , the 2 S-3 P 

transitican (5016 A) is pumped with a laser pulse of width (FWHM) 5 

nsec, spectral width 0.1 A, and relative population densities of 

the 3^P and 3^D levels are monitored at two pressures, 0.4 torr 

and 5 torr. The electron densities amd electron temperatures at 

10-3 

these two pressures respectively are n^ = 2 x 10 cm and 3 x 

11 --S 91 

10 cm , and kT^ = 10.34 eV and 3,45 eV. The calculated and 

experimental population densities at both the press\ares are shown 

in Figure 4,8. 

90 

In the experiment of Gauthier et al , the triplet transi- 
3 3 ° 

tion 2 S-3 P (3888.6 A) is pumped by a laser pulse of power 2 KW, 

■ o 

width (FWHM) 3. 5 nsec, and spectral width 0.1 A. The plasma para- 
meters are T^ = 600 °K, n^ = 2,9 x 10^^ cm“^, T^ = 300 °K and the 

17 -3 

ground state population density = 4.6 x 10 cm . The fluore- 

3 3 3 

scence decay of 3 P, 3 S and 3 D levels are monitored as a function 
of time. In this experiment the electron temperature (kT =0.0517 
eV) is q\xite low. At kT^ 0.1 eV, the population densities from 
the present collisional-radiative model are rather uncertain 
because of s\ibstantial xander flow in the process of computer 
calculations. Hence, in order to make a rough comparison, the 
absolute population densities as a function of time are calculated 
at higher electron temperatxires kT_ = 0.517 eV and also at kT^ = 1 
eV and the results are presented in Figure 4,9. 

In none of the above laser induced selective excitation 
experiments, the absolute population densities are provided. 

Hence to compare the calculated absolute population densities^ ^ 



Popylo t f on Density ( retot ive units ) 
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Figure 4,8. 


Relative population densities of 3^P anu 3^0 states 
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calculated, only one rate coefficient (3^P-3^D) non 
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with the experimental results^ the calculated absolute population 
densities are converted to relative population densities by normal- 
izing with respect to one of the experimental population densities. 
In Figures 4.6 and 4.7 enhanced population densities N (N = 

normalized with respect to respectively S’*"? and 4 D peaks 
are presented. In Figxare 4.8/ relative population densities 

normalized with respect to 3^P peak at 5 torr are shown. In 

. —3—1 

Figiore 4 . 9 , the nimber of photcxis emitted cm sec (calculated 

absolute population density x transition probability/ normalized 

3 

with respect to 3 P peak) against the elapsed time are plotted. 

The calculated absolute population densities corresponding to the 
data presented in Figures 4. 6-4. 9 are presented in Appendix II. 

85 

The derivation of Einstein coefficients B involves the 
use of Boltzmann distribution law. ThuS/ strictly speaking/ the 
term B.p(v) gives the rate of induced emission or induced absor- 
ption for equilibrium situations. However/ even for the non 

51 

equilibrium situations , it has been used in the recent literature . 
It has been found that a considerable variation in the numerical 
values of B (by an order of magnitude) does not change the rela- 
tive population densities presented in Figures 4.6-4.9/ and also / 
the peaks of enhancement curves on the time axis do not shift much. 
Cnly the absolute population densities increase or decrease with 
an increase or decrease in the numerical values of B. 

A comparison of the magnitudes of electron impact and 
heavy particle collision terms in the rate expression of excited 
state population densities (Eq. 2.46) shows that in the presstire 
range of the laser induced selective excitation experiments 
disctissed above/ heavy particle collisions and the radiative 
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processes, and not the electron impact processes, play dominant 
role in populating or depopulating the levels studied. This is 
presented in Table 4.8 where the numerical values of all the terms 
of rate equation (2. 46) are shown for the experimental conditions of 

O Q 19 

the experiments of Catherinot et al (kT^ = 5.17 eV, n^ = 10 

0 . 0 

—3 42 

cm and P = 1 torr) and of Yasumaru et al (kT^ = 3.45 eV, = 

11-3 

3 X 10 cm and P = 5 torr) for some levels of Hel. 

The calculations on temporal enhancement and decay profiles 

of population densities show that with rate coefficients for 

neutral-neutral excitation, deexcitation, ionization as well as 

79 

recombination taken from Drawin and Emard , the calculated rela- 
tive enhancements with respect to the pumped levels for the condi- 

39 

tions of the experiment of Catherinot et al and of Yasumaru et 
42 

al are much higher (shown by dotted lines in Figures 4.6, 4i*7 

and 4.8), than the experimental values. They also do not agree 

with the experimental peaks in terms of their positions on the 

time scale. It was found that the rate coefficients for heavy 

79 

particle collisions given by Drawin and Emard are about 1-2 
orders of magnitude higher than that needed for an agreement with 
the experimental results and the reduction of the factor 128 (in 
Eq. 2.24) to 2.3 helps. Results of calculations with the factor 
128 in Eq. (2.24) (Eq. 23 of Ref. 79) reduced to 2.3 are shown by 
dashed lines in Figiores 4. 6-4. 8. A definite improvement in the 
calculated enhancements is noted. With this modified expressicai 
for heavy pairticle collisional rate coefficients, the deexcitation 

rate coefficient KN(3^P-3^D) turns out to be 1.96 x lO"^*^ cm^ 

-1 39 

sec for the conditions of the experiment of Catherinot etal . 

This reduction of heavy particle rate coefficients, however, is 
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still not adequate to reproduce the decay patterns. The calcu- 
lated decay curves are less steep and this departure may be attri- 
buted to the low values of the optical escape factors , used in 
this work. This is so because the curves represented by dotted 
lines and dashed lines have been obtained by using partially 
optically thick conditions and are broader than the experimental 
cvirves. CSi calculating the time dependent population densities 
under optically thin conditions (shown by dot-and-dash lines in 
Figure 4.6), the enhancement cxxrves are found to be very sharp, 
in fact sharper than the experimental curves. Thus the optical 
escape factors used in this work are smaller than what are required 

to explain the experimental enhancement curves of Catherinot et 
39 

al . The optical escape factors used in the present calculations 
are calculated by using Bq. (2.32), This equation yields for the 

_3 

conditions of the experiment of Ref. 39 A . . = 3.32 x 10 

3 P-l S 

By using higher (about one order of magnitude) OEF's, sharp decays of 

all the foxir curves (shown by dash-double dot lines in Figtare 4.6), 

in near agreement with the experimental results, can be obtained. 

It may be added here that the optical escape factor for 3^p-l^S 

41 

transitictti experimentally estimated by Dubreuil and Catherinot 

10 n —3 

at 0,4 torr, in the n^ range 10 -10 cm and T^ = 350®K, is 

A , . = 0,014. The estimated values of a - . and t ^ . 

3 p-l S _2 3-^P-l-^S 3 P-l-'-S 

in Ref, 41 are respectively 2.6 x lO” and 106. Under similar 

■ *-5 ' ** ' 

conditions Eq. (2.32)yields a . . = 1.99 x 10 , t , = 

3 P-l-^S 3 P-l-^S 

107.41 and A , = 2.75 x 10”^. Obviously, the experiment of 

3 P-l-^S 

Catherinot et al^^ requires optical escape factors about one order 
of magnitude higher than those given by Eq, (2.32), 
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In Fig\ire 4.8# using the reduced rate coefficients for 
heavy particle collisions# and taking the 3^P peak at 5 torr as 

the normalized reference point, the calculated results are compared 

. 42 

with the experimental results of Yasumaru et al both at 0.4 torr 

and 5 torr. Under the conditions of this experiment# for the 

3^P-3^D excitation transfer# the value for the electron impact 

rate coefficient K(3^P-3^D) used in the model is 2,19 x lO”^ cm^ 

sec ^ and for the heavy particle rate coefficient KN(3^P-3^D) is 
—10 3 -1 

1.45 X 10 cm sec # both of which are close to the correspon- 

42 

ding experimentally estimated values of Yasumaru et al , Their 

-5 -10 

corresponding values respectively are 3.5 x 10 and 3.5 x 10 
3 -1 

cm sec .In order to explore the role of heavy particle rate 

—10 3 —1 1 1 

coefficients# the value 3.5 x 10 cm sec for KN(3 P-3 D) 

. 42 

given by Yasumaru et al along with all other KN's equated to 
zero was used. The resulting enhancement curves for 3^P and 3^D 
at 5 torr are shown by dot -dash lines in Fig\ire 4.8. The ctirves 
thus obtained are broader than the experimental ones as well as 
the curves calculated earlier which considered additionally other 
transfer processes. This indicates that in the decay of the 
population densities of enhanced levels# several states make 
significant contribution. 

It has already been mentioned that since the electron 

90 

temperature in the experiment of Gauthier et al is very low 
(0.0517 eV) # the time dependent population densities were calcu- 
lated at two higher electron temperatures 0.517 eV and 1 eV to 
have a rough idea of the matching of calculated and experimental 
enhancement curves# as the present calculations are rather 
uncertain at low electron temperatures (kT^ < 0.1 eV) due to 
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substantial londerflow. Since lander the conditions of this experi- 
ment the contribution of electron impact processes is relatively 
\animportant ^ it is expected that the results at relatively higher 
kT^ should not deviate much from those calculated at lower electron 

3 

temperatures. Relative intensities for 3 S and the decay patterns 
3 3 

for both 3 P and 3 D are significantly off. A change in the 
optical escape factors in the singlet system has imperceptible 
effect on the calculated curves . Incorporation of partially 
optically thick condition for the 2^S-3^Py 2^S-4^P and 2^S-5^P 
transitions increases (not shown in Figure 4.9) the relative 

3 

fluorescence intensities of both 3 S (by a factor of 2) and 

3 

3D (by about an order of magnitude)^ keeping the peak shapes 

3 

practically imchanged but flattening out the decay of 3 P, 

Triplet states with n = 3 do contribute to the associative ioniz- 
ation process^ a process which is not included in the present 
model but is expected to affect materially the triplet state 
population densities at high press\ares. This could be a reason 
of the discrepancy between experimental and calculated cxirves of 

Figxjre 4.9. The required increase in calculated relative popul- 

3 3 3 

aticffi density of 3 S (in canparison with 3 P and 3D in Figure 

4.9) is ccxnpatible with the estimates of associative ionization 

3 3 ^ ^ 

cross sections in the literature: (3 S) = 0.06/ Q^(3 P) - 

Q^(3^D) = 2.9 (Ref. 92)/ and o(3^S) = 0/ a (3^) = 1.6/ o{3^D) = 

4.5 (Ref. 93)/ all in vinits of 10"^^ cm^. 

Thus from the application of the model to a range of laser 
excitation studies in he lixam plasmas / it may be eoncluded that the 
heavy particle rate coefficients from Ref. 79 used in this work 
should be lower by about 1-2 ordersof magnitude and optical escape 
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factors from Ref, 83 should be larger by about one order of 
magnitude. 

This latter conclusion if applied generally to Eq. (2.32) 

(Ref. 83) would cOTie in conflict with its suitability already 

noted in connection with the applications to steady state he lixjm 
11 34 3 5 38 

plasmas. * * * The latter experiments were at low pressures 

—2 0 

with optical escape factors in the range 10 -10 . The implication 

here may be that whereas Eq, (2,32) is satisfactory in the low 
pressure range/ at high pressures/ the relative accuracy of the 
optical escape factors from Eq. (2.32) (where absolute values are 
small) is poor and need of appropriate corrective factors is 
indicated. 

4.7. DETERMINATION OP ELECTRON TEMPERATURE FROM LINE INTENSITY 
RATIO 

Both spectroscopic methods and electric probes have been 
used for determination of electron temperatvire T^ in plasmas. 
Spectroscopic methods have an advantage over electric probes 
because T^ estimates can be obtained without disturbing the plasma. 
Among the spectroscopic methods the one using line intensity 
ratios which has bearing in the collisional-radiative model has 
been controversial. Since the present work involves a 

detailed examination of the collisional-radiative model of he lixam 
plasmas we examine here/ using the results of our calculations as 
well as some experimental data/ reliability of the method, 

Cunningham^^ suggested that in non-thermal plasmas the 
electron temperature could be determined from the ratio of two 
spectral lines if they have similar excitation energies and the 
excitation cross sections of the emitting levels have very 
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different dependence on electron energy. These criteria aire inet 
in the case of singlet and triplet lines of Hel. Cunningham 
applied this method to Hel lines 4^S-2^'P (4713 A) and 4^D-2^P 
(4921 A) . 

97 

Drawin and Henning also used the saine pair of lines. 

96 

However^ in addition to processes considered by Cunningham they 
considered ionization from the excited levels 4^D and 4^S and 
collisions of excited atoms with other plasma particles/ parti- 
cularly the sing let -triplet excitation transfer. They found that 

the line intensity ratio shows less temperature dependence than 

96 

that observed by Cunningham / and the intensity ratio depends 

more strongly on gas pressure than on electron temperatxxre . They 

concluded that this method is suitable for determination of 

electron temperature in non-thermal laboratory plasmas only if 

the lines lie far from one another. 

98 

Sovie / on the basis of a simplified collisional- 
radiative model/ calculated the line intensity ratio of Hel lines 
4^s-2^p (5047 A) and 4^S-2^P (4713 h) . For plasmas with n^/n(l) < 
10 / he considered excitations from the groxand state and meta- 
stable states to excited levels and their radiative decays. The 
effect of cascading was also considered. He used the following 
relation to calculate the line intensity ratio: 

X = - = [k(1^S-4^S) + B*K(2^S-4^S) + D*K(2^S-4^S) ^ 

° ^4713 . 

+ I R(r^P) A(r^P-4^S)] A(4^S-2^P)/ 

n(l s)n r>4 
e ■ — 

[K(1^S-4^S) + B'K(2^S-4^S) + D'K(2^S-4^S) 

+ . — J R(r^p) A(r^P-4^S)] A(4^S-2^P) 
n(l'^S)n r^4 


(4.6) 
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where K represents electron impact excitation rate coefficients; 

n(l^S) and n^ represent gro\ind state population density and 

electron density respectively. B* and D respectively are n(2^S)/ 

1 3 1 — 

n(l S) and n(2 S)/n(l S), andA(a-b) = A(a-b)/A(a, tot) where A is 

the spontaneous transition probability. R(a) is the rate of 

production of state a. The ratio given by the above equation is a 

13 

function of n(2 S) and n(2 S). To accoxint for the cascading 
98 

effects y Sovie suggested that the quantity to be experimentally 

determined can be obtained by subtracting the product of observed 
11 — 1 1 

r P-4 S transitions and A(4 S-2 P) from the observed number of 

o 2 1 

5047 A (4 S-2 P) transitions. On the basis of this premise and 
Eq. (4.6 ),one obtains 

X = ^K(1^S-4^S) + B'K(2^S-4j‘'S) + D‘K(2^S-4^S) ] A(4^S-2^P) 

[k(1^S-4^S) + B'K(2^S-4^S) + D'k(2^S-4^S) ] A(4^S-2^P) 

(4.7) 

Frcxn Eqs. (4.6) and (4 ,7 )y Sovie evaluated the valxoes of X_ and X _ 
for kT^ range 4-20 eVy both with and wi^-hout the contribution from 
the metastable states . The ratios were found to be dependent lonly 
on T^. Howevery this method is valid only, when population of 
excited states by radiative recombination and depopulation by 

electronic collisions are unimportant processes. 

99 ' ' ' ' 

Schieber et al found that by using different cross 

sectionsy different estimates of T_ are obtained and concluded 

e 

that with the available knowledge of cross sections y determination 
of electron temperatiore by this method is of little value . ^ ^ ^ ^ 
Eastl\and et al^^^ modified this method by suggesting an averaging 
technique to evaluate experimental data. This technique uses both 
relative and absolute line intensities. Howevery they do not take 
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into account the reabsorption of resonance radiation y excitation 
transfer and cascade processes. 

101 

In a recent work, Brenning has examined the limitations 

of the method of determination of fron line intensity ratios y 

for various pair of lines. He foxind thaty in general y secondary 

processes such as excitation from metastables, excitation transfer 

collisions, cascading etc., greatly influence the method and the 

7 -3 

method is restricted to only low density (n^ < 2 x 10 cm ), 

short duration (t < 5 x lO”^ sec) plasmas. He suggested a 

method for T determination from absolute line intensities, 
e 

Thus the usefulness of this method has been viewed 

96 98 

diversely. Whereas C'unningham and Sovie support this method 

97 

with certain qualifications, Drawin and Henning , Schieber et 

99 100 101 

al y Eastlund et al and Brenning either reject it or 

recommend it for a very restricted use. We examine here the 

validity of this method for determining T , by using the colli- . 

sional-radiative model calculations of the present work* 

Using Sovie 's simplified model (Eqs. 4.6 and 4.7) and 

the rate coefficients used in the present model, we calculated the 

ratio both with and without metastables, at two 

12 14 -3 

electron densities n^ = 10 and 10 cm . In this, for the 

electron impact excitaticxi l^S-4^Sy we used, in two separate set 

68 

of calculations, cross secticxis given by Raan et al and by ^ ^ ^ ^ ^ 
65 

Drawin . The results for partially optically thick case are 
presented in Figure 4.10, We find that the line intensity ratio 
^5047'^^4713 strongly depends on the cross section for the process 
1^S-4^S , On uising Raan et al's^® experimental cross section, the 



5047 / ^4713 (Eq. 4 * 6 ) 


zo 


18 



Flgtire 4.X0. Variation of calculated (Sovle's model®® and preaent rate 

coe££iclenta) relative nximber of 5047 a and 4713A transitions 
12-' 14 -3 

with at » 10 and 10 cm for partially ofticaliy 
thick case y . eatperlmantal cross section^®, metastables 

included; — , experiiirental cross section**®, metastables 

6 IS" ' ' ' ' 

exclutted; empirical cross section , metastables inclu- 
ded; empirical cross section^^^ metastables excluded. 
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intensity ratios first decrease and then gradually increase. 
Whereas on using Drawin's empirical cross section, the curves 
obtained are more steep. Whether Eq. ( 4 . 6 ) or Eq. (4,7) is used 
the results do not vary significantly, on excluding the excita- 
tions from metastable states the line intensity ratios regularly 
in^trease with increase in and are independent of n^, but still 
they depend on the choice of cross sections. With metastables 
included the line ratios are also a fianction of electron density 
as is evident from Figure 4.10. Fig\are 4,11 presents the results 
of ratio ^ 5047/^4713 (obtained by Sovie's method) for an optically 
thin plasma. By comparing the results obtained with partially 
optically thick conditions (Fig\ire 4.10) and optically thin 
conditions (Figiare 4,11), we find that absorption of resonance 
radiation also affects the line intensity ratios. Similar results 


are obtained for the ratio ^72 qi/^ 7065 Presented in Figure 4.12. 
Thus these results show that Sovie's line intensity ratio method 


(neglecting secondary processes) is not valid for determination 


of T 


Sovie's method does not take into accoxint radiative 


recOTibination , electron impact deexcitation, electron irrpact 

icmization and recombination, excitation transfers, and absorption 

of radiations by plasma, and therefore it does not describe a 

realistic situation of non-LTE laboratory plasmas. Using Sovie's 

model and oior rate coefficients we also calculated the absolute 

3 13 1 ^ ^ ^ ^ 

population densities of the levels 3 S, 3 S, 4 S and 4 S for the 

34 

conditions of the experiments of Johnson , and Johnscxi and 
35 

Hinnov . Table 4.9 compares the population densities obtained by 
Sovie’s method with the experimental results. Population densities 



kTe («V) — — ► 

' 1 ^' ' ■ ■ 

ri^ur* 4 . 11 . Variation of o«lcalat«d (using sovltt'i Modal fOKl jpMMnt 

' '121 ’ - a . ' 

rata coafflelawts) ^5047^4713 with at • 10 ■ 

for optically thin casa; , axperinental cscoas aaction^^, 

Mtaatablas includadi — aKparisental cross aactions . 

BWtastablas axcludad# *«.•,-# aMpdrlcal cross oaction*®, nwta-> 

■ ■ 65 

stablas includadi satirical cross sactioo , aetastahlas 

• awJludad. 
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flowm 4 . 12 . VmrlmtJLvu ei «»toa«t4id (So«l«*a and pra —o t rat* 

coaffiaianta) zalatiaa iraabar «E 728lA and Toa^ ixanaltlana 
12 •» 

at « 10 CM and £or partially optically tMck caaat 
iMtaatablas Ineladafli nataatablaa aoDeludad..' 


Coscode corrected ratio 
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obtained from the present model are also shown In Table 4.9. 

From this table ^ it is clear that Sovie's method explains the 
experimental population densities of only 3^S level. The popu- 
lation densities of other levels are very far from experimental 

1 ■ 

results. 

The collisional-radiative model used in the present work 

takes into consideration deexcitation by electron impacts electron 

impact ionization and recombination ^ radiative decays excitation 

transfer between excited states and absorption of radiation by 

plasma. From this model, taking all these processes into accoxint, 

we have calculated the line intensity ratio ^ 5047 /^ 4713 * line 

intensity is taken to be the product of the calculated absolute 

population density and spontaneous transition probability of the 

upper level (e.g. I 5047 = n(4^S) . A(4^S-2V))* The plasma 

parameters used in these calculations are kT = 1-20 eV, n_ = 

10^-10^^ cm ~ 300°K, n(l) = lO^^-lO^^ cm”^. Both optically 

thin and partially optically thick cases were considered. Figure 

4.13 presents the calculated line intensity ratios for 3 ^ 5047/^4713 

against electron temperature. The cxarves represented by conti- 

14 -3 

nuous lines correspond to n(l) = 1.28 x 10 cm and a partially 
optically thick case. 

9 -3 

At low electron densities (_< 10 cm the intensity ^ 
ratio is a f unction of T^ alone but at higher electron density it 
becomes a function of both T^ and n^. As the electron density 
increases the intensity' ratio in general decreases. For a parti- 
cular electron density, the variation of line ratio with electron 
temperatiire is irregular. Figure 4.13 also presents line inten- 
sity ratio data for optically thin case (shown by dashed lines) 



5047 'H713 



9 11 

kTe (eV) - 


figur* 4.13, 


Variation o£ tlw calcul*t«<3 lla« iSikm»S*-y ratio ^ 5047/^4713 

with alactron tanparatara at •• « partially 

optically thick, nd) • 1.28 X 10^* c»“^r — optically thin. 


optically thick. 0(1) • 1.28 X 10^^ . optically 
n(l) - 1.28 X 10^* ca*^! partially optically thick, 
n{l) • 1.0 X 10^* 
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and for high pressure n(l) = 1.0 x 10^^ cm”^ ( shown by dash-dot 
lines) case at n^ = 10 ^^ ^^-3 ^ ^^14 ^^-3; ^ comparison 

of these curves with curves mentioned earlier reveals that line 
intercity ratios, in addition to n and T , are also functions 

G 0 

of gas pressure and extent of absorption of radiation. 

Thus, for determination of frcan the line intensity 
ratio method, a knowledge of n^, gas pressxire and optical thick- 
ness is essential. Moreover, a particular valve of line ratio 
at a particular n^, may correspond to two or more values of 
electron temperature. Thvis the line ratio method, even with a 
knowledge of n^, n(l), and optical thickness may not give 

correct values of T^. However, at low electron densities 
9 **3 

(n^ <10 cm ), can be estimated from this method in the 
kT^ range 2-12 eV. 




CHAPTER 5 


CONCLUSION 

The pxirpose of the present work was to improve the 
collisional-radiatlve model for low temperatiire non-LTE helium 
plasmas and gain an insight into the complicated coupled processes 
which contribute to observed phenomena. This was attempted through 
a ctetailed examination of the atonaic processes in the model and 
also comparisons of predicted population densities with experimen- 
tal observations. It is believed that the effort, within the 
inherent limitations of collisional-radiative models, has yielded 
some positive results. 

Introduction of experimental cross sections for a nxamber 
of electron impact processes, used for the first time in colli- 
sional-radiative models for helium, is believed to have signifi- 
cantly contributed to the near agreement with the results of 

34 

several Hel, Hell experiments, viz., those of Johnson , Johnson 

35 38 11 

and Hinnov , Otsuka et al , Hegde and Ghosh , involving popu- 
lation densities of 3^P, 3^D, 3^S, 3^P, 3^D, 4^S, 4^P, 4^D, 4^S, 
4^P, 4^D, 5^S, 5^P, 5^D, 5^S, 5^P, 5^D states. The population 

coefficients r^, r^^ were coit^red with recently available values 
52 51 

of Fujimoto and Hess and Burrell . r^^ values are in close 

agreement; the source of differences in r^ at low electron densi- 
ties could be settled if experimental population densities are 

9 —3 

available at low ground state densities ( < 10 cm ) or at 
extremely low electron temperatures. The relative lack of agree- 
ment at high electron densities indicate need of lower values of 



148 


ireconibination coefficients and higher values of electron 
impact excitation rate coefficients. A detailed examination of 
the direct ionization-excitation process shows that it is domi- 
nant in populating Hell levels at kT >2 ev and n < 10^^ cm“^ . 

e e 

The process is more effective in populating Hell excited states 

at high gas pressures. Fran a thorough examination of the line 

ratio ^method of determination of electron temperature/ we find 

that in the n range 10^-10^^ cm”^, kT range 1-20 eV the line 

ratios are not independent of n^ and n(l), and we conclude that 

9 —3 

the method is not dependable at n^ 

From applications of the present collisional-radiative 

moctel to recent laser induced selective excitation experiments 

39 42 90 

of Cathetinot et al / Yasumaru et al and Gauthier et al / 

we find that a good agreement between the experimental enhance- 
ments and decay patterns with the calculated ones can be obtained 

only if the heavy particle collisional rate coefficients from 

79 A 

Drawin and Emard are by about 1-2 orders of magnitude 

83 

and the optical escape factors given by Drawin and Emard are 

higher by about one order of magnitude . The latter conclu- 
sion and the fact that optical escape factors from Ref. 83 yield 

satisfactory population densities in steady state helium ^ 

11 34 35 38 

plasmas ' ■' ' might appear to be conflicting. The latter 

experiments/ however/ were at low press\ires with optical escape 

—2 0 

factors in the range 10 -10 . The Implication here may be that 

whereas Bq. (2.32) yields satisfactory optical escape factors in 
the low pressure range/ at high pressures/ where absolute values 
of optical escape factors are small/ the relative accuracy of 
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the optical escape factors from Eq. (2.32) is poor and there is 

need of appropriate corrective factors . 

An interesting study in the present work has been the 

58 

evolution of a systematic procedvire to sort out the important 

rate processes used in a collisional-radiative model. Though 

the method has been applied to heliiam plasmas, it is of general 

applicability. It is shown that by application of this method, 

it is possible to leave aside as much as 70% of the atomic 

processes employed in elaborate models and still predict, in the 

kT^ range of 4-20 eV and n^ range 10^-10^^ cm”^. He I population 

5S 

densities within 20% of the predictions of elaborate model. 

The method yields statewise population mechanism and also 
focusses those processes where knowledge of experimental cross 
sections would be particularly valuable in enhancing the 
predictive properties of a collisional-radiative model. 
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APPENDIX I 

POPULATION DENSITIES WITH CC^RECTED WEIGHT FACTORS 

—3 

Population densities (cm ) of Hel levels with the weight 

2 2 

facstors used in Ref. 56 (g = 4n for n < 6.^ g = 8n for n > 6) 

n n 

2 

and with the corrected weight factors (Table 2.1, g^ = 4n ) are 

presented in the Table -I.l. Plasma parameters are kT_ =10 eV, 
12 —3 14 -3 

n =10 cm , n(l) = 2i0 X 10 cm , T = 300 °K and optically 
0 g 

thin case. 


f 
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Table I.l. Population densities with weight factors used in 
Ref. 56 and with corrected weight factors 


Levels 

Population densi- 
ties with weight 
factors used in 

Ref. 56 

Population densi- 
ties with correc- 
ted weight factors 

- 

Percent change in 
population densi- 
ties 

2^S 

1.472+11 

1.472+11 

0 

2^P 

6.754+07 

6.754+07 

0 

3^P 

6.901+08 

6.905+08 

0.06 

3^D 

1.243+08 

1.243+08 

0- 

3^P 

5.389+07 

5.390+07 

0.02 

4^S 

1.770408 

1.770+08 

0 

4 D 

3.471407 

3.467+07 

-0.1 

4^P 

3.067407 

3.068+07 

0.03 

5^3 

4 .088+07 

4.110+07 

0.5 

5^S 

2.591407 

2.598+07 

0.3 

S^P 

3,146+07 

3.23540'7 

2.8 

5 D 

4.245+07 

4.366+07 

2.9 

5^D 

1.548+07 

1.590+07 

2.7 

5^P 

2.186+07 

2.194+07 

0.4 

n=6 

4.^10 i07 

4.331+07 

> -1.8 

7 

2.577407 

2.554+07 

-0.9 

8 

1.851407 

1.839+07 

•»0#6 

9 

1.587407 

1.579+07 


10 

1.552407 

1.545+07 


11 

1.698+07 

1.691+07 


12 

1.933+07 

1.924+07 
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APPENDIX II 

ABSOLUTE POPULATION DENSITIES FOR FIGURES 4.6-4. 9 

In Figiires 4, 6-4, 9 y the noirmalized relative population 
densities were presented. In the following tables y we present 
absolute population densities corresponding to the data given in 
these figures. In these tables y for a particular curve and for a 
particular level y the absolute population densities have been 
listed for the initial time (t = Oy first value in each box) y for 
peak (second value in each boXy time in nsec mentioned in paren- 
thesis )y and for t = 40 (75 in case of Figure 4,9) nsec, third 
value in each box). 


—3 

Table II, 1, Absolute population densities (cm ) at t = 0, peak, 
and t = 40 nsec corresponding to data presented in 
Figure 4,6 


1 

t 

Curve } 

Time 

I 

1 Absolute popialation densities 

1 ^ ^ ^ ^ ^ ^ 

(cm“^) 

i 

1 

1 

1 

(nsec) 

* 1 
{Level 3 P 

I 

1 • 1 
; 3D 

...1 .. ... 

* 3 

; 3D 

1 _. .' 

* 3 

{ 

. -I . 

dashed 

line 

0 

2.8+10* 

3.3+10 

1.4+10 

3.4+10 

peak 

4.1+12(4.4) 

6.1+11(17.2) 

1.9+11(18.6) 

2.0+11(35.2) 

40 

6.9+11 

3.7+11 

1.3+11 

2,0+11 

dotted 

line 

0 

6,0409 

1.6+10 

3,7+10 

3.1+10 

peRk 

1.2+12(3.2) 

1.5+12(4.4) 

2.9+12(5.0) 

1.1+11(23.6) 

40 

4,9+10 

1.2+11 

3.5+11 

1.0+11 

dash- 

0 

8.7+08 

3.0+09 

4.1409 

1.6+10 

dotted 

peak 

4.1+11(4.2) 

2.8+10(7.8) 

1.2+10(7.8) 

2.1+10(9.0) 

line 

40 

9.1408 

3.3+09 

5.1+09 

1.9+10 


* The absolute population densities of level 3^P at time t =0, at 


peak (which occtirs at 4.4 nsec), and at t — 40 nsec respectively 
are 2,8 x 10^^, 4,1 x 10^^, and 6,9 x 10^^ cm”^. 
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Table 11.2. 


Cxxrve 

T 

i„ , , 

{ Time 
! (nsec) | 

JL i 

dashed 

0 

line 

peak 

40 

dotted 

0 

line 

peak 

40 

Table 

.11.3. 


Absolute population densities (cm ) at t =0, peak, 
and t = 40 nsec corresponding to data presented in 
Figxire 4.7 


Absolute population densities (cm ) 


3.9+09 6.4+09 2.8+10 1.1+10 

1.0+12(13.4) 2.5+11(25.7) 8.8+10(26.7) 2.0+10(38,6) 
1.5+11 1.9+11 7.7+10 2.0+10 


3.2+09 

2.6+-ll( 

7.8+10 


9.5-409 
7.7+11 1 
2.4+11 


6.0-409 

7.8+09' 

7.5-409 


2.0-409 
1 . 1 + 11 < 
4,3+10 


Gas 

pressure 


Curve 


t = 40 nsec corresponding to data presented in 
Figure 4.8 


Absolute population density (cm ) 


3^D 



0.4 torr 


5 torr 


dashed 

line 

0 

peak 

40 

2.7+09 

1.2+13 

6.4+12 

(6.5) 

1.1-409 

2.5+11 

2.4+11 

(29.3) 

dotted 

line 

0 

peak 

40 

3.0+08 

3.8+12 

1.5+11 

(3.7) 

8.9+08 
4.3+12 
4 .0+11 

(6.4) 

dashed 

line 

0 

peak 

40 

7.8+09 

7.8+12 

7.9+il 

(5.6) 

1.4+10 

1.5+12 

7.1+11 

(16.0) 

dash- 

dotted 

line 

0 

peak 

40 

1.4+10 

8.9+12 

1,6+12 

(5.6) 

1.7+10 

3.2+12 

1.6+12 

(15.0) 



159 


Table II. 4. Absolute population densities at t = 0^ peak, and 

t = 75 nsec corresponding to data laresented in 
Figure 4,9 


Curve 1 

Time 

CO 

0 

(I) 

—3 

population density (cm ) 


(nsec) 

I Level 3^S 

1 

t 3 

• 3 

; 3D 


0 

5.1+02 

1.2+03 

2.4+02 

dashed 

line 

peak 

4.1405 (51.5) 

8.2+06 (7.0) 

1.3+04 (43.0) 


75 

3.84+05 

3.5+06 

1.1+04 


0 

1.7+05 

1.3+05 

1.1+05 

dash- 

dotted 

peak 

Q.QaOI (51.5) 

1.7+09 (7.1) 

2.9+06 (43.1) 

line 






75 

7.8+07 

6.9’+08 

2.4+06 
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APPENDIX III 

SNOBOL PRCX3RAMME 


This programme, written in SNOBOL, prints the inverse 
of a square matrix in algebraic mode. In the programme, N is 
the number of rows of the matrix. The dimension of arrays A 
and C should be made (N, 2N) and (N, 2N, 8) respectively. The 
output obtained contains the elements of inverse matrix in 
algebraic form in various levels of indexing. The programme 
itself taxes care of index referencing and defines the 


variables \ased. 



0 u r i 

01 0 4U 

0 G f: ■' 
0 'J i'» b 
nijn /!'■ 

Oiff' iid 
OuOyr! 
O'j 1 ']<■ 


Y = 'A' 

= 

7 = ■> 

M = *4 

Mr'!*? 



beta 

nEf.TA 

AA 




0490 

0500 

a540 

f§dO AJ 


I tt 0 
I » T + 1 
LeUf H) ;F(AA) 
MECi, K) JPUJ) 

•J = ■« - I 
1 3 J' + 1 
rM(J, ^J iP(AG) 

TOPNfCAh K *0n 
lUfNTCAtT.JJ . *0n 
ACirJ} » 'r A 

ACTiKJ « *0' S(a 5) 
FULtj « 'TRUE' S(PR) 

guTe*wr •* ■ - 

output "s ■' c -ris- 

OUTPUT JS *• c ■-Ilf 
|UT-PUT « ' ^ 

I « 5,+,, 

t'- ■*' 


'WMM 




*CAHJ 


00700 
00710 
Oijl'iO 
OU7JO 
00740 
0075)0 
007t,f» 
0U7 70 
0U7W0 
00790 
OUROO 

OMlO 

n'(r?4io 

00930 

OufiOO 

ooasi) 

OOgoO 

00870 

00880 

00890 


; '' ’ \ - 


mW-- 

'VfiA !-h'- ‘I- 


■ ' ' 






1 i '' 

Pi 

k Uol3i'' 
’ 0 U 1 4 /' 

001 5f> 
t)U 1 bO 
I 0U1 70 
' OUIB-' 

001 9') 
0U?'>0 

I ou?n 
' OU77IO 
00230 
0U240 
> 0U26U 
' oy?o.) 
00? H) 
ou?ao 
OU290 
00100 
003i0 
OiiUo 
00130 
1* 


THKTA 

ZETA 


RU 

CAMMA 


alpha 


J ')' icznisj 


A = ak^’AYI. '4,8*7 
: ~ AhPAYI '4,8,8') 

T 3 T + I 
Led. U; ;F(RO) 

,1 = 5 

T = 1 4- 1 

r.e(J, M) JFITHETA) 

HdU = 'A^ '(* I 

T 3 t + t 

LECl, U) :F£neLTA) 

.1 3 S •■ 

J 3 J + t 

Le(J, M) :FtGAMMA) 

K = J - PI 

AitfJJ * ,1? sCAIiPHAI 
K S 0 , 

^ 9 |C’ I, 

,;P<S|2|(A£1,(N + 1)J), oO) .:SCAZ) 

.¥., * .. 1 

• 0 

■>!»• P ♦ 1 

IfCPjt M) JS(AZ) 


>» I ■ 0 + 1 

'itiOt :S{CE) 

=.Arp,Q 3 

‘'“%4t?I:ull,^^trur' ? :^(Kh ^ 


fCFj 


pfK, 
a K - 1 

,T s ,T + 1 
t.E(J, >1) JF(AD) 
NE(J, K) ;F(Ab) 


TaFNT(AlK.JJ , ' 0 ') :S(AE) 

nhv> i -‘o- ''■ 


OU^u'.' 

00‘>1 

•10 

Ou* i > 
DO'l'j ^ 
00 Ob"' 
OUl / o 

0 . -1 y A 

OiOu "* 


niiWyj 

^ L 
PD 
Rrd) 


r,tfi, jj ;F(Ru) 

I = -.t 


T i I'.::R = 

T = T + 1 


r.KfJ, 'j 

TjR; 

T jT’r.R = 

= AtT,.TJ 
lu^VlT = * 

nurb>Oi' = 


Ft AO) 

BriA>lrO 


T P’RR 


:S(DOwN) 
AC1,J3 ' 
CJ - 
'X* I '=' 


*' *B' tJ - 
NJ 

THTRR t(Ail) 
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APPENDIX IV 

PROGRAMME TO CALCULATE P0PU1ATIC»I DENSITIES 


This programme, written in FORTRAN, calculates the 
population densities of Hel and Hell levels for both optically 
thin and partially optically thick conditiais. The input 
parameters (besides those defined in the programme ) are as 
follows: 

WL Wavelengths of Hel (l^S-n^P, n = 2-6) and Hell 

2 2 

(1 S-n P, n = 2-6) transitions. 

OS Absorption oscillator strengths for the above 

transiticaas 

EAE, EBE, Discrete electron energies (ev) at which the cross 
section data for various transiticaas are given 

EEE,EFE 

CSl Cross section data at electron energies EAE for 

transitions 1^S-4^P, 1^S-5^P, 1^S-5^S, 1^S-4^S 

CS2 Cross secticaa data at electron energies EBE for 

transitions 2^S-3^S, 2^S-3^P, 2^S-3^D, 2^S-4^S 


CS3 Cross section data at electron eneirgies ECE for 

transitions 2^S-3^P, -3^D, -4^S , -4^P, -4^D, -4^F , -5^S , 

_3 ^3 ^3 ^3 i 

-5 P,-5 D,-5 F,-2 P 
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AF 


PA 


EA 

AP 

H 

BP 


EB 

PB 

E 

a 

AAP 


AAfi 

CP 


DP 

DE 

EF 


EE 

PE 

PE 



Constant a (Iki. 2.17) for transiticais 1 S-n^F 
(n = 2-5)/ -n (n = 6-12) 

Absorption oscillator strengths for the above tran- 
siticMis 

Ejccitaticsn energies (eV) of the above transitions 
Roots of Laguerre polynomial 

Wei.ht factors of Laguerre polyncxnial 

Constant a for transitions 2^S— n^P (n = 3-5 ) , 

1 PaI 

2 S-n (n = 6-12) 


Excitation energies (eV) of the above transitions 
Oscillator strengths for the above transitions 
Energies (cm”^) of Hel levels 
Weight factors of Hel levels 

Oscillator strengths for transitions 4^S-5^P/-6^Pj 

11 1 1 11 11 3 13 

4 P-5 S,-6 S^-6 D; 4 D-5 F^-5 P,-6 F; 4 F-5 D,-6 D; 

3 3 3 3 3 3 3 3 3 3 3 

4'^S-5'^P,-6'^P; 4 P-5-^S,-5'^D^-6'^D; 4'^D-5"’p,-5'^F,-6"’p; 

3 3 3 

4 F-5 D,-6 D 

Excitation energies (eV) for the above ^transit ions 
Oscillator strengths for transitions 2^P-n^S (n = 
3-5), 2^S-n (n = 6-12) 

Excitation energies of the above transitions 


Absorption oscillator strength for transitions 


2^P-n^D (n 


= 3-5), 2 P-n (n = 6-12) 


Excitation energies (eV) for the above transitions 


Oscillator 


strengths for transitions 3^S-n^P(n = 3-5), 


3 S^n (n = 6-12) 




Excitation 


Oscillator 
3^p^ (n = 


eneirgies fev) i'fOgg-.'.'the .aboVe -trehSitions ■ 


r 3 P-n S (n = 4 #5)^ 

'.fi< -li ' 


Excitation 

Cteoillat'or, 

1 

3'‘.D^v.ijn..# 




-,1^ ! 


'Kiilil 


tte '-abo^ /transitions / 
t^itiobs P (n =^3?-5),, 


f hr "the above tr^siticns 
.^B^^li#i"’il^-%ransiticns -5^F,-n 


(n = 
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HE Excitaticai energies (eV) for the above transitions 

1 11 

OF Oscillator strengths for transitions 3 P-4 D^-5 D,-n 

(n = 6-12) 

OE Excitation energies (eV) for the above transitions 

3 3 

PF Oscillator strengths for transitions 2 S-n P (n = 

2- 5), 2^S-n (n = 6-12) 

PE Excitaticxi energies (eV) for the above transitions 

3 3 

QF Oscillator strengths for transitions 2 P-n D <n = 3-5), 

-n (n = 6-12) 

QE Excitaticn energies (eV) fjsfir the above transiticais 

3 3 

RF Oscillator strengths for transitions 2 P-n S (n = 

3- 5), 2^P-n (n = 6-12) 

RE Excitation energies for the above transiticais 

3 3 

SF Oscillator strengths for transitions 3 S-n P (n =3-5), 

3^S-n (n = 6-12) 

SE Excitation energies for the above transitions 

3 3 

UF Oscillator strengths for the transitions 3 P-4 S, 

-5^S,-n (n = 6-12) 

UE Excitaticai energies for the above transiticais 

3 3 3 

VF Oscillator strengths for transiticxis 3 P-4 D,-5 D,-n 

(n =6-12) 

VE Excitation energies for the above transitions 

3' 3 ' ' ' ' ' 3 '' 

WF Oscillator strengths for transitions 3 D-4 P,-5 P,-n 

(n = 6-12) 


Excitation energies for the above transitions 

3 3 

Oscillator strengths for the transiticais 3 D-4 F, 
-5^P,-n (n =6-12) 

Excitation energies for the above transiticais 
Values of 4Q _ in Eq. (2,19) for transitions 1^S-2^S, 

*1 "I ■ 1' ' 1 

-3 S,-3 D,-4 S,-4 D,-5*S,-5 D,-6 D 

Excitaticai energies of above transitC^^'‘;,,.v.t - '-m 

Excitaticai energies for transitions":, I, 

Values of ^ ^ 


3 1 

3 P-3*P: 
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CRQ Val-ues of 4Q for the above transitions 

SFE Electron energies 

SFF Cross section data at electron energies SFE for 

transiticais 1^S-4^F,-5^F 
ES Electron energies 

QQQ Cross section data at electron energies ES for 

transitions 2^S-3^S/-4^Sy^-5^S/— 2^P#— 3 P/-4 P, 

-3Hz -4^Dz -5^Dz -4^F z -5^F 

APQ Spontaneous transition probabilities for Hel 

transitions 

ETT Electron energies 

QQ Cross section data at electron energies ETT for 

3 

transitions from level 2 S 



-■ no ooo nom rso o 'joooooooooorior:»'ioo'3000o- 


I .!-.A:--i c-'ilrulrites t. he oonulatiori if'a': i t- les ( 

!| ! lovf'ic f^vr hotn onticaUv fiin -ini 
HIV’ ronaillons. Tmoortaat ter-ris 

; tell nv.s ; 4 a >V 

R t.'=>ci‘ i h(! uemoereture in ev ID i 

;i eci. run densi tv 

: dnpuiatijn liensltles ot Hel and neiT levels 
jprt'/ of levels in inverse centlTieters 
M nt f ect drs 

ie I iro'Jivi state nof’ulatlon density 
. a s L e ni o e r . 1 1 u r e 1 n e V 
lotlral esrane rectors 
;=n^(lj/nf+ 

Electron impect excitation and deexcitation rate coefj.*;, 

: »eutral-nputrai excitation rate coeffs, 

Ki,rp s Central neutral Ionization and recom, rate n,)rt^s 
noontaneoas transition probabilities 
: bailatlve recombination rate coefts 
Tscti lator strenoths 


V^r 

Pc:^P 


laolos with IntPqer naniPS are 
•<::,K!iPc,K'''iCP,K(^PO,KNsn(M!,wC 


oeclar 


real 


Kb^,KP0SJ'1,''P, iO,KGDv'»<<t<PO,Ti>,KP,f':pX 


.. h imens 1 uns declared tor subscripted variaoies 
Pi'^fclMST JM TtFrC 1 ) ,uRFf 10 J , wT,t 1 0) ,nsf lO J 

i) ,H:Ar 11 J , APC1 11 ,-U 1 il ,HFf n J l ) J 1 j 

4r4t‘ i Vh ^ TU ^ J 1 2 ^ f GF ( 1 0 ) , CE ( 1 n j , H'-’ 1 1 u ) , lit f .1 0 j , jp ( 1 u ) , m-: f i ) ) , 
if Cf /} fat*!?) , AAf (20) f AAE(2o) ,Ct f 10 ) , CP ( 1 0 1 , dp f 1 0 1 , DP ( 1 0 ) , 
iPPf i 1 j ,PR(U ) rOtliO) ,gPClO) ,»Pf 10) ,kP( 1 u) ,.*?KiO) , SPC 1 Ji ,UKC in i , 


r I i w f r ^ I / ^ vr i M / , Vl'. i u j , «r I JL ^ , Kf. V ' u J ». 1 O » 01* l ! V I , ' r (, I I » 

4 tJFtl 0) # vrc 10 J . VEC 10) ,'dFf lO) , Ul) ,Xi* ( 1 0 j , A'V( 1 ^ ^ 1 u) , OPf-’c t C i 

STEEt 1 n ,0rP( t n ,CKEt t 2) ,CRQl 1 2l ,StT(1 2) .PPf 1 0 J ,00of 1 ?, 141 ,CK( 1 2) , 


6Pr(32i,K!<PQCl?),00Ct7,l2) ,PTT(12) 

DTMP^SiOd KPtl2) f SKF(1 2,14) 

niMEdSlUM PFP( I2,i2) ,KOPC(:)2) ,XDrp( j9j ,KMPOc =<2,32) j2) 

DTdE^Sirju HEIfl^i ,Hai(15) .ilAi^OHh, I*' ) , d<PO ( 1 b , l 5 ) , -iXP^ (1 o 1 
1 ,HMPC15) ,d^0(l5} ,HKPSUd(1 'j) ,M0V>Slid( l*^) ,UKAT (47 ,47 J ,uPXSuM( 47) 
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< i’ I 1 b , I i 

• 1 f M f- A , I J 

4 , I (C.5 I. f I 

b'-«4 , ( rc5^; f .i 
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-"u 6ft 7 J=1,5 
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■’J 6ft6 
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n=:5i(i,j’)4-i,nfc;-u2 
O.j 6ft9 .1 = 1,21 
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7 J ->79 I = l ,2 4 
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7J 67 1 1=1,2 

o;j 671 1 = 1,1 J 

T , T j=C,54(l , J) ♦I ,0fe;-02 
OJ 57^ 1=1,13 
7;i5(T,1 )=:55f 
f)J 22 T = 1 ,M 
f' J 2 2 .1 = 1,1 

APOCTf T)~kPjf I ».1)^1 , Jt+w-ft . 

IHj 6fti3 r = 1 , 32 
nAP2a,i)=4pa( i.u 
Frj«,^iAj'r iK,i jFi3.5) 
5eft.)2*,(R:i'T(T),l=l,l2} 
FJR=4Ai'U2F5.t J 
DO 26 1=1,17 
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DTS(r)=^2Tt£(l)^1 ,Dt;-j4 
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000(6, J) = 000C6,J)^‘1 .Oc-0? 
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T=9 

TU51C11 JsCSUMPf S'] -'-0 + 7(0)3 41,0/3,0 
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SU'U = sn''lX+2,-) + )HJJ 
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1=18 
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K;X’=S.'"ti;(i5)-KCK(l5) 

Ton + rn j = U)t. ibUOU 5) 

roM T n=Tj41 ( I t ) + 1'!) ,(21. T i 1 +.T0'' jC ( 1 ) + r'f .. 
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t , V/7.i •:> ^5*-r:t if J ) ) /I'f r^-p+fl .S j 

^i<c'i 7 .dtiou rate coetfs. arf> cdiriU=»tPo jv 
,- rro&K sertton Pxr'ressioa ur.ln.i uajuefr" 

xF r i" V ? It'? Stored In different arr= 5 vs . 

iilonTs?,ii 

^.i,i=£At 1 ) /' 3'7 

■ = 7 o? /3 + Af { 1 ) f Fa{ l)rP:xP(-Ar )/{. Ai,*< Te;'^‘P*n , 5 i V 

r'j "e.; T=1,1t- 
J)/Au 

Ab=l ./, 5 r(Z + l.U) 

6 F.I i = 7 .f Ai/.TdfAGj/tZ + l .nj 

AS:! -I^ATyrf !l( T J 

5 ll''<='>u '1 + A isd 4 . 

ClJ'.' CT.V tj 

AXf A 1 =Su 7 f A:; 1 „ 

DO 76 1 =l,li 
ALsfcifUT }/tp: 

RCO < = d^ 41 o 2 / 3 fRK( 

80 '' 7 b'’j 3 i , J I 
ZrAPfT)/<vr, 

AM= 1 . 38 »(’Z+i .01 

8 Sfr«l = ^*Ar,[jr,(Ad)/(Zfi.O) 

FUM = n 5 '.rin»clC J ) 

SUf^ = 5 UM + Pi)M 

C.OMTrMUK 
DKf J, ) 3 ;SUM» 920 («i 

8 ^ 3108677 . 0 
Fi=l 8 M 3 l 1 .0 
ACt)=:(RSf/Et ) 1 '* 0.5 
Du 80 Ts?,M 
A(T) = (R)f/(EA-b(i)))f*o .5 

cDunMue 

00 80 TsU ,25 
DO yB J= 27 ,M 
ftNaiid ) 

Ajsi.U/C AN+AM) 

A 23 l. 0 /t AM*A aJ 

Aisl.U/( C(Ai-A 7 jtAd+A .,1 t*i) 

FP 0 (T, JJ = l.y 6 uUAJ/t Ar'♦A^‘) 

r’OMri.jtJi*: ■ 

no J 05 i=i ?,25 

no 185 J=k 7 ,M 

ALsCKC TJ-e(X) j/i’rE* 8064 ,r)) 

A|jr:fj 835 , 4 if>^ 73 *FqQf t r J 1 ♦ FXP t -AT. 1 / ( Ab* ( T'FMpr’t 1 , 5 1 > 

TAM. rTMrtAP,H,AL, su'd 
KPOtT. J)=AbCOU*Sin'l 

no ilo 1 = 26 , f'l 

, 1)0 120 .1 = 1 , M 
IF tl.SO.Jl GO TO 120 
AN»Ad ^ 

AM«ACJ^ 

ASsi.O/CAMtANJ 
A 6 s 1 , 0 /(Am*AM) 

A 7 *t-. 0 /{ f I A 5 -A 6 j 
FPDCIx .d = l ,960 3 *A 7 /CAf«*'AN) 

.ft 6 »iriJl::y(I))/C 0 O 64 .O*TE) 

CAbb FT«TCAP,H,AL,SUH) ■ ■ ■' 

-■KPOCTf 7 ^»ALC,^snM 
cgMiiwuti: 
coMtiMiie 

no 213 lsl ,20 ' . - 

AbRAACCU ' ■ 

AbsAt/fi; 

AAC 0 ^» 5 . 4 | 6273 '*iAa 5 (ig*€XPC*AL.)/(Ali+C!rEMP^* 1.5 1 ) 

"SSiltc 52l5SSiiS4§ii'!f ? 
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923 
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;S T 
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U 
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r*» .0%-, 

W . 1 

f, ==; 


in 4 1 A * C 


'■'xr,. 

, 1 

T 

1 A P , 

i'i f 

9 b 

n iv n 

i ^=2 


i ' 2 1 

c;k 


f'U < 

, .7 r 

=: 1 

, i " 



f 1,1 •=, 

) i T 

) / 

1 



t/ 


m 

m ^ 

5*727 



C‘-\ra 

. r’T 

I 

1 1 r , 

a , 


r^v 

.r 1 1 

1 

J'" i-n 

C ’1 


n/j * 

n ^ 

q; ! 

, 1 U 



f i.» * J* ’ 


) / 

f f-' 



1 , 

, = . 

1 152 / 14FF 

91 

CAfL, fl 

41 " 

t A t> , 

•'I , M 

£KiT )~ 

t f*; a 

J -i 


1 


1 

= 1,9 



Au = ,- 

''•'t r 

j / 

IT 



1 


1 11 


1 n 


121 


133 


143 


151 


163 


173 


183 


193 


203 


303 






>.’' 1,1 / rAf,*f j .b) i 


i .3 1 j 


» / 


r-’,iP) 


.“i J ) 


FC 1 . i3o2 / 3*sJXr>(-Ht,) *KF{ T j / ( ( TF 'P** 1 .«S n 

'"Af.i, ^ TvT(A?,ii,A.l,,F-,jrj 

F';\f n==',j!rtt.-;- 'If, 


s • 


1 i 1 r-1 , 1 


)i 


) ) 


s n 


ALsuFU ) / '.IF 

r,C"l . 1 3b27 3 + flK( 1 1 <^ ( - i\l. ) / t Al* I 'T'FP'P* ♦ 1 

TArn, f r.prcAP.H, a 1 ,,FuF) ^ ^ u I .. 

f:k( iisr.uF»Gan« 
no i 1 .3 r = i ,q 
Ab=HFtr)/rF 

HCq.'< = 6 , 1 362 /llMf- ( n *FaIH-a1j ) / ( A (,♦ c 1 

FAr.u f T,'jTtAP,H,AfMHiiF) 

fin 123 lsi,9 
At.=ne(i5/T£ 

nciri = b,4.3b273»Ob( 1 1 ♦FXP t - hT ) / 1 ^ 1,* C H-'* * 1 
CALL FI^U’(AP,H, AL,unn 

no uj 1 = 1 ai 

AL=PF(TJ/JF 

Pcnt^ = 5 . 43527 31'Pf r n^FxPl-Al.J /(. Ah* 1 ,S j i 

call FiMT(4P,H,AL,PnFl 
PK(T JaPOF^PCON 

no 143 r=i,io 
ALsOKlU/rF 

OCriLsb ,43b2 /340F( JL ) TPxnOul.) / (. Ai-> t vk a' J ) 
CALL FrimAP,H,AL,OUFj 
OKCi)=OOF40COH 
DO 153 T = ia5 
ALsRFtll/TF 

RCnA) = 5, 43627 3*RF( 1 1 +eAP I -ALI / C AL* (Tfc'‘.*P*» t . 5 } ) 
CALL FfMTCAPrH, AF./POF) 

RfCCDsROFl'RCO'l 
DO 163 i*iao 
AL=SF;^T)/fE 

SCOM = b* 4362 /3 4St' ( 1 ) ♦FXP t - A?.- ) / ( A 5 * (Tirvp^ 1 1 s> 1 
CALL FI,a( AP,ii,AL,50r) 

SXf llsSUFtSCON 

no i7i 1=1,9 

AL=UPtTJ/rK 

f'C lb = b.43b2ni‘Ut- C,l)»FA'M-fJ ) / ( AL* (Tr/''P'1^ + 1 ,5)) 
CAf.L Fli'jri AP.H, ALfOUFj 
flKf U=70F'i'UC0Ai 
no i 83 ‘ 1 = 1,2 

AL = VEUJ/TF 

V : 0 w = 5 , 4 3 6 2 7 3 4 V F ( JL ) t F A P t - A T . ) / ( A L * ( T K M P * * 1 . 5 ) 1 
f’Af.L FT^TC AP,H,mL, VUFJ 
V((f n=i^jF+7:oft, 

DO, 193 1 = 1,9 

ALswF(I)/IK 

WCOAi = b,43o27 3*WF(I)4EXP(-AL)/CALl'(Tfc.MP**l .5) ) 
CALL Ft'HT(AP,H,AL,W’aF) 

DO , 203 1 = 1 1,9 

XCoS«5,436273fXF'(l)*EXPC"AL)/(,^L'»l:tT'EMPt*l .5)) 
CALL FfNTCAP,H,AT,,XOF) 

Xt(( naX0F4XC0N ■ 

DO 3.93 T = 1.8 
'ALaSF£tI)/tE' 
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SSCR=7,3S8721|-3*SEECl')4‘SEE(I)*OEECl)*FXPt-Al.)/(AL*(TFMp*4l ,L j) 
CALL,3Tr'jTCAP,H,AL#SUH3 

SKKcnssiiMi'ssc^ ■■■■■■- -' 

00 Moaflrll , ■ ' ' 

AL»fFB C 1 } /ti . 

T£Ca^»7.4587ai'S'A3''fTEg(I3*TEECI)*0tE(T)'''EXp{-AL3/CAL4CfFF)Pl‘#l.S) ) 

' 

..■■ ^ I M-■■■liili^i^ll■^■■llll 





'* \;4 


4;) 7 


1 

>' I.J ( 1 , 1 1 = 

^ i , i j 
.• *" I. , /i'M 
'■'J i.'*/ r- 

T I = I f ^ J 

K i , t T I 

'<'P''H ^ , 5 j = r^ 

’ I 1 = 

K-p'^l 1 , i :; ) = 
KH7i ■-* = 

'''J ’ I 

I i = I + I 


^ ^ f / f. } 

’ r: 'I ?.'/ ) 

™ tTK f j I 

= !r , f / 1 

S , J 1 

=:.V'< U' ) 

C f ) 

'<CK(^S P 

<<('+) 

■" ) f i. I 


1 

1 1 

= Hfl ) 

7 ia 

:: K 1 ! ) 

1 

i 1 

=D.n2) 

7 

0 5 

='’<f n 

1 

n 


1 

■n 

= 7< r2 ) 

n 

3^ 

=D<f i) 

7 


1 = i , r> 


OJ 


’•'PIHH, fAls" 

Kk>-)il , 1 I lar^i^Q ) 

Kpii 7,?o)=F:^f,n 







4 0*^ 


41S 


417 


419 


421 


nj 1 J 7 Xa 4 ,i.,» 

U = I t- 2 > 

Ki?7(7,U’‘='^N( 1) 

K?7tU ,nj5:fKit J 
KPTtU , 2 r)A=:f<C?l 
’’■P 7 (i 0,1 nsr,s<(n 

K£i :iU ), l 9 )s:.jK '(2 j 

^P 7 ( 1 J, 2 '^ 

fPO I I 0 , 1 7 ) =HK ( 1 ) 

KP ">(1 U, 24 )a,iK(?) 

no ns i=i,y 

T 1 = I 4-21 

KPatlO,lT)a,lK(T) 

KPOCll ,l 6 )=r)Ka) 
KPOCit , 23 )=a(v( 2 ) 
na 417 £=i ,9 

11 = 1*23 

t^PyCil ,t.n=nK(i) 
KPy£ 2 , 4 )=HCKC 20 ) 
KPC)( 2 ,§)=rtC.<C 9 ) 

KP0C7# H) = HCK(121 
KP 0 £ 2 r 2 l)=HCKC 16 ) 
no 419 T=S,ll 
ri=i* 2 i 

fCP 0 ( 2 ,tl)=O(C(l) 
xpocnpjsQKcu 
KP 0 t 4 , 15 )=*K( 2 ) 
KP,;C-i. 22 l = yK (3 j 
OLl 421 T=niO 
Tl = I »-22 

XPOC'l,!! 1 = 3f£(l)tRK(I) 
(<Pjf 4.o)='Uf i ) 

Kpac ni 2 ) = 7 *< 1 ( 2 ) 
KP 0 i 4 , 1 9 >=RKf 3 ) 
!fP 0 { 7,9 J=SKC t ) 
(vP\;(o,l 4 j = 3 KC?i 
Kp7(3,n) = S.<( 3) 
KpHS,2ol=.SrC(4) 

KP 9 tP f 1 21 =inci) 

Kp^C^f 1 9 ) = 'nC 2 ) 

KPOCcJ, iSls^KU} 

KPUf df 22 ) = \?K (21 
DO 429 Ts 3,9 
Tlal *23 

KPaCBnTJ-VKdJ 
COMITilUe ■ . 

KPDtPpI 4 )e' 2 i£f " ' 
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streaatns tot ipi transitions, re-iu -nj star 
are restored in arra/ of ^Inonslon 


f- M )=0,.)<,74'> 

15,l(Sj=5.n9'Sfc:- 

ppPt I ,57=Fnti ) 

,insFAC2) 

F-f^Ci ,iH)=FAf J) 

PF;M 1 ,251=Fa( 4) 
nj 05ai 1=5,1 i 
Ti = U21 

Fef>t1 ,Tn=FAf 1) 
pf'Ptd, rn=Fdf u 
PFP(2,Il)=nFf l) 

Ue 

FKO(d,5)rFHt t ) 

PFPO, 1 n=FBf 2 ) 

FKPtd,1d)=Pt5(3) 

Pr°C 1,25)='^b(4) 
FFP(1 J,25)=AAl*(l) 
PFP( 1 3,26j=AAFf 2) 
«'FPti a,20j = AAFf.i) 
PFPaB,26)=AAFf 51 
PFPUor24)=AftF(b) 
rFPUb»25)=:AAF(?) 
PFP(16,26 JsAAt («) 
PFP<.17,2 3) = AAF C9l 
FFPa2,2t )=5AAFru ) 
PfP(.l2»26t=AAF( IP) 
FFPII 1, IP)=AAF( t 3) 
PFPd 4,2? ) = AAP( 14) 
pfP(14,26)=:aAF( 15) 
PFPt15,2). )=aAK 16) 
fF!^tl5,21)=AAF(17) 
PFPtt5,26)=AAF( IR) 
FFPtt7,2?)=AAF(lP) 
PFOt I 7 ,26) = AAFC2U) 

Pi- Pt5,'^ ) 

PKP15,1 i)=CF{2) 
PF^t5,7y IsCFCII 
DU 65 B? T = '1,10 
T t=if2? 

PFPt4,ni=?F(l) 
FFP(.5,TJ 1=l)F(i) 
Pi-Pl6,n)sSF(l) 

pP“t7,Ti)=FK( 1) 

."OdriMUF 

PFP(5,1u1=nFf 1) 
PFP(5, 1 6)=nF{2) 
PFP(5,?i)=DF(3) 
FFPt7 ,1 i 1=PFCi ) 
PFPt7,18)=EFC2) 
PFPt7,25l»EF(3) 
PFPO lajlsFFU) 
PfPtl 1,20)=FF(2) 

DU 6583 1*3,5 , 

Ifpu i'!n3fQrci| 
rrpnof iP'Siircp 


6581 


6f63 
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r 

r, 


c. 


467 

4o6 

C 


I 7 

K’f 1. 7 
t, 7 
1 1 
7 ( .1 
^ I 

“'r 1. “S' 

r ;.•!-> ■) 

"r t * 

7,. M'-. 


7 i=. 
1 i 1 

70 

'! i ■= 

> /- ) 

I i ^ 

1 ) 

q j = 


'■TO 
f .}) 
-'FT O 
'•’ ( 1 ) 
■OT:,i> 
■'.T J ) 
--'(1 ) 

( /•) 
( if 
-'Ll) 


Pf-' 

n 

dr 

1 

j ) 


(■2) 

rr./ 

n 

16, 

P 

i 1 

.iSP 

( i ) 


r> 

d, 

1 

/ 1 


( J 1 

li"* 

!*> 

L >' , 

I 

o 


(2) 


1.7 

d, 

1 

'> 1 


( 1 1 

iri- 

o 


y 

* ) 

= 'Jf( 2) 

Fr’ 

r> 

<• 

1 

41 

= ^'F’ 

f 1 1 


:> 

1^0 

7 

1 1 


(2.) 

iri/ 

P 

I'^r 

1 

1 1 

=xr 

( 1 ) 

rrpp ^ 

y 

n 

= :<F 

(.') 


'a7i7.iOne is called io caicuia'’.e tOiit- ra j -i^utrai 

dUoa, uni'^atinn „ni rac^ i,h lOrtf ion rat? co^tfs 
r J H A, OfcJTTG^TF.FF'P.KnPC.is -C ^ , *< .'i f-’U I 

,, ‘<jO;,s are sa:n-''ed 
aj tf'p T = l,i2 
RU’.'=:n.;) 

PJ UW j-iO? 

TFfl.iO.IJ GO n 467 

Kdsii'U r )=sdq 





40 3 


,, Pe'naidlaT electron irr^pact exr. rate coetfs. ar» sLor«n In 
arrav KpO 
Kpn( 0^i)=S6K (7 ) 

^POl 1 f?)-riC<(2^) 

KPJ(1,4)=HCF(21 ) 

KPOfI ,6j=rFK( J) 

KPOl j ,l2)=Tfe;K(6j 
KPOtt , 9)=T£KC9) 

KPOn 

KPOd ,9j=lC^i(53 
KPOCt ,12)=TEKCft) 

KPOd ,14)=TfiK(73 
KPOtl /15) = Tt*«CS) 

KPOCI ,19)=reK(9) 

kp6i,i #20=TeK(i 0) 

Kpacl ,22)=T£K(U1 

no 403 1 = 1,12 

AL=CREC n/T£ 

CRKCT J=7 . 35a72n>J*EXpf 
K'P0t2,3}=COKf n 
.KP,/f2,7)=CRK(,?) 

Kpni7,13)=CRKt3) 

KPD(6,? J=CR(<f 41 

»fP0(6,1 3)=r!KK(S) 

RPOf o,2(0=CnK(o) 

Xppd^.,! 3j = :Ri<(7) 
xpod 2,2P3=:RK(d) 

RPv>f 4,5)=Ct<K(9j 
Kpoc 4,1 I)SCRK(IO) 

KPS.T4 , lR3=CPKf 1 1 3 
Kp7lf1,1 1 IsCRtCl^l 
KpO( 1 , ?d1=KpO( f ,2b?4-SKK(8) 

('0R = 7.36«72d>03/(T6Mp^>!<1 .5) 

DO 4S Isl,?. 

SIIMsO.O 

on 46 1 = 1,11 

rui+j 

PX=sP”c-Ciri-sFfi(i) 

POa-SFElD/Xe 

TERH«SFECl)»PK+SFP(I,a)»EXP(PO) 

SfTMsSOivitrGRM 
SPFf(j)5Sii!4#criR 
KPoh d7)»SP!<n ) 


,5)) 



3 7 
3t) 


p .=./T ) r J n f J ) . 

H.n / ’’F, 

^ ^ 2 ~ ^ * I ,i V ^ t *’ ' ■ *•'■ ) * W: P t T , J J 

V;'’ l = \:* '!’ S/i'^k .^ ) 

r ' ' i / f "! ' ' j S. ;‘ .;> \J_ ( J I 

^’,5 j=f.3pn( 1 J 


b) 

>*’ ;i P f o ) 


P •> ( 7 , 1 i V 

!(n='S'\;f .3) 

Kp't ( '> , ' b '>='<o!^pf ^.) 
•^’bi?,-;>n=/S‘->v’(ioj 
i 2 J - n 4 j 

rPH?,’ t^='^:!^(la)fPCKtl9) 
:<■ P b I ? , ^ j = I r ^ ) 
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‘‘ r ■:> , UM = JJV ^ . J , 

‘•'P’H ? , 1 = 1 j ' ( I h) 


r. 

c, 


4b 5 


450 

C 

c. . . 

c 


BHH 

c 

^ • I 


73 S 


0 -» t ^ I 3 ) 

i' K H ^ , 7 J : 

^'i ?t?, / /1^ ?C/1.1 /) 


ii n 

iCKf b) 


179 


. '’i.lba e 
J > b ■) 


uHlbrian pbpuiaLion densities 

i - 1 » 


(nor'T’rtl i7ei j ire 




U ' 5 = 1 . I i 7 5 C 0 it - ’ b / ( T (•; •' i- ♦ ♦ 1, - 


') 


750 

240 

C 

C • » I 


470 
460 
C . 

c 


T p = ;,J / 3,-|n * . ,) 

I ^ ^ • I i ^5 r 'j 4 55 

Tp = (^i-t( I IJ/dOo'! ,'■> 

'U^C !.)=7il 1 *:n:,si*txpf iP)/4,t 
"Pf rl = 'ip( [ IM .On-t-fo ' 

... ‘tiect'rnn Impact leeACitation rate coeffs. are cal 
ns 101 oriel ole of ietaliea halanrlnp 

U 1 b H j = I , 

oj r=i , T 

Kpn( T, T J=kO jci , J) ♦‘ip( t V5'Qfu) 

* • electron i’^'c-'act tonlr.ation rate coetfs, 

K = P,1 / d')o4 , 0 

OJ MO X = l,;i 

rt'f I 00 TO 2i5 

x;=rr:i-=:(T) VMU64.0 

"Dni’Tiinfc: 

K = X/Tc; 

COMST = t0.yOot}/tKXPC-X)/ 1X4^1 TP!, 

SU ’ = 0.0 

nn 250 J=1,11 

7, = 4P(J)/X 

ni = (Z/fz+i.ou*4i .s 

n2s:( 1 , 0-( i . 0/ t 2 .04 (Z+ 1 .0 > J ) J4^i.i) /i. ij 

r)i=2. /+z**o.b 

n5=Ar.or,.(ni) 

DbsHCJl+Dl + tl .04-02*05) 

StIM=S'J:i + D6 
KPC(i)=SiJM4C0MST 

• t. For 5 oartlc'Jiar level KOC's are sn'umed 

, 00 460 l = i ,fl 

SUMsO.O 
00 470 Jal,N 
IFtI.EQ.J) GO in 470 
SUWsSUvi + KPQCT , J) 

2aNiT,xrte 

KPOSftvK I)=.StN 


MUacp'i i)V 




5051 


Paiiatlve recomal nation rate coeffs, 
Cl = l5.0*6.7o02f':T04 

F(1) = iP8n!.0 
C2 = 3.')E4-1u 
00 “^iuSl 

XM=1 ,0/f (MT)443)*C24etI)) 
rtRl HT )=Cl ♦X Al/r TFMp4*U.5) 
nj ‘^052 t = 
voo=c2tt':(n 
7KK=4.l75Ftl5 
V01=V00+XKK/CAf 1)**3) 


for Hef and^ He IT calciuatea 





Vll J = X0l*( ( VKil/700)**0,5)/VOn 
P = “*0.o6b6&Gb6 

5052 GntT)=l,0-0.2l«*(XlT)**P) 

. no 505i Tst.^l 

5 050 RETA(U=BeTACl)*GG(X 3 *l , 0 E +10 
DO 251 Isl,15 
lI=T+32 

pas J ■ ' -1 . 

XCl)s157aOQ,0*4.0/tTEMP*tPf4<2)) 

5UM«0,0 ‘ 

00 j 5 a^j«i|U 

soMis.saMi'M ■ ' 



r 

C, .. 


tpctors tor rosonanre transitiw'is of r. 
‘*’1 = I'^n LPinOerat iir e , •■■!, = ..av? Th-ij*- 
■IN ^ jsol LI ator - ■ 

•r|=,r D 


f a.sr 


s tr en.*!^' a f Xi. - ni t* 


Hi 
3 = 1 f I 


^ fc I. = j. ,.,5 
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,A 


t.> 


T,NC,.tfl ,1 , Xr! ,us j ,UP(-') 


«24 


«25 

Pi 3 

C 

r 




c 

C » • 


490 

C • • < 


164 
C 
C 


r= 

X ■>, 

1 ) 

. 1.1 

TAl.u 

3?b 1=1,5 
1 = T t 5 

PrT 

P = P4 1 ,!> • 

T = T j, ■ , 

V ,', -w. ■ 

VLls^JLfj) 

osi-isf.n 

Zhhu OKFAC c P ,1 , -VC , '.’1. 1 , Xr , US 1 , 1 'tPF 1 

iFFiT)=jPr 

’!u'^rTfL'’F 

PRi-'i «}i,cnL';F( 1 j,i = i ,10) 

>='OR''inT(lX,1uFn.4/) 

For DdrLl;ili 7 ooticaliv tnlcF case, transition or or* ioi l l 1 es ari^ 
nylCtpllei oy optical, oscare tectors. For ooticaljv tnin case 
tnpse s taceiTsenf. s stioiUc* te r*‘iTioved 
?\P0(5,1. )=AAPQt5,t l^UFFd) 

APOtl 1 , i )=AAP^f 11 , t )*uFf (2) 
ftP0l.la,l) = !\APv(lR,1 )tuFKf j) 

AP0L?5, n = AAPw(2S, 1 )*UKM 

Ap0t?U, I ]=AAPv( 26,1 )*UFMt. V 

, Transitioa pronabll ities are surisr-.e^ ■ 

00 l = l,r« 

SUM=0.0 
00 49o J=l ,1 
SUM^SU^ItAPOC T, J) 
coNni'j'JE 

ftPOSIdCDaSUM 

iSooroutine Heri is called to calculate energies# wlegot fectors, 
electron Impact exci tat ion , 1 oni zat Ion rate coeffs., and 
transition probaoiitles for ”pTi 
OEF 1=OFF(6) 
oeF2=LiFF<:? ) 

' OEFisUFFCS) : 

nEP4=OGFC9) 

0£F5=aFF(lO) 

CAft. HFlT(Tt;,HFI ,HGl,HAP0,WKPy,HKPC,Hl4P,Hfj0,HKPS(J?t,HAPSJ'-1,nt;Fl , 
1uEF2,OEF3,DEK4,nEf55 
00 16.4 1 = 1,15 
T 1=1 

FcinsHtriF) 

GCTl)=lGTtX j 
KPGlT.t)=HKPC(l) 

Mpf I T J=HiIP(l) 

MOn. T)=HN0(.l) 

KPOS'JMCIT JrlJKPSUMlT) 

ApOSiJin(iT)=HAPSUM(T) 

00 164 J=1,15 
TJ = J + i 1 

KPDtTI, JJl=HKpOa, J) 

APOCTX, JJ)=HAP0(T, J) 

CUtllTNUE 






Direct ion.l»at ion'-excitatlon rate coeffs. 
f|i. l= 0 .u 
00 166 I = r,M 

Sll'^« 0,0 

DO 468- 0 = 33, MHE ' . “ • 


calculated 







tS*i‘ 



181 


0 

601 

620 

621 

6 U 0 
6001 


■f 1 = 1 - i 

Tt' ( j .I-:''. ) j u' Tii O'' I 
Oi! oi j T=1 , 1 i 
«:;(i , J 1J=-•*'C♦^p,(: I , jjj 

r J'l i ] ii .0:. 

’'i' f i. . i ■) = 0C * I C t T J + KPgSi)ru T J 1 + ApOb" .. ( I 1 
T I " I t i 


IF 


<' i 'll] rn 621 


Ajf i ,.j n=-o:i'6r\;C i,j.T)-->P'.'fu.i,T}i'.'.o(. !j)/, 'pi'll 

C'l.jT I Mti*-' 

ROf I . I )=o;;;*^PL•{ 1) 

0..ri,^)=PKT6Cl)/'P( 1) 
rj'iiiK'tt 

itij o'Jji l = t,^2 

^n( i , i )='lCnKPC( I) ItiFXSu^'lT j4f.PvSu H T J 1+f P0^n,vi( .( 1 
‘010.1=!, 12 


I 1 =T -1 

’■p'f l,i.D.1 ) 


C[j ID 617 



6 J.Q 

617 


6 2^ 

6Jf. 


1 7 6 

r 


n 

r 

c 


7bn 

807 

7 00 


157 


i , 

1 1 


1 , 

u ) - 

= 0 ( ,T J *K Gu 


i ) 1 

' GtfKivf’Cf 1V+ 

1 

U 1 

'1 ■’") 

7' U t, 1 5 

* f 

i , 

J i. 

j y ( J H K .. Pu ( u 

13 

f i H 



T)/. 'MT 1 
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ri i I .■! 1; 

'".jr i. , j)=A.jf 
f'3" I 

Jf i , i )=.’' j 
T j = i + i 
1 I’’ ( J . O 7 . / 1 

I'O n"> i » = 

A ..f f i , J } = A J 

7;)-' 

0) |7o J: 

7i I I 7 o i s I , ;j 

J ( J , ,l Is- iP( T j «ii?-:4l( T ,.T) ''i'i 7(.l j 

.. ’^ro'i iiOtrtx (\j, ^ new T.atrlx AOt”-' ( ■^' 0 X 4 ^ j t s tor’^al 'Viico 'Jg^s 
nut ::jnt„iiin first ani tnitrv ttilr i ro^^s, ana tJrst ana r.n;i.ri.v 
tnU i rruljTins of. .43, cilinllarlv first and L^irtv talri ro'?*; u'^ 
'natci.x al are removed and ti.e renal aim matrices are stored 
in 0 ae V array OOKa' l 'iSXlS ) . ^i-^p and tnlrfy' tOiri coin, >'3 

ot nl are nronnot In thlro ann rourtn columns nfc •* f r i 
n 3 / 1 ) J [ = 1 , M 5 
P:J 7on J = l,iS 
7isl+i 

JJ=J+1 ■ 

a\fr,vf I , J) = AJ(I r,.TJ) 

no 8u7 

OMP.-jf I, .0=63(11,31 
o,M!r-vf i,4) = :’».a 
n:gt.:wC l,i) = -An(U,l) 

Cn.gT 

no }52 r=32,n8S 

no 152 t=u,nhs 
Ti=I+i 

JJ=j'f2 

MgR«( 1 f J ) = Arj( II , jj ) 

DO 750 1=12, OHS 

,ii: = i + «{ 


2 , 7 > 



OEtlO 

i.OE<fIO 

35*MPXC1I) 

CII))/HPX(10 


is arff calculated aoct storeti in 
to^ ttie copulation aensity of ; ; 


A./ , /r" 

i^|./cUofi+l#ioeGReEf 2 .0) > 


B;'jp;v«ci, n=Bjri I, u 
B.VR/g(l,2l=4JCn,2) 
RNR«(l,J)=-.\nCTl,l) 
Ri'lE3(X,4)=-A3(TJ , m 
75,3 COWflMUE 

DO 4uS J=32,7rtS 
DO 408 1 = 1, tgS 

JJ=J+2 
IToI+l 

408 4M<i:W(J,U=ADtJJ,lI) 


716 


7222 

71S 

C 

c.,,. 


815 


Subroutine FSOLVE Is called to solve tne svstem at linear 
simultaneous equations, Tne arrav fi'd contains tne final solution. 
Population Coeff. ru = (FS( 1, 1 1+FS( l ,2) ) # rl « FS(T,.3), r2 9 FStIf 
rne elements of first rov^ of FS corresDoori to. level h and tne 
elements of 32nd ro« correspond to level 34 fipvei 1 = ^ of I'elll. 

fjSt 

DD 71 5 J = l ,4 
DO 7t6 i = l ,i4HS 
l)=B.')FA'(i,J) 

7Atl,. FSJtvfRt ANErt,dH,F,NEP.v','gHs,L) 

DO 2227 I=l ,4 hS 

i,J)=‘^(l) ■ 

I, = 7 

RO'Jf I.g.OE 

, f^ooulatlnn coefts, prlnteo 

OKT.gr 8i5,C('S’SCi, J) ,T = 1,MHS) ,J = li3) 

PiJPvSA r( I <,1 y'^13.4/1 X, iOLl 3.4/lX,iaFl3,4/lX, lOEl 3,4/lX,5eiJ.4/) 
pour 815, (FStI, 41 , 1 = 32,NHS) 

, Frot! population coetfs. tne pooulation densities of HeT leve.ls 
are calculated and stored In array OPOP, OPO*’(li corresporids to 
tne oooalatlon density of. HeT level Itl, vPx are tne actual. 

Sana equlilorlu® population densities, 

no 5D5> 1 = 1,31 


4) 






Sun 

5 U 00 


,, Pon J 1 -jf 1 Oil (lf?iisities ?q, noouliit inri u?;jSi^i.os of 

.L?V«LS nrlnted 

p,-( r . r s / 

’ '.M f i u, '‘’:iPt>ij.'rnn i urr-.,.*;! n p-; np Uj-^T 7 * vi- !!r.Tl LK^-'.u'-.' / i 
(uPijPi T i , l = i 

i "'I T *:>)'■■ 

. P.jPUf<4Tl.rt„ -uF ■. S i »' l 1 

If- np, Ci.p/tc 1 ) , 1 = 1 , 4 o-» 

C'l ■■ 

^li'J I'T, -J'!;:, 

BI’IP 

'-■’I'jn 


, . stiPr out me FSQfjvF solves a svste'^ o,£ linear cuole'j s 1 i.oi *• a'^yf 
ponyr lOQS by Gauss el iinlnat ior* 'necho<1, linoaF n .-i r -3 mef e’'s are 
<' = Mdfrix of .'Umension ^'X?' 

= rolufno vector Of ill liieos ion i.* 

N = Ma'Ttoer of equatloiiS 

t = 1 for first expcutior, 2 tor tortaer execatiois 
F = Tutpnt oaraTieter contalnlnn tne final solntma 
SUBKHtl Pi fit FSni//l£(H,H,r’f 
B.Mr.fSpru,’^ A(45,4b),B(4'j) 

nil t ' ^ C i ?.|r 4 3 )’ 9 fj 

Ml-jmi 

no 4 r=i ,Mi 
X = D.O 
K = 0 

no 2 . 1 = 1 , ■'! 

At =RCJ,r) 

Tf f AHbf Al J ,[,F.X1 GU xn 7 
XsAnSXAD 
K. S J ' 

CDMl’lN'-TE 
M £ W a f i ) = X 
no 4 J = i,"» 

Y*ft(T , IJ , 

A(T,.l)=AfK, J) 

ACX, U=X ■ 

H = l + 1 
no 4 J=xl 

A(J,l)=At T,T)/A(I,i) 

DO 4 (vl = Il 

ACJ,^)=Ar . 

DO 5 T = 1 ,141 
NMasNERMCl) 

acNnisBCu 

B(x 5 ?y 

lisit-i 

00 s jisu.N . ■ . 

BCj)=BtlJ-A(J,i)FB(T) 

P(l 4 j=a 04 j 

00 A H = l,i 4.1 

1 = ii 1 •! i + I 
lt=t+1 

UG 7 , J=ri ,il 

tif i)=B(l)-AC I, Jll'Fl.J) 

Fa)=ii(i)/AtT,n 
R E r d R U 

EffO 

'lunroutines PINT, SlNX, tint, Tina evaluate toe t.nteorais by 
Gauss ouadrature metOod (see page 37), Input paraneters are; 
Ap = Roots of Laouerre polvnomial 
H = height factors of liaouerre polvnomial 
AL = ED, ri/Te • 

SIiaRQtli’IdE FiNl'(AP,H.AL,SnM) 

OlME..\!Sin;M AP(in,H(ti) 

StJM=0,0 
DO iu J*lrll 
Zs=AP( J)/Ali , 

(:! 4 al, 25 i‘C£ 4 f|, 0 ) ' 



7 = At->( J) / -'i. 

F>jM i + , / I f 

Sif.^sS ! .i + K” ;T»-U >1 


Sli '.*='■ ,') 

1 = 1, 1l 
Z = AMf J )/Au 
I11=Z/ C?.+ 1 .0 j 
n2=lO.')4. (Zfi 
Obsiif J)<ni t=A 

PKf.iPri 



10 


50 


70 


^ , 7 J = , j 

T , Tj=u.') 

>- p-n i , Tj 

^ ^ C T , T ) = ' ■ , ) 

«“( 1 i = i 1 .0 

'■’051 = /, 

Psr 

ttir n ( 1 . !-( 3 .u/(D>^P) ) )4|,;f n 

'\j X'» T-1 ,'J 

i-"= I 

' r ) = /. , 0 -S' I-- S' p 
M T j^LCfll/HOfvuU 
P.l=^.7.f « + F( 1 ) 

M =!::r/Ru04,i» 

T = i,!-. 

^ 0 ™ C ^ S '• / 1 C 1. 1 ) / ][* ^' 

y-'C f i ) = t . fGt r ) »(-:kp( -aij ) / f , 51 ) 

r'j v' = 4, i ’■^pOoK-i u/(. rr;.v + *1 .S) 

rp = l Ai - vTi > )/'lF. 

•iPf 1 l=r;t p *CDf,t(.;ypf iPj/2.u 

'ipf I )='!pf 1 . 1 ,op>10 

l)='ip( I, ), 

P.M = R.n*l. ! 4*i. J 4M.R(i ^♦4.R0 7/7 7.i24 

H.IM stilil*! ,09677*1 ,l)qb77 + l ,v)F + Om 

0X1 =6 4.()/( 7.0*SwRrf i,(})*4,0»ATApf 1 .vJl ) 

OiJ 40 1=1, M • 


185 


no ys 
P = T 
0= T 

T F r 1 , t; 7 


r = T , ,M 


1) G1 TO 95 


95 

90 


C., 


56 

55 


50 


470 

460 


490 

480 


AX1 = (l**0)*'<'2/tP*P-0*Q) 

(■;PO(I , 1 J=i ,0-0, I72fi*( ( 1 ,0/AX1 1**0, 7 V«. ( i,7.r,/ (y»1;) ) J Kaxi - I , n 

FPOIT , Tj=tiKl*l1 .0/{2.U*(0*0)'))*CAX1 **i’)iGpn(T,.7l/tf f’*>j)** 

A P 0 I ! , n = 1 0 / P 1 ♦ * 2 * f 1 . 0 / A X 1 1 * P u 7 * F P 0 i T , T J ♦ 7, * 7 Z * , 7 
Apr)(J, J:)=APoh,d) 
roMi'i .NOii: 

CHnTIMUE 

,, For opticaiiv totn case foHowlaa five statements suould be removed 
Apn(2,i j=APQ(2,i)*[)t';Fi 
AP0C7,1 )s:APQ(i,l)>t'DEF2 
AP0U,1 )=APQ(4,l)*0EFi 
AP0t5,1 J=AP0(5,l)*nhF4 
AP0t5,n = AP0(b,l)*nEF5 

no T= 1 ,N 
no 55 J-=I,4 
P=I 

aE=(A( Jl-ACi))/'fF 

TFf I,£0. J) sn Tu^56 ^ ^ : 

KPOCT # J)*4. ?5E-05*P*P*0*O*F.XP(-Anj/C0*0-P*P} 

KPOCT, JJ=KPOCl,vJ)*FPOCT,J}/(ZZ*C'^EP.P**O.F j) 

COM'rTNOE 
CO^IXINIIE 
00 80 .Isl 

no «o ib 1 ,j 

KPatJ,T)»KP0fI,O)*Np(l)/Mw(J) 

DO 460 1»U!H 
Sii^sQ.O 

no 470 J«i,4 

TFCl.EO.J) SO TO 470 
SUMsSysi + KPQd^J) 

CDMTTNge, „ , 

KPOSOMCDaSUM 

no 4 Bo 

§8’'?§3“a.. 

.sffMwso'af 

iip‘ 





p,,r> 


1 0 
C 

1 1 


1 1 


^‘.‘'2*’ j‘':(-'AC c-ilrul rites optical escanef tact^rs 
Taraneteraare: 

O S '' J P ,;j r> j- Q f 1 g y g j 

r - "as t s'^oofaLuro 

c. 1 e 0 1 »■ 0 0 .-i e n 5 1 t V 

'i-uve Ipnjt ii 
s'^oUtator strength 

oooulation oeasity 

i)KKACCt>,T,.'.C,wt, 

. 1 'i ‘a 

.0)/tSgOT( r) t.Wb) 

,'"•'‘3^ f I’p 1 tf 1 , ;+(Pi* t ALP+ALD ) J . 

= 4(.rJ/(iRAl 
= i.it-agA2 

.1 j e:-J^tXiJ*Xh*il3*vir.*Rf<A/fSOHT('^) 1 

. ' + ( tn/c^.u+'i'ott^) ) 

* ul / f KU + si}rtT(Pi ♦ inj 1 

'-’IJ 

M.P* Pjtii^u.O 

r'J,V n . 0+ i 0 U . a* ( SuPT ( AT P 1 >► tu + Af.P^.SORT ( TQ) ) 1 

.o+pi^AT.p^*^) 1 ) 

1 i)5/ J C P 1 ^ f 1 . 0+PT-» Alp t *>! 1 ) 

Sv»^MfPT*T[)6) 
ro7^ rj4!itTu2 

rai^/(i^y+'lo*SgRT(pi>)‘ALnuri .u+iOj) j 

o'- , 

!'a2 


M, - 
T,S = 
V,, = 

n 

Pi'.Al 
p 1 = J 
Xi- = u 

PiTA 1 
Pf'A .i 
PKAi 
R K A - 
r,i=i 
T..li = 
'n,:i? = 
TP t 
TlI^ = 
■f'U-T = 
TgS = 

IMPS 
^ 1)7 = 
T|)H = 
Tjt) = 

npr = 

r,\ 'i r 

nt*r = 


TplUPK.G T.!, ) DPF=1.0 

PtlTURiM 

FnH 


T'lnijt 
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(. _ t?iK.pn trom TB'' Su'^routjne p^c“'-l'’s» d .•.v?!,'*-:; 

‘'rsL^jr^ei iitterent i^l Ha.iiminj's ;n" o Hi <-• r, n / cor'^''.' 

'‘2*7 tlis -letaou is not self stHrr.iuo, t''S st.-irtlo’ 

--JlC'UHtP.I DV H Sfjectnl RijPqe-'^'lttH "ietriO-i. '’’T® lonoor. 

ii'erera are \s follows: ■ 

ill I = iiO'*ier hoiini nt the tnt^^rval 
rr/l = ^poer Hoiin.i of the interval 

ilfi = Jniitiai increment of tne InOeoenu'' it vari-^ulo 

I' i) = 'Jpoer error bound 

fonont vector ot initirii vaiiies, ueter on erriv v ro".CT.i''>, 
foe values of denendent varinoles ruf o la te j al i or er i ore 
^ oa1 lit X. 

V = Vector of error Heionts, f.ater on roncelns the voines 

of ier tvdtlves 

o = iiiitner nt eouatlons in tne svste.ft 

i' s An external suDroutlne wiilcb rofripMt.es the matrix ot 

dimension hum \ .'iDld at -jlven values nt inn^o^oHoni. varlooi 
* An external subroutine to ro'runt'* tfie vector oo toe rin,i»: 
rioot hand sloe of the sYste"! at ilveii values ot in .le j«.sOt>n t 
vnrlaoLe, 

P = /\n ex:teraai subroutine '•riich orints toe caj.ci)aten les 
dt nlven values of Indenenaent varianle. 


SUOci'HJTiMii HPCblPK*'ir,r ,uEKY,MUT«, lUu5',AtCr,FCl’, Jifi^,AUX,Al 

hi iKMdTaM Pu’vf(n,Yiumf^HY(i 
dbi=t 

•JU'? stands for the number of steps after »*!birn tne results are 

desired to or orlnted 

••IU7525 

'10 n lyo 

rhhu hFCT(x,a) 

CA'.U FCrC ft/UFRY) 

OJ i 4 = 1 ,.1bM 

r,b=d-t'jbid 

d3=ii.i' 

no i 0 = 1 fdOiM 

r j (j “ u fi t d 0 1 d 
HS»dS + Mr.L)#Y(L) 
neOf (NUHS+DEHYtM) 

r;u i'lCl05,2a2,204,206,ll5,125,125,5UR,3i?#^^Y»3i'l.tiidJ . 

n=i 

Idbt’sO 

XsPKdXCi) 

t-fsOKdTfS) 

PRdi’C5)ei0,0 

no i0| |ai rHOlH ■ 

AUK(16#I)=0, 

AUX(15,I)»D|RYf I) 

TFfiUlPa5lC2)-X))i0J,lO2,iO4 
THf,.r»l2 ■ 

qo fo ,ia4 

IgSSSi’ 

da 10 4 .. f:,:* 

aAr*r'S!»i#i)i.!f >t)ER>tOHJ<F,HBI«,Pfl«T,NLt .ML?) 

ijKrioi 



i88 


^ > 1 ^ ^ 

•* , ■ 

1 " il 'I ,x 

c , 1 

1 Z “5 ,« 

: « 1 

J f ■( 



r " V i 

»* ^ 



" . i 

•1 ) 



/ z f -t 

1 



T *' 

S 



r 

1 I 







'1 1 1 

^ 'i 

l-i ^ . 

’) 1 ' ' 

1,]X t 


\-i(A 

) 

Attii 

r 

} K 

Y(TJ 

j ,-v 4j 1 

“ i 



■rt ' |i 

^ , j 1 

1 1 



'> M M* 


%j 


i 

1 i 

l^l , w 

or? 

t ^ 

5,}^ 

bft A,, 

YtlS 


,i)MbS(Y(T>-Ailx( j,t) j 


I* ( i) T - 

I f ‘ f i 'I u t i 
! 

t I 

’lii ll tal: 

K s X t '^”1 


1'’ 

) ^ ^ ^ r i '7 'J ^ 1 M 


r i> "1 1 = 4 ' 
lo n 2 0') 

■!=) 

ifsXfO 

r6Y2=; 

":tj n J 

XsPRMlfl) 

Ou i 2t) i = JL fijOlf'i 
.,, ’'ii'f t \1 ; i '^nhPYCl) 

J “‘Vj'ji i t*,uy fB.t-Yip , T J 

f^XClO,Tj+u,04j6of''o7 + nf';P)(fi ) j 

y s X f H 

'J = Mf 1 . , 


FiiW2 = l2 

'^CJ n 1 

'U.l sMul f 1 

rFfPRMTC?^ )10Ki?9a07 
Tie C.’«-4ll lo ,300, JOO 

00 in l=*,NniM 

^UXCM,|)=vf£) 

\dX( «+7 . tIzDERYCI) 
rrf.^-3U.n, 134,300 . 

00 I 33 i = 3 ,!'ir>XY 

Y( TJ=.tv!IUl,l > + 0.3 333333 ♦H^fAUXCB, n + OEr.f + A'fXr 10, T>) 

? U 1 J 1, 2 / 

oij ns i=i,NPiM 

OE{,rsAOAC9,l) + AaX(.10,l) 

OKnaOFbT + DFra' + DELT 

y C I > sftoxci ,n+o, oYS + M+iAUxtB, i) + nEbr+Anx( it , T) ) 

^0 n 127 



nu 201 Xai,NDIM 
XsH+AUXCMt7,I> 

AUVC5,r)aX 

Y(T)*A!tXCN,l>+0,4»X 

xsz+on#'-! 

TSW2»a^ 


s 1 AU3r ( S I J ) f 3 * S s 0 *3 6 5 ♦ a U X ( 6 , U ♦ 3 . *1 3 ' 2 'fi to 



rr-’f^'i-') j '■. ,.i'/, ii4 
nj ,=/,7 
■'^S'Ml=:i,.')r' 

V. M M - I , r ^ -Srt'lA ( v , O 

+ f J=‘VJ A (,»♦'/, D . 

■'S ! t 1 

"<1 :} '*:> I =l ,^*'lj, .! 

^uv I'!-! , r , = y,u 1 

■'<j V i, , r .) Y( T ) 

/ = ■/ ♦• I 

tor.^PsT;., !. 

'''',' .4 1/ l = i,f'iniv. 

ortrj. i- ^ nx C , 1 - t , I U 1 . 3.n .Hi + ti<= t ^ iiX t " + ^ , T ) + aha f -,4 o , X 1 - ii n ■• 
1 ai.k i‘;4 X, f J+Hl'xf I) J 

'M '■/ =I)'‘’U '-'ySoI HR 4.x tjxxl n 

'•1'' 1 1 o, n=Df,f,i 
T b ■.' (i - d 

ID n 1 

'■'D jij i. = ij,-.7ni'-s 

onKI'.i =1.' . l AuXi'i-l , t l-At'X f l l f l , ♦X* ( 11 4 X IT 

1 A ij X r-i 4 *4 , 1 j-A"xr,.+b,i ) j 1 ■ 

X( r r4X,.r/-*j6.)l/*Ai‘XriA,T J 

'll) i ! 0 isi 

HDl, I SbR u'l’t MJX 1 1 ^ , 1 ) + AtiSl A ux (1 o , u J 
!>■ f A ) ) 3 1 1 , ^4;4 , 3 X i 

'D,' I'D 1 


189 


r A . T j 4 


irx *'1,1 =’•'1^ +1 

" A [,U J'l i'-n 1^ , if , r)t:f5 i , ,1 HLP' , '<L. T ... , pP u 1 , b? 7 

T4 C .(bl J NM = 1 

f P f pI?..(T(«i) 1 

3H T). flMb«'-ni 

314 RtTuRiii 

315 Tf rH*(X-PKf-xC2)n.31P,il4,ilA 

3l^i Tp(ABS/X-PRnT(.2j)-0.14XbStH}li14, 317,11/ 

317 TFfL)!f:bT-0.02*PRf‘!T(4j)3lH,31y,30l 

31R n (iHb«’l301 ,301 ,310 

310 TF^^-?)301 ,370,320 

320 TFflBl’RP-4)3')l ,321,321 

321 T.vnosiSTKP/X , 

TF(lSl'EP-lMon-TMnuljfil ,32 2,301 

322 OsH+H 

Trtr.F = lHi,P-! 

Ts'»’fc:P=o 

no 323 i=i,waiM 
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